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SUMMARY 

The objective of the present study was to experimentally invest- 
igate the near wall flow parameters in the blade end-wall corner re- 
gion. The blade end-wall comer region was simulated by mounting an 
airfoil section (NACA 65-015 base profile) symmetric blades on both 
sides of the flat plate with semi-circular leading edge. The initial 7 cm 
from the leading edge of the flat plate was roughened by gluing No. 4 
floor sanding paper to artificially increase the boundary layer thickness 
on the flat plate. The initial flow conditions of the boundary layer up- 
stream of the corner region are expected to dictate the behavior of 
flow inside the corner region. Therefore, an experimental investiga- 
tion was extended to study the combined effect of initial roughness 
and increased level of free stream turbulence on the development of a 
two-dimensional turbulent boundary layer in the absence of the blade. 

The measurement techniques employed in the present investi- 
gation included, the conventional pitot and pitot-static probes, wall 
taps, the Preston tube, piezoresistive type pressure transducer and the 
normal sensor hot-wire probe. The pitot and pitot-static probes were 
used to obtain mean velocity profile measurements within the bound- 
ary layer. The measurements of mean surface static pressure were ob- 
tained with the surface static tube and the conventional wall tap 
method. The wall shear stress vector measurements were made with a 
specially constructed Preston tube. The flush mounted piezoresistive 
type pressure transducers were employed to measure the wall pres- 
sure fluctuation field. The velocity fluctuation measurements, used in 
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obtaining the wall pressure-velocity correlation data, were made with 
the normal single sensor hot-wire probe. 

The changes in the initial flow conditions significantly influ- 
enced the development of a two-dimensional turbulent boundary layer. 
The increase in free stream turbulence level, in the presence of initial 
roughness, had an opposite effect on the variation of different mean 
flow characteristics relative to the case when initial roughness was ab- 
sent. 

At different streamwise stations, in the blade end-wall corner 
region, the mean values of surface static pressure varied more on the 
end-wall surface than on the blade surface. The observed large changes 
in the values of mean surface static pressure on the end-wall surface in 
the comer region were mainly caused by the changes in the curvature 
of the streamlines. The magnitude of the wall shear stress in the blade 
end-wall corner region increased significantly in the close vicinity of 
the comer line. The maximum value of the wall shear stress and its 
location from the comer line, on both the surfaces forming the comer 
region, were observed to change along the comer. These observed ch- 
anges in the maximum values of the wall shear stress and its location 
from the corner line could be associated with the stretching and att- 
enuation of the horseshoe vortex. The wall shear stress vectors in the 
blade end-wall comer region were observed to be more skewed on the 
end-wall surface as compared to that on the blade surface. The differ- 
ences in the wall shear stress directions obtained with the Preston 
tube and the flow visualization method were within the range in which 
the Preston tube was found to be insensitive to the yaw angle. 
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The values of wall pressure fluctuation intensity in the blade end 
-wall comer region were found to be influenced by the changes in the 
strength of the horseshoe vortex. The spectral content of the wall pre- 
ssure fluctuation increased, at any stream wise station, as the blade sur- 
face was approached. The correlation of the wall pressure fluctuation 
with the velocity fluctuation indicated, in general, higher values of cor- 
relation coefficient in the inner region as compared to the outer re- 
gion of the shear layer. Also, the values of wall pressure-velocity cor- 
relation coefficient in the blade end-wall comer region decreased and 
increased in the presence of favorable and adverse pressure gradients, 
respectively, in the streamwise direction. 

Based on the measurements of various flow parameters, it be- 
came quite evident that in the upstream region of the blade leading 
edge the effect of the presence of the blade diminished in a distance 
of approximately 8 and 6 times the blade leading edge diameter in the 
streamwise and lateral directions, respectively. 
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CHAPTER I 

INTRODUCTION 


1 . 1 Statement of the Problem 

The comer flows, generally turbulent and three-dimensional in 
nature, are encountered in many engineering applications. Flow along 
a streamwise comer in a non-circular duct, at the junction of wings 
and tail planes with the fuselage, in the root and tip regions of the 
turbomachinery stator and rotor blades, at the root of flns with the 
heat exchanger surface, and at the junction of a bridge pier in a river 
bed are a few examples of the typical corner flows. The flow at the 
junction of stabilizing fins and a submarine's hull is another practical 
problem of interest where the noise caused due to two body juncture 
is of great concern. As the turbulent wake produced by the juncture 
flow is intercepted by propeller blades a pulse of sound is generated. 
This noise can lead to early detection of the submarine by enemy’s 
Sonar. In turbomachinery some of these types of flows are termed 
blade end- wall flows. 

An understanding of the blade end-wall boundary layer in turbo- 
machinery is particularly important because the total pressure loss as- 
sociated with the end-wall flow is much larger than the blade profile 
loss. The available literature shows that the secondary flow losses form 
approximately half the total losses in axial turbines. In the blade end- 
wall corner region, mainly two types of secondary flows have been 
found to exist. These are called secondary flows of the "first kind" and 
the "second kind". The secondary flow of the "first kind", also termed 
"skew induced secondary flow", leads to the formation of a horseshoe 
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vortex which wraps around the blade. The secondary flow of the 
"second kind" is generated due to the lateral gradients of the Reynolds 
stresses in an anisotropic and inhomogeneous turbulent flow and is 
also called "stress induced secondary flow". The "skew induced sec- 
ondary flow" can occur both in laminar and turbulent flows, whereas, 
the "stress induced secondary flow" occurs only in turbulent flows. 
Moreover, the "stress induced secondary flows" are expected to domi- 
nate in the near wall region. A schematic of the flow in a blade end- 
wall comer region is shown in Fig. 1.1. Numerous experimental stud- 
ies have been reported to understand development of various types of 
secondary flows existing in axial flow turbomachinesll-24]. Based on 
the informative review, Sieverding(25] observed that the experimental 
research on the secondary flows has resulted in a fairly detailed de- 
scription of the complex end-wall flow pattern in a turbine blade pas- 
sage. However, it was remarked that the information concerning the 
turbulent structure of the end-wall flow and the mixing processes on 
both the blade suction side and downstream of the cascade are not 
available. 

To examine the influence of the horseshoe vortex on the end- 
wall losses a simplified blade end-wall comer flow has been studied by 
many investigators[26-47J. The horseshoe vortex effect was investi- 
gated either by mounting a circular cylinder or a constant thickness 
body with a semi-elliptic leading edge on the flat plate. Also investiga- 
tions have been made of a blade with an airfoil section (representing 
one component of a cascade) mounted on the flat plate as in the pre- 
sent study. With this arrangement streamwise curvature effects could 
also be considered. However, most of the previous investigators con- 
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Fig. 1.1 Schematic of Flow In the Blade End-Wall Corner Region 



centrated their efforts to understand the evolution of secondary flows 
and their interaction with the velocity field only. Very little emphasis 
was given to the understanding of the wall shear stress, fluctuating 
wall pressure field, and their interaction with the fluctuating velocity 
field in three-dimensional flows. The interpretation of these quantities 
requires the understanding of turbulence phenomena near the wall 
and comers. 

An interesting fact on the understanding of the near wall turbu- 
lence structure was revealed by Townsend [48]. He hypothesized that 
the turbulent motion in the inner region of a turbulent boundary layer 
can be divided into an active and an inactive part. The active part is 
responsible for the shear stress (x) and scales with the local flow pa- 
rameters (x and z), where z is the normal distance from the wall. The 
inactive or passive and effectively irrotational part is of large scale and 
is associated with the pressure fluctuations of the outer region. 
Bradshaw[49] supported Townsend's hypothesis and noted that the 
only noticeable effect of the inactive motion is an increased dissipation 
of the turbulent kinetic energy into heat in the viscous sublayer. This 
turbulent kinetic energy is supplied from the outer region to the re- 
gion near the wall through the diffusion process. This implies that the 
outer and inner regions of a turbulent boundary layer are interdepen- 
dent. 

The interdependence of the inner and the outer region of the 
turbulent boundary layer was also reported by Thomas[50]. It was 
noted that the large scale structure influences the sweep and burst 
mechanisms (small scale motions) in the wall region. Thomas also 
concluded that the pressure field of the large scale structure did not 
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directly influence the generation of small scale events. The reason 
suggested was that the pressure field of the large eddy did not have 
sufficient intensity to change the velocity profile. 

It can easily be seen by examining the equation for the pressure 
fluctuation[51] that the pressure fluctuation at a point is influenced by 
the entire flow field. Willmarth[52] suggested that the pressure at a 
given point will not be highly correlated with the velocity fluctuations 
at any one neighboring point. It is not difficult to observe from the 
fluctuating pressure equation that the fluctuating values of the wall 
shear stress and the surface pressure are indirectly correlated. Both 
play an important role in the understanding of the near wall turbu- 
lence structure. The wall shear stress is also used as the velocity scale 
in the inner region of a turbulent boundary layer. In a three-dimen- 
sional flow an accurate determination of the wall shear stress requires 
the knowledge of both its magnitude and direction. The wall shear 
stress vector variation in the presence of the two types of secondary 
flows, mentioned earlier, needs to be examined in detail. 

For a two-dimensional turbulent flow the available literature (see 
for detailed review reference 52) shows the existence of a relation be- 
tween the RMS value of the wall pressure fluctuation (p' w ) and the 
time averaged value of the wall shear stress (t w ) in the form 

P’w = C ^w 

where C is a constant of proportionality. The value of C is found to be 
weakly dependent on the Reynolds number and a strong function of 
the pressure transducer diameter (attenuation due to spatial resolu- 
tion) [521. The average value of constant C reported in the literature is 
3 for flows with Mach number less than one[52]. However, very little 
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information is available about the fluctuating wall pressure field in the 
three-dimensional flows. 

The correlation between the wall pressure fluctuation and the 
fluctuating velocity field is also important from the point of view of 
gaining insight about the near wall turbulence structure. The knowl- 
edge of these correlations can be used in identifying the regions of the 
shear layer which contribute to the wall pressure fluctuations. The 
available measurements[49, 53-55), in general, are for two-dimen- 
sional flows. In a two-dimensional turbulent boundary layer the mech- 
anism leading to the production of intense pressure and velocity fluc- 
tuations due to stretching of the vorticity produced by the viscous 
stresses in the sublayer was suggested[55). Bradshaw[49] reported that 
the correlations of wall pressure fluctuation with the velocity fluctua- 
tions are also influenced by the presence of the pressure gradient. The 
present measurements were intended to gather such data in the blade 
end-wall corner region. 

There is another important effect which needs to be examined. 
This is the combined effect of initial roughness and increased level of 
free stream turbulence on the development of a turbulent boundary 
layer. In most of the turbulent boundary layer experimental work some 
transition promoting device is used to trigger the transition and also 
artificially increase the boundary layer thickness. The two most com- 
monly used transition promoting devices are a thin wire and dis- 
tributed roughness strip. It has been shown that a boundary layer de- 
veloping downstream of such devices reaches a equilibrium state if 
compared at the same momentum thickness Reynolds number[56). 
Also the streamwise distance required to achieve the equilibrium state 
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differ for different devices. However, very little information is available 
on whether similar results are achievable if the free stream turbulence 
level is increased, especially when the distributed roughness is used as 
the transition promoting device. 

The above discussion clearly shows that detailed measurements 
of the wall pressure fluctuation intensity and the mean wall shear 
stress vector in the blade end- wall comer region are required. Also 
required is the information on wall pressure-velocity correlations in 
three-dimensional turbulent flows in general, and in a blade end-wall 
comer flow in particular. It is hoped, that this data will help in under- 
standing the reasons leading to the interference losses, flow induced 
noise, and structural vibration problems in the field of turbomachines. 

1.2 Objectives of the Present Investigation 

Based on the discussion in section 1.1, this experimental study 
was undertaken to gain insight about the near wall turbulence struc- 
ture in the blade end-wall comer region. The blade end-wall corner 
region was constructed by mounting a symmetric airfoil (NACA 65-015 
base profile) on a flat plate at a distance of 620 mm from the plate 
leading edge. The specific objectives of the research undertaken were 
as follows : 

(a) To study the effect of initial conditions (distributed rough- 
ness of certain streamwise extent placed downstream of the plate 
leading edge) and increased level of free stream turbulence on the de- 
velopment of a two-dimensional turbulent boundary layer 
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(b) To obtain detailed measurements of time averaged mean 
value of the surface pressure in an upstream region and in the blade 
end-wall comer region 

(c) To obtain detailed measurements of the wall shear stress 
vector in an upstream region and in the blade end-wall comer region 

(d) To measure the wall pressure fluctuation intensity and its 
power spectral density in the blade end-wall comer region 

(e) To measure the correlation of wall pressure fluctuation with 
the streamwise fluctuating velocity component in the blade end-wall 
corner region 

(f) To obtain measurements of steps (b) and (c) in the presence 
of a separation bubble in the blade end-wall corner region. 


CHAPTER n 


LITERATURE REVIEW 


The available literature relevant to the present research work is 
reviewed in four sections. In the first section, the effect of free stream 
turbulence on the development of a two-dimensional turbulent bound- 
ary layer is reviewed. In section 2.2 the present state of knowledge 
about the understanding of secondary flow phenomenon and their as- 
sociated end-wall losses in a blade end- wall comer region is discussed. 
Much emphasis is given to the experimental research work done in 
the last decade. Various available wall shear stress measurement tech- 
niques are briefly described in section 2.3. The literature on the wall 
shear stress measurements in three-dimensional turbulent shear flows 
and difficulties associated with these measurements are also included 
in section 2.3. In section 2.4, literature related to the wall pressure 
fluctuation is reviewed. The available literature on the correlation of 
wall pressure fluctuation with the fluctuating velocity field is also dis- 
cussed in this section. 

2.1 Combined Effect of Initial Roughness and Free Stream 

Turbulence on the Turbulent Boundary Layer Development 

A number of techniques are available to trigger the transition 
and artificially increase the thickness of a turbulent boundary layer. 
Various devices to promote transition and artificially thicken the tur- 
bulent boundary layer such as the circular wire, normal plate or strip, 
slotted or toothed strip, distributed roughness, screen and air jet have 
been used[56,57]. A detailed investigation on the development of a 
turbulent boundary layer downstream of different types of transition 
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devices was reported by Klebanoff and Diehl[56]. Mean velocity profiles 
were found to become independent of flow distortions introduced by 
the transition device at a certain distance from the plate leading edge, 
the extent of which depended on the type of transition device used. 
The comparison of non-dimensional mean velocity profiles showed 

self- similarity for different transition promoting devices at approxi- 
mately the same value of momentum thickness Reynolds number Reg. 

Among the different transition devices tested, the distributed rough- 
ness in the form of a sand paper strip was found to give good results. 

The effects of free stream turbulence on the characteristics of a 
two-dimensional turbulent boundary layer have been studied by many 
investigators (58- 73], Kline, Lisin, and Waitman(58] investigated the 
effect of free stream turbulence on the development of a turbulent 
boundary layer with zero pressure gradient. The boundary layer was 
tripped using 0.009 in. diameter wire located 3/8 in. downstream of 
the plate leading edge. Transition was found to occur within first 3 in. 
of the plate leading edge. The free stream turbulence intensity (Tuoo = 
u'oo/U e ) without the presence of a turbulence generating grid was 
0.5%, where u’oo and U e are the RMS value of the streamwise fluctu- 
ating velocity and the local free stream velocity, respectively. The 
measurements were made at a free stream velocity of 100 ft/s, and 
Tuoo varied between the range 0.005 < Tuoo < 0.2. It was concluded 
that the increase in free stream turbulence increased boundary layer 
thickness, created a fuller mean velocity profile, and increased the 
skin friction coefficient. Also, the law of the wall and the law of the 
wake were found to be influenced by the increased level of free stream 
turbulence intensity. 
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The influence of free stream turbulence and pressure gradient 
on the boundary layer transition was reported by Hall and Gibbings[59]. 
Various available methods for predicting transition including their 
advantages and limitations were reviewed. The mesh size of the grid 
producing free stream turbulence was found to influence the transition 
position. Free stream turbulence intensity and its length scale both 
were reported to influence the boundary layer transition position. An 
empirical relation 

In (R e e tr - R e ei) = ■ 103 TUoo + 6.88 (2.1) 

was proposed to predict the transition Reynolds number. Equation 2.1 

represents a best fit of the available data and the experiment[59], 
where Regj has a value of 190 and Reg^ is the momentum thickness 

Reynolds number corresponding to the transition position. This equa- 
tion is valid only for the zero pressure gradient case. Free stream 
turbulence effects on turbulent boundary layers in the presence of an 
adverse pressure gradient were studied by Evans[61). A distinct de- 
crease in the wake component was observed with the increased level 
of free stream turbulence. This decrease in the wake component was 
attributed to the increased fullness of the outer profile because of in- 
creased mixing with the turbulent free stream. A modification to the 
Coles’s profile was proposed. A reasonable collapse of the velocity-de- 
fect profiles on to a single curve was obtained by multiplying the de- 
fect velocity by a factor of (1 + 5TUoo). Therefore, the velocity-defect 

law in functional form, taking into account the effect of free stream 
turbulence, could be written as 

(1 + 5TUoo)(U e - U)/u* = f [z/5] (2.2) 
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where U e is the local free stream velocity and 8 is the boundary layer 
thickness, and u* is the friction velocity. 

The first detailed investigation of the effects of free stream 
turbulence on a turbulent boundary layer behavior was presented by 
Huffman, Zimmerman and Bennet[62]. Sandpaper of 3.5 grit size hav- 
ing a streamwise extent of 1 in. was used as a boundary layer trip to 
promote the transition. The structural properties of the turbulence 
were observed to be more strongly affected by the changes in the free 
stream turbulence level than the mean boundary layer properties. As 
the free stream turbulence level was increased, the boundary layer 
thickness (8) and skin friction coefficient (Cf) increased, whereas, the 
shape factor (H) decreased. It was concluded that the inner layer 
shear stress producing mechanism is of universal nature and is largely 
unaffected by external conditions. It was also found that the eddy vis- 
cosity increased rapidly with increased level of free stream turbulence 
and its value at Tuoo = 0.05 was approximately an order of magnitude 
greater than the value at TUc* = 0.01. It was also noted that the bound- 
ary layer became highly non-equilibrium in nature mainly because of 
the advection of the free stream turbulent energy and the non-equal- 
ization of production and dissipation. 

Robertson and Holt(63] studied the effect of free stream turbu- 
lence on a turbulent boundary layer in the range of Reg from 400 to 

2000. The free stream turbulence intensity was varied from 1% to 
10%. The skin friction showed a linear variation with Tuoo up to its 
value of approximately 7% and became constant for Tuoo > 7%. The 
maximum increase in Cf value was approximately 20% at Tuoo = 7%. 
The value of 8/0 increased with an increase in Tuoo. The maximum in- 
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crease was approximately 25% at Tuoo » 7%. The value of shape factor 

decreased with an increase in free stream turbulence level. The values 
of different parameters (Viz., Cf, 8/0, and H) at approximately zero free 

stream turbulence level, i.e., in the absence of grid generated turbu- 
lence were evaluated using well established empirical relations 


Cfo = (4.4 + 3.8 logio R e 0 ) 

(2.3) 

(8/0)0 = 8.22 + 0.475 logio R e 0 

(2.4) 

HQ = 1.09 + 0.40 [0.72 logio R e 0 ' 10 1 _1 

(2.5) 


where the subscript 'O' with the different parameters refers to their 
values at zero free stream turbulence. 

Chamay, Comte-Bellot, and Mathieu[64] noted that both the in- 
tensity and the length scale of free stream turbulence influenced the 
characteristics of a turbulent boundary layer. Compared to a turbulent 
boundary layer evolving in an undisturbed flow (i.e., in the absence of 
the free stream turbulence), the influence of increasing the free 
stream turbulence was to increase the growth rate of the boundary 
layer and the wall shear stress. It was also found that the law of the 
wake did not follow the relation applicable to the undisturbed flow 
case. A modification to the velocity-defect law was proposed in which 

(U e - U)/(u* - cqu’oo) = f Iz/8] (2.6) 

where the quantity cqu'oo takes into account the rearrangement due to 

the external turbulence. The measurement of turbulence quantities 
showed that the turbulent shear stress and turbulent kinetic energy 
production increased with an increase in the free stream turbulence 
level. 
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Green[65] analyzed the existing results of the previous studies 
and proposed some empirical relations. The influence of free stream 
turbulence intensity on the fractional change in Of and H is given 

below. 

At constant Reg 

ACf/ CfO= 4.8 Tuoo, (2.7) 

AH/Ho = - [2.4 - 0.25 H] Tuoo, (2.8) 

and at constant H 

ARe 0 /(Re 0 )o = - [0.27/Cf - 0.97(2/C f ) 0 - 5 ) Tuoo (2.9) 

An analytical method to predict the effect of free stream turbu- 
lence on the turbulent boundary layer development was presented by 
McDonald and Kreskovsky[66]. By allowing the entrainment of free 
stream turbulence into the turbulent boundary layer and performing an 
overall turbulent energy balance, it was found that satisfactory quanti- 
tative predictions of the free stream turbulence effect on the turbulent 
boundary layer behavior could be made. Furthermore, a 30% increase 
in skin friction with a 5% increase in free stream turbulence was pre- 
dicted. This result compared well with the experimental data. 

The influence of wind-tunnel turbulence on the development of 
a turbulent boundary layer was investigated by Meier and Kreplin[67]. 
The results were obtained for a turbulence intensity in the range 
0.06% < Tuoo < i%. Based on the experimental results, the following 
empirical relations for the change in skin friction were proposed: 

ACf/CfO = 190 Tuoo 2 Tuoo < 1% (2.10) 


I 
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(2.11) 


ACf/CfQ = 3.8 Tuoo - 0.019 Tuoo > 1% 

The Influence of free stream turbulence in the range 0.3%-6% on 
a turbulent boundary layer development at high Reynolds number (Reg 

= 10^) was studied by Raghunathan and McAdam[68]. The influence of 
free stream turbulence on integral parameter 8/0 was found to be non- 
linear and an empirical relation was proposed: 

(8/0)/(8/0) O = 1 + 500 Tuoo 2 (2.12) 

which fitted the available and their experimental data well for 0 < Tu,» 

< 0.04. It was also pointed out that the results of Robertson and 
Holt[63] showed a linear variation of (8/0)/(8/0)o with Tu,*, which could 

be represented by 

(5/0)/(5/0)o =1+5 TUoo (2.13) 

However, this linear relationship appeared to be accurate only up to 
Tuoo = 0.01 beyond which it underestimates the effect of Tu oa . The in 
fluence of TUc* on change in shape factor was found to be well repre- 
sented by an empirical relation 

AH/Ho = - 2 Tuoo (2.14) 

The influence of free stream turbulence on the wake strength 
parameter (11) was found to be well represented by an empirical 
relation 

An/rio = - 0.9 (u'oo/u*) (2.15) 

The change in skin friction coefficient due to an increased level of free 
stream turbulence was found to be well represented by an empirical 
relation in the range 0 < Tu^ < 0.04 
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ACf/CfD= 12.6TuJ- 4 (2.16) 

The measurements of Raghunathan and McAdam also showed that the 
streamwise turbulence intensity (u’) within the boundary layer was 
more sensitive to free stream turbulence level especially in the outer 
region, whereas, the streamwise turbulence intensity did not change 
in the inner region. It was also observed that the wall pressure 
fluctuation intensity increased with the increase in free stream turbu- 
lence level. For example, the value of p' w /c[°o changed from 0.65% to 
2.05% when the value of TUoo increased from 0.0034 to 0.063, where 
p' w and q«M are the RMS value of the wall pressure fluctuation and free 
stream dynamic head, respectively. 

Some investigators[for example ref. 64, 65, 68] reported that 
both the free stream turbulence intensity and its length scale influ- 
enced the turbulent boundary layer characteristics. However, no 
systematic investigation on the effect of free stream turbulence length 
scale was attempted. 

Recently, a detailed experimental study of the effect of typical 
free stream turbulence length scale ratio (i.e., the ratio of the free 
stream turbulence dissipation length scale, Le. to the boundary layer 
thickness) on the turbulent boundary layer characteristics was re- 
ported by Hancock and Bradshaw[69]. In general, two turbulence 
length scales are commonly used in the literature. These length scales 
are the integral length scale and the dissipation length scale of the 
free stream turbulence. Hancock and Bradshaw preferred to use 
turbulence dissipation length scale which can be easily calculated once 
the decay characteristics of the free stream turbulence downstream of 


16 


a turbulence generating grid Eire known. The turbulent kinetic energy 
equation applied to free stream flow reduces to 

Uoo[d(u’oo 2 )/dx] = - u'oo 3 /Le (2.17) 

and the variation of grid generated free stream turbulence can be 
written as 

[u’2/U 2 loo-°-8 = Cl (x/M - C2) (2.18) 

where M is the mesh size of the grid, x is the stream wise distance 
from the turbulence generating grid, and Cl and C2 are constants 

which can be determined experimentally for any grid. Using Eqns. 
2. 17 and 2. 18 one can easily calculate the value of Le. A wider range of 
turbulence length scale ratio, 0.7 < Le/8 < 5.0, was covered in their 
study. It was found that the effect of free stream turbulence on a 
turbulent boundary layer increased with the increase in free stream 
turbulence intensity. On the other hand, the effect of free stream 
turbulence decreased with the increase in its length scale ratio. It was 
also reported that the skin friction and other boundary layer parame- 
ters have a nonlinear dependence on the free stream turbulence in- 
tensity. The change in skin friction at constant Reg with the change in 

free stream turbulence intensity could be made to fall on a single curve 
if the effect of length scale ratio was also included. The parameter for 
abscissa for plotting change in Cf selected was TUoo/(Le/5 + 2.0), 

where the denominator was selected by trial and error procedure. 

More recently, Castro[70J studied the effect of both the free 
stream turbulence intensity and the length scale ratio on a low 
Reynolds number turbulent boundary layer. This investigation was 
made at a momentum thickness Reynolds number in the range 500 < 
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Reg < 2500, TUoo < 0.07, and 0.8 < Le/8 < 2.9. Direct Reynolds number 

effects were observed to be significant and it was shown that such ef- 
fects diminished as the free stream turbulence level was increased. A 
modification to the empirical correlation, suggested by Hancock and 

Bradshaw[69], relating the fractional change in skin friction coefficient 
at constant Re 0 , to a free stream turbulence parameter containing a 

dependence on both the turbulence intensity and the length scale was 
proposed. At low Reynolds number in the range 600 < Re 0 < 2000, 

the correlation between fractional change in skin friction and modi- 
fied free stream turbulence parameter was written as 

(ACf/CfO)l 1 + (10/F2) exp (- Re 0 /4OO)1 = f [F] (2.19) 

where F = 100 Tuoo/(0.5 Le/8 + 2.5), is a modified free stream tur- 
bulence parameter. 

Effects of free stream turbulence and pressure gradient on mean 
velocity profiles and heat transfer were investigated by Junkhan and 
Serovy(71]. The measurements, obtained at low local Reynolds number 
(4x10* < Re x < 4x10^), indicated no influence on the heat transfer 

with the increase in free stream turbulence intensity. Simonich and 
BradshawI72], however, argued that if Reynold’s analogy were to be 
exact, one would expect the percentage change in turbulent heat 
transfer because of increased free stream turbulence intensity to be 
the same as percentage change in skin friction. Based on the available 

measurements, the change in skin friction due to free stream turbu- 
lence, at constant Re 0 , was written as 

Cf/CfD = 1 + Ax TUoo (2.20) 
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where Aj = 3 and depends strongly on the free stream turbulence 

length scale. For turbulent heat transfer a relation similar to Eqn. 2.20 
was written as 

St/Sto = 1 + Aq Tuco (2.21) 

which gives the Influence of turbulence intensity on heat transfer, 

where St is the Stanton number. The measurements showed the value 
of Afl = 5 suggesting that the Reynolds analogy factor increases slightly 

in the presence of free stream turbulence. Quantitatively, it was re- 
ported that the Stanton number increased by about 5% for each 1% 
increase in free stream turbulence intensity. Also the Stanton number, 
at a given free stream turbulence intensity, decreased as the free 
stream turbulence length scale increased. 

A more detailed investigation of the effect of free stream turbu- 
lence on a turbulent boundary layer was recently reported by Blair[73]. 
This investigation showed increase of approximately 14 and 18 per- 
cent for skin friction and heat transfer, respectively, for a 6% increase 

in free stream turbulence intensity. This study covered Reynolds num- 
bers, Reg, ranging from 1000 to 6000, whereas the free stream tur- 
bulence intensity was varied from 0.25% to 7%. 

The influence of free stream turbulence on the turbulent bound- 
ary layer characteristics was recently reviewed by Pal[74] while study- 
ing the effects of free stream turbulence on the development of wakes 
produced by a thin plate and also the wake-boundary layer interaction 
phenomena. 
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Based on the discussion presented above, the quantitative effect 
of the free stream turbulence on various turbulent boundary layer 
parameters of importance is presented in Table 2.1. 

From the available literature reviewed here it is quite evident 
that very little information about the influence of free stream turbu- 
lence on the development of a turbulent boundary layer developing 
downstream of initially distributed roughness is available. 

2.2 Secondary Flow Effects in End-Wall Comer Flows 

In this section recent advances made in the understanding of 
the development of secondary flows in the blade end-wall corner re- 
gion are described. It is important to mention that no attempt is made 
here to present a detailed review on the corner flows. A detailed re- 
view of existing literature on the corner flows up to mid-seventies was 
presented by Shabaka[36J. Recently, a good description on corner 
flows with the systematic classifications of various types of corner 
flows of practical interest was given by Hazarikal47). However, in this 
section experimental results on the nature of secondary flows and 
end-wall losses obtained by various researchers in the last ten years 
are briefly discussed. 

Various empirical and semi- empirical correlations developed by 
the late fifties to predict secondary losses in cascades and axial-flow 
turbomachines were compiled and discussed by Lakshminarayana and 
Horlock(75]. In this review article, different types of secondary flows 
encountered in axial-flow turbomachines and cascade were described. 
Based on the available experimental results it was remarked that the 
end-wall losses resulted primarily from the separation of the boundary 
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TABLE 2.1. Free Stream Turbulence Influence on Mean Flow 
Boundary Layer Parameters 


Parameter 

Tuoo 

(%) 

% Variation in 
Parameter 

Remarks 

Investigator 

(Year) 

Boundary 

11.0 

14.3 



Kline et al. (1960) 

Layer 

3.9 

50.0 



Huffman et al. 

Thickness 





(1972) 

(5) 







7.0 

45.0 

Reo = 400-2000 

Robertson & Holt 






(1972) 


1.94 

16.7 

Ree = 2113 

Evans (1974) 


3.42 

25.6 

Ree = 1760 

Evans (1974) 


1.0 

10.0 

Ree = 1.0 x 10 4 

Raghunathan & 






McAdam (1983) 

Displacement 

11.0 

4.0 



Kline et al. (1960) 

Thickness 

3.9 

- 9.8 



Huffman et al. 

(6*) 





(1972) 


1.94 

20.7 

Ree = 2113 

Evans (1974) 


3.42 

5.2 

Ree = 1760 

Evans (1974) 

Momentum 

11.0 

10.0 



Kline et al.(1960) 

Thickness 

3.9 

- 1.6 



Huffman et al. 

(0) 





(1972) 


1.94 

25.0 

Ree = 2113 

Evans (1974) 


3.42 

8.0 

Ree = 1760 

Evans (1974) 

8/0 

I 

7.0 

35.0 



Robertson & Holt 






(1972) 


1.94 

- 6.7 

Re e = 2113 

Evans (1974) 


3.42 

16.0 

Ree = 1760 

Evans (1974) 


3.14 

49.0 

Re e = 1.0 x 10 4 

Raghunathan & 






McAdam(1983) 

H = (8*/0) 

11.0 

- 5.5 



Kline et al.(1960) 


3.9 

- 8.5 



Huffman et al. 






(1972) 
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TABLE 2.1 (Continued) 


Parameter 

Tuoo 

% Variation in 

Remarks 

Investigator 


(%) 

Parameter 



(Year) 

H = (8*/6) 

7.0 

- 5.0 

Ree = 400-2000 

Robertson & Holt 






(1972) 


4.1 

- 8.7 

Reported by Pal 

Tsuji & Iida 




(1981) 

(1972) 


1.94 

- 3.5 

Ree = 2113 

Evans (1974) 


3.14 

- 8.0 

Ree = 1.0 x 10 4 

Raghunathan & 






McAdam (1983) 

Wake 

1.94 

- 81.0 

Ree = 2113 

Evans (1974) 

Strength 

3.14 

- 55.0 

Re e = 1.0 x 10 4 

Raghunathan & 

Parameter (n) 





McAdam 

(1983) 

Skin Friction 

3.9 

16.9 



Huffman et al. 

Coefficient 





(1972) 

(Cf) 

4.1 

10.8 

Reported by Pal 

Tsuji & Iida 




(1981) 

(1972) 


1.94 

5.3 

Ree = 2113 

Evans (1974) 


3.42 

10.5 

Ree = 1760 

Evans (1974) 


3.14 

9.9 

Ree = 1.0 x 10 4 

Raghunathan & 






McAdam (1983) 

Ree 

1.94 

25.0 



Evans (1974) 


3.42 

4.0 



Evans (1974) 


3.14 

47.0 



Raghunathan & 
McAdam (1983) 

Mixing 

4.1 

96.0 

Reported by Pal 

Tsuji & Iida 

Length dm) 



(1981) 

(1972) 

Eddy Viscosity 

4.1 

233.0 

Reported by Pal 

Tsuji & Iida 

(e) 



(1981) 

(1972) 
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layer on the comer formed by the end-wall and the blade suction sur- 
face. An interesting explanation was presented for the observed higher 
losses in turbines and relatively smaller losses in compressors due to 
relative motion between the wall and the blade. It was noted, that the 
effect of a moving wall was to impart a rolling motion to the air in the 
vicinity of the blade leading edge. Moreover, the effect of moving wall 
was observed to be favorable in the case of compressor and adverse in 
the case of turbine. The above mentioned factors explained the differ- 
ences in the influence of the relative motion between the wall and the 
blade in turbines and compressors. 

Dunham[76] reviewed different analytical and experimental 
studies on secondary losses in a turbine cascade. Using the available 
results, an empirical relation was proposed which correlated the ef- 
fect of blade loading and incoming boundary layer thickness on sec- 
ondary flow losses. A list of various parameters which influence sec- 
ondary flow losses were provided and their effects, based on the avail- 
able data, were analyzed. Also different mechanisms of the secondary 
flow losses as envisaged by numerous investigators were discussed. He 
also noted that the secondary flow losses were approximately half of 
the total losses in a turbine. Based on this study he concluded that 
even though the flow pattern associated with secondary losses is un- 
derstood the reasons for the magnitude of these losses are not well 
understood. Therefore, it was suggested that more detailed studies of 
the flow inside the blade passage would be required before accurate 
secondary loss prediction schemes could be developed. 

The suggestion of Dunham did not go in vain as the number of 
detailed experimental studies to understand three-dimensional flows 
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in the blade passage came into existence in the literature in late sev- 
enties. Langston, Nice and Hooper[2J reported a detailed experimental 
study of the three-dimensional flow in a large scale turbine cascade 
with particular emphasis on the end-wall flows. Mean flow and surface 
pressure measurements on the blade and the end-wall, including flow 
visualization studies, were made. The effect of the change in incoming 
boundary layer thickness was also investigated. Three-dimensional 
flow separation was observed to be an important feature of the end- 
wall flows. The position of the saddle point of separation and its asso- 
ciated attachment and separation lines were found Important In 
defining the flow field near the end-wall. The saddle point was ob- 
served to depend on the pressure field in the blade passage and not on 
the change in incoming boundary layer thickness. However, it was 
noted that in the case of symmetric flow (such as cylinder mounted on 
a flat plate) the location of the saddle point depended on the state of 
the incoming boundary layer (laminar or turbulent), the boundary layer 
Reynolds number, and the cylinder diameter[4]. The size of the pas- 
sage vortex was found to increase in the downstream direction en- 
training fluid from the main stream as well as separated fluid from the 
end- wall and the blade suction surface boundary layers. The passage 
vortex was found to significantly alter the pressure field from that of a 
potential flow. It was concluded that any attempt to model analytically 
all or part of the flow field in the cascade passage must in some way 
take account of the three-dimensional separation process. 

Using the results of extensive mean flow measurements and flow 
visualization studies made by Langston, Nice, and Hooper(21, a model 
(Fig. 2.1) for the secondary flows in a cascade passage was proposed by 
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Fig. 2. 1 End -Wall Flow Model in Cascade by Langston(51 


Langston[51. General mechanisms leading to the generation of differ- 
ent secondary flows in the cascade as described by Langston and also 
later observed by many other investigators will now be discussed. Basi- 
cally, there are two mechanisms which influence the generation of 
secondary flows in a turbomachinery passage. First of all, secondary 
flow is produced due to rotation or curvature which causes the flow 
with an initial distribution of normal vorticity to generate a stream wise 
component of vorticity. This phenomenon, it should be noted, also oc- 
curs in curved ducts. Secondly, as the end -wall boundary layer en- 
counters the blade leading edge, a horseshoe shaped vortex is formed. 
A majority of the fluid from the incoming boundary layer gets trapped 
in this vortex which has two legs. These two legs are known as pres- 
sure- and suction-side legs. During this process the normal com- 
ponent of the vorticity associated with the incoming boundary layer is 
transformed into the streamwise component of the vorticity. This 
phenomenon is a dominating feature in simple blade end-wall flows 
(for example wing-body junction). 

The interaction of the pressure-side leg of the horseshoe vortex 
with the crossflow from the end-wall boundary layer formed within the 
cascade and the entrained fluid from the mainstream flow in the cas- 
cade passage leads to the formation of a vortex which is generally 
known as a passage vortex. The suction-side leg of the horseshoe vor- 
tex remains in the suction-end-wall corner region and is known as a 
counter vortex. The counter vortex has been observed to have a sense 
of rotation opposite to the passage vortex and is much smaller in size 
than the passage vortex. These vortices explained above are shown in 
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Fig. 2.1. Langston also noted that the counter vortex may be dissipated 
by viscosity. 

The loss mechanisms and the behavior of secondary flows down- 
stream of a large scale stator and rotor by Hunter[6] and downstream 
of a stator by Binder and Romey[7] were investigated. Hunter found 
that the presence of spanwise pressure gradient in the annular 
geometry enhanced a migration of boundary layer fluid towards the 
hub end-wall via the blade trailing edge and on the suction surface in 
the laminar flow regions. It was also observed that the secondary flows 
in the rotor blade row, mainly at the high turning hub section, were 
considerably more complex than in the case of stators. The compari- 
son of data obtained in a single stage turbine with the cascade data re- 
vealed that the behavior of the end-wall flow in the stators and cascade 
was almost similar. A major difference was the existence of apprecia- 
ble radially directed flow in the case of stator. Binder and Romey 
noted that the mixing losses existed far downstream from the trailing 
edge. They also noted that some vortices existed even far downstream 
of the trailing edge which will affect the flow in the following rotor 
passages. 

The secondary flow field at each downstream plane was found to 
be dominated by a passage vortex which decayed in its strength be- 
cause of mixing occurring in the flow. It was also noted that the in- 
crease in total pressure loss was accompanied by a decrease in the 
secondary kinetic energy of the flow. Furthermore, they found that the 
spanwise component of secondary velocity contributed more than 75% 
to the total decrease in the secondary kinetic energy. Finally, it was 
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remarked that more than one-third of the total secondary losses were 
found to occur downstream of the blade trailing edge. 

Sieverding and Van Den Bosche[8] studied secondary flows in a 
turbine cascade using a colored smoke flow visualization technique. In 
this technique a tungsten wire of 0.25 mm diameter was coated with 
Waxoline dye mixed with oil. The dye evaporated as the wire was 
heated and produced smoke filaments which could be obtained in dif- 
ferent colors. The main conclusion, besides identifying various sec- 
ondary vortices, was that the horseshoe vortex and the passage vortex 
did not exist independently of each other but were observed to be a 
part of the same vortex structure. 

Moore and Ransmayr[10] experimentally investigated the effect 
of the blade leading edge shape on the overall losses in a large scale 
rectilinear turbine cascade. Two leading edge shapes (circular cylinder 
and a wedge) were studied. The influence of leading edge shape on 
different measured parameters viz., surface pressure on the blade and 
end-wall and the total pressure loss and the flow directions down- 
stream of the cascade, was found insignificant. The loss distribution 
measured downstream of the cascade showed regions of high loss near 
the mid-span irrespective of the shape of the blade leading edge. In 
fact Moore and Ransmayr expected suppression of the formation of a 
horseshoe vortex with wedge shape leading edge. Sieverding[25] pre- 
sented an interesting explanation on the leading edge effects observed 
by Moore and Ransmayr. He observed that, in the case of cascade, the 
stagnation streamline is curved and the incidence angle to the sharp 
nose will be different from zero. Therefore, the flow "sees" a sharp 
leading edge as a blunt leading edge. Hence the size of the horseshoe 
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vortex will greatly depend on the incidence angle of the incoming 
flow. Sieverding also observed that the shape and position of the pas- 
sage vortex would not be much influenced by the reduction or 
elimination of the horseshoe vortex. Furthermore, the rotational speed 
of the particles near the vortex center might be affected by the 
strength of the pressure-side leg of the horseshoe vortex. 

Sieverding[25] reviewed the existing results of secondary flow 
research and presented a detailed description of secondary flow 
structures and their effects on end-wall boundary layer characteristics 
and loss growth through a rectilinear turbine blade cascade. He noted 
that the vortex model in a turbine cascade proposed by Kline 
(reference No. 8 of Sieverding) and Langston[5] differed with regard to 
the interpretation of the suction-side leg of the horseshoe vortex. The 
main difference in these two vortex models was as follows: Langston 
suggested that the suction-side leg of the horseshoe vortex continued 
in the suction surface-end-wall comer, whereas Kline suggested that 
the suction-side leg of the horseshoe vortex gradually dissipated as it 
came in contact with the passage vortex. Also, Marchal and 
Sieverding[3] in their flow visualization study found that the suction- 
side leg of the horseshoe vortex was observed in the trailing edge 
plane on the mid-span side of the passage vortex rather than in the 
suction-side end-wall comer. Among the various loss contributions the 
comer losses in the suction surface-end-wall corner was found to be 
more significant than the end-wall losses in the pressure surface-end- 
wall corner. In conclusion, two important remarks were made: (i) a 
clear description of the merging process between pressure-side leg of 
the horseshoe vortex and the passage vortex and also of the interac- 


29 


tion between suction-side leg of horseshoe vortex and the passage 
vortex has yet not emerged, and (ii) the experimental research on 
secondary flows has resulted in a fairly detailed description of the 
complex end-wall flow pattern in turbine blade passage. However, 
more information concerning the turbulent structure of the end-wall 
flow and the mixing process on the both the blade suction-side and 
downstream of the cascade is desired. 

At the same time the Sieverding review became available, an in- 
teresting flow visualization study on the spatial development of 
streamwise vortices in a turbine cascade was reported by Sonoda[12). 
Based on many cross-sectional views normal to the flow, the evolution 
of the horseshoe vortex into streamwise vortices and the generation of 
a new type of streamwise vortex pair on the suction surface near the 
end-wall was described in detail. A brief description on this recent 
study on the secondary flow development in a turbine blade cascade is 
presented here. More detailed description can be obtained in Sonoda's 
original paper. 

At x/C a = 0.3 (section 1 in Fig. 2.2) the pressure-side and suc- 
tion-side legs of the horseshoe vortex were observed to exist as a vor- 
tex pair, where C a is the axial chord length of the blade. As this vortex 
moved downstream, at x/C a = 0.85 (section 10 in Fig. 2.2) it rotated by 
90^ and then was observed to rotate by another 90® at the blade trail- 
ing edge. This rotating motion of the vortex pair was attributed to the 
increase in circulation of the pressure-side leg (Hp) of the horseshoe 
vortex due to the transverse pressure gradient as compared to the 
circulation of the suction-side leg (Hs) of the horseshoe vortex. Also as 
the vortex pair moved downstream it was observed to be divided into 
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two parts at approximately x/C a = 0.54 (section 5 in Fig. 2.2). One part 
corresponded to the "leg" of the vortex pair (Hs and Hp), whereas the 
other part evolved as another vortex pair (termed as III and IV). It was 
noted that the migration of the root of the "leg" of this new vortex pair 
on the suction surface corresponded to Langston’s separation line on 
the suction surface. This study also concluded that the pressure-side 
and suction-side leg of the horseshoe vortex did not exist indepen- 
dently but as a pair of vortices. Furthermore, this vortex pair showed a 
mutual vortical induction motion which in the literature is known as a 
passage vortex. 

Moustapha, Paron and Wade[13] reported experimental mea- 
surements on the secondary losses in a large scale annular cascade of 
highly loaded turbine rotor blades. The blades had a turning angle of 
128.50 and aspect ratio of 0.88. Distribution of the losses and flow an- 
gles revealed the presence of two large vortices which occupied a 
major portion of the trailing edge plane. A large high loss core was ob- 
served in the center of the blade passage and it coincided with the re- 
gion of maximum flow undertuming. 

More recently, Govardhan, Venkatrayulu, and Prithvi Raj [14] in- 
vestigated the effect of the incoming boundary layer thickness on sec- 
ondary losses in large deflection (120®) annular turbine cascade. The 
incoming boundary layer thickness was systematically varied both on 
the hub- and the casing-walls. In general, the loss core near the hub- 
wall was found larger than that near the casing-wall. With the increase 
in inlet boundary layer thickness on the hub-wall, the magnitude of the 
hub loss core increased due to the increased feeding of the inlet 
boundary layer fluid into the loss core by the pressure-side leg of the 
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horseshoe vortex. However, with the increase in inlet boundary layer 
thickness on the casing-wall, the magnitude of the hub loss core was 
again increased, whereas the magnitude of the casing core loss re- 
mained more or less unchanged. In the latter case increased radial 
flow from the tip to the hub of the blade was found to be the cause of 
the differences observed. The measurement of pitchwise mass aver- 
aged flow angles showed that the flow overturned near the walls and 
undertumed in the region of loss cores. The degree of underturning 
increased with the increase in inlet boundary layer thickness. 

Production and development of secondary flows and associated 
losses in stator and rotor blade rectilinear cascade were measured by 
Yamamoto! 16]. In the case of the stator, a significant contribution of 
the secondary flows on the loss production process was observed only 
near the blade suction surface downstream of the cascade throat. 
However, in the rotor blade case larger turning resulted in a much 
stronger contribution of the secondary flows on the loss development 
mechanism. It was also noted that the evolution of overall loss started 
quite early within the blade passage and also the rate of the loss 
growth was much larger in the rotor blade case compared to the stator 
blade case. 

Very recently, the growth of losses, secondary kinetic energy, 
and the streamwise vorticity in a highly turning turbine rotor blade 
cascade was experimentally investigated by Gregory-Smith, Graves, 
and Walshf 18]. The effect of a change in inlet boundary layer thickness 
on various parameters was also considered. A new technique was de- 
veloped to calculate streamwise vorticity component from the mean 
flow velocity data. The new method consisted of fitting a least square 
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bi-cubic spline surface to the data points and then the gradient of the 
surface was taken in the required direction at each data point coordi- 
nate. It was noted, that the negative vorticity associated with the pas- 
sage vortex agreed with the predictions of classical secondary flow 
theory in the early part of the blade passage. However, towards the 
blade passage exit, distortion of the flow by the secondary velocities 
resulted in inaccurate predictions by analytical methods. The positions 
of the loss core, end-wall loss region, and vortex center were observed 
to be alike with the change in inlet boundary layer thickness. The 
main difference caused due to change in inlet boundary layer thick- 
ness was that the above mentioned features became broader and less 
intense with an increase in inlet boundary layer thickness. For exam- 
ple, the maximum values of positive and negative vorticity were found 
to be much lower in the thin boundary layer case compared to the 
thick boundary layer case. This observed change was attributed to the 
high values of vorticity in the thin boundary layer which are con- 
centrated near the end-wall. Moreover, areas of positive vorticity were 
observed to be associated with the transfer of losses into the bulk flow. 

The first measurements of turbulence characteristics related to 
the secondary flows in the blade passage of a turbine cascade were re- 
ported recently by Zunino, Ubaldi, and Satta[19], The secondary flow 
pattern was observed to be dominated by the presence of a passage 
vortex which caused high losses near the suction side. The passage 
vortex was found to play a double role: it was found to transport low 
momentum flow towards the suction side and changed the mean ve- 
locity and turbulent stress distributions. The action of Reynolds 
stresses on the shear flow was then to extract energy from the mean 
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flow and produce losses within and downstream of the cascade. It was 
concluded, that the secondary loss generation in cascades is related to 
the production of turbulent kinetic energy in the vortex flow and its 
viscous dissipation. 

Using the available studies in conjunction with the flow visualiz- 
ation studies by Sharma and Butler! 1 7 j, an interesting interpretation 
of the mechanism of secondary flows and associated secondary losses 
in a turbine cascade evolved. Sharma and Butler agreed with the mec- 
hanism of the formation of the horseshoe vortex as described earlier. 
But, according to them, the streamwise component of the vorticity 
due to the horseshoe vortex is generated as soon as the fluid enters 
the cascade before most of the turning of the mainstream flow has 
been achieved. However, it was noted, that the existing secondary flow 
theories used for the prediction of secondary vorticity relies on the 
assumption that the inlet boundary layer entering the cascade exper- 
iences flow turning in the passage. Therefore, the assumption made in 
the theory is inconsistent with the revelations made by their flow visu- 
alization studies. Moreover, their flow visualization studies also re- 
vealed that the fluid particles closest to the wall in the cascade inlet 
boundary layer did not become part of the horseshoe vortex. These 
fluid particles, in fact, were observed to move towards the suction-side 
of the blade where they climbed the blade surface and exited the blade 
on the top of the passage vortex. 

Moreover, the pressure-side leg of the vortex was observed to be 
immediately influenced by the blade to blade pressure gradient as it 
entered the passage and moved towards the suction-side, meeting the 
surface near the minimum pressure point. Since all the fluid particles 
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from the inlet boundary layer either became part of the horseshoe 
vortex or moved towards the suction side, a new boundary layer 
started at the end-wall downstream of the separation line defined by 
the horseshoe vortex system. As the pressure-side leg of the horse- 
shoe vortex moved across the passage it entrained low momentum 
fluid particles from this newly formed passage end-wall boundary layer 
and it grew to become the passage vortex. This passage vortex was 
noted to be the most dominant feature of the flow. This entrainment of 
the passage end-wall boundary layer fluid significantly affected the de- 
velopment of the passage vortex and is a key mechanism in the gener- 
ation of the secondary flows and end-wall losses. Since this mecha- 
nism is governed by viscosity, inviscid secondary flow theories are not 
expected to provide accurate magnitudes of the secondary vorticity or 
end-wall losses in a cascade. 

Furthermore, in the close vicinity of the leading edge, the suc- 
tion-side leg of the horseshoe vortex was observed to follow the con- 
tour of the airfoil surface. Also, the suction-side leg of the horseshoe 
vortex was observed to remain near the end-wall until the point where 
the passage vortex interacted with the suction-side of the blade. Near 
the minimum pressure point in the cascade passage, the suction-side 
leg of the horseshoe vortex was forced off the end-wall as it interacted 
with the passage vortex. As the flow proceeded downstream, the suc- 
tion-side leg of the vortex moved around the outside edge of the pas- 
sage vortex. This observation, it should be noted, is similar to that ob- 
served by Marchal and Sieverding[3J. The location of this suction-side 
leg of the horseshoe vortex in the cascade exit plane depended on the 
size and vorticity associated with it and the passage vortex. 
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Based on the flow visualization and the discussion presented 
above, a schematic of the flow in turbine cascade as proposed by 
Sharma and Butler is shown in Fig. 2.3. The height (Zte) of the sepa- 
ration line S2 at the trailing edge of the blade approximates the 

diameter of the passage vortex. It was noted that the knowledge of the 
penetration height of this separation line (S 2 ) at the trailing edge on 

the suction surface can be used in estimating both the end-wall losses 
and the secondary flows. Based on the available data it was reported 
that penetration of the separation line towards the mid-height of the 
blade suction surface was found to give good estimates of the sec- 
ondary flows. It was also found from the available literature that the 
region adjacent to the suction surface between the separation line and 
the end-wall was one of high loss and high external heat load associ- 
ated with the passage vortex. The information of separation line and 
its penetration height, especially at the trailing edge, was used to de- 
velop an empirical relation to predict end-wall losses. 

The earlier measurements showed that the net loss generated 
inside the passage was approximately same for both the thick and thin 
inlet boundary layers. Using the foregoing observation, the total loss 
was divided into two parts; namely, the inlet loss contained within the 
cascade inlet boundary layer and the passage loss. The passage loss was 
further divided into the profile loss and the end-wall secondary flow 
losses. An empirical equation was suggested to predict end-wall losses 
and was found to agree within ±10% compared with the available ex- 
perimental data. 

The measurement of mean flow and turbulence characteristics 
inside the hub and annulus wall boundary layer in the trailing edge and 
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near-wake regions of a moderately loaded compressor rotor blade was 
reported by Lakshminarayana and RavindranathI20J. Interaction of the 
wake with the annulus-wall boundary layer, secondary flow, tip leakage 
flow, and the trailing vortex system resulted in slower decay and 
larger width of the wake. The end-wall boundary layers and the sec- 
ondary flows were found to have a significant influence on both the 
decay characteristics and the profile of the wake. 

A detailed investigation of three-dimensional flow in a low speed 
axial compressor stage with an aspect ratio of two was reported by 
Cyrus[22]. The flow in a compressor stage was found to be considerably 
affected by a large area of low total pressure fluid arising near one end- 
wall of the rotor and stator rows. The origin of large high loss regions 
in each blade row, rotor and stator, was found by means of a diffusion 
factor. 

The experimental investigation of the blade and the end-wall 
boundary layers in a subsonic compressor stator was reported by Gallus 
and Hoenen[23]. An increase of the separated region of the hub 
boundary layers and reverse flow at the blade suction side were ob- 
served with the decrease in flow coefficient. 

In an axial compressor excessive diffusion was found to cause 
substantial regions of separated flow near the junctions of the blade 
suction surface and the end-wall[24]. This leads to deterioration in 
compressor performance. Dong, Gallimore, and Hodson[24] observed 
that on the hub of the rotor, where the diffusion factor was 0.56, sig- 
nificant secondary flows existed. However, no suction surface-end-wall 
separation was found. But at both the ends of the stator, suction sur- 
face-end-wall comer separation was found to exist. Also, the sepa- 
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ration was more prevalent on the hub. Surface flow visualization stud- 
ies revealed the existence of large radial movements of the surface 
flow within and outside these regions. It was reported that the intro- 
duction of clearance between the stator blade and the hub resulted in 
a significant reduction in the extent of the suction surface-hub separa- 
tion with an attendant decrease in loss. 

It is quite evident from the discussion presented so far that the 
passage vortex is a dominating feature in cascades and axial turboma- 
chines. Also it is the existence of the horseshoe vortex which leads to 
the formation of a passage vortex. Therefore, it would be of interest to 
first understand the effect of the horseshoe vortex alone on secondary 
flow generation in the end-wall comers. It appears that the above 
mentioned fact led many researchers, in recent years, to study the 
structure of the horseshoe vortex itself. It is often easy to study simple 
fluid flows to explain and predict the behavior of more complex flows. 

Many investigators (for example, references 29, 30, 32) studied 
the cylinder end- wall flows to examine the horseshoe vortex structure. 
Belik[29J reported that the saddle point location depended on the 
state of the incoming boundary layer (laminar or turbulent), the 
boundary layer Reynolds number and the cylinder diameter. Han, Ma, 
and Rapp[30] found heat transfer in the vicinity of the leading edge 
near the end-wall plane higher for a thick inlet boundary layer than for 
a thin boundary layer. 

Ishii and Honami[33] generated the horseshoe vortex effect by 
impinging a confined jet on the plane wall placed normal to the jet 
flow main direction. It was reported, that four vortices existed in the 
detached region which were generated by different processes. The 


40 


detached region could be divided into two parts (region I and II). Re- 
gion I included the horseshoe vortex and the comer vortex, whereas 
region II contained the detachment point and the separation and 
counter vortex. A new model of the flow field associated with the de- 
tached region was suggested. 

More recently, Eckerle, and Langston[32] investigated the hors- 
eshoe vortex formation around a large diameter (D/8 =10) circular cy- 
linder centered between the sidewalls of a wind tunnel. The main aim 
of their study was to obtain measurements in the vortex region, es- 
pecially in the region of vortex formation (l.e., between the saddle 
point and the leading edge of the cylinder). Only one end-wall saddle 
point and a single horseshoe vortex were observed. However, Moore 
and Forlini[31] reported a multi-vortex system implying the existence 
of multiple saddle points. Eckerle and Langston also noted that the 
horseshoe vortex was fully formed between 5® and 25® planes and not 
on the streamwise plane of symmetry (i.e., on the stagnation stream- 
line). Furthermore, the horseshoe vortex was observed to be confined 
to a region less than one boundary layer thickness from the end -wall. 
Also the scale of the separation process at the saddle point was found 
to be much smaller than the scale of incoming boundary layer thick- 
ness. It was also observed that the vortex core moved radially outward 
and slightly away from the end-wall as the vortex moved around the 
cylinder. It was concluded that a key implication of these measure- 
ments was that the saddle point region did not exhibit swirl. 

In other studies the end-wall comer region was generated ei- 
ther by mounting a constant thickness body with semi-elliptic leading 
edge[36-41] or an airfoil section body[42-47) on a flat plate. Shabaka 
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[36] made extensive mean and turbulent flow measurements in a simp- 
lified wing-body junction. The wing had an elliptic nose with constant 
thickness in the downstream direction. It was found that the atten- 
uation of the skew-induced secondary flow by the Reynolds stresses 
was very slow. Also, the skew-induced secondary flow dominated over 
the stress-induced secondary flow even at a distance from the wing 
leading edge equal to 50 times the boundary layer thickness at the 
wing leading edge location. The turbulence structure was observed to 
be far from equilibrium. A logarithmic region was found to exist in the 
mean velocity profiles in the first 10% to 20% of the shear layer 
thickness, except very close to the comer (-20% of the shear layer). 
It was also reported that the eddy viscosity attained a negative value in 
the comer region. The detailed measurements of the surface pressure 
and wall shear stress on both the surfaces in the corner region were 
also obtained. The wall shear stress measurements were made with 
the Preston tube; however, no information about the direction of the 
wall shear stress was obtained. 

Oguz[37] and Kubendran[38] reported mean flow and turbulence 
measurements in a simplified wing-body junction formed by a constant 
thickness body having an elliptical (1.5:1) leading edge mounted nor- 
mal to a flat plate. Kubendran also investigated the influence of the 
leading edge shape on various parameters. The strength of the sec- 
ondary flow vortex in the corner region was found to increase as the 
leading edge of the wing was made more blunt. 

The effect of the wing nose shape on the flow in a wing-body 
junction was also studied by Mehta[39]. The effects of three types of 
wing nose shapes (elliptic, super- elliptic and wedge-elliptic section) 
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on the size, position, and strength of the horseshoe vortex were ex- 
amined. The size and strength of the horseshoe vortex was found to 
increase with nose bluntness. This study revealed that the circulation 
in the horseshoe vortex decreased by approximately 40% as the nose 
shape changed from super-elliptic to a wedge-elliptic. Hawthorne[26] 
also demonstrated the effect of wing leading edge shape on the flow in 
the end-wall corner. He performed a simple experiment where the 
scouring of sand about a strut on the floor of water channel was taken 
as the strength of the vortex. Two leading edge shapes (bicuspid and 
elliptic) were examined. He found that the sand about the strut with 
the bicuspid leading edge was less influenced than that about an ellip- 
tic leading edge. 

Chu, Rios-Chiquete, Sarohia, and Bemstein[40] developed a new 
probe, based on the same principle as the five-hole conrad tube yaw 
meter, to measure the mean velocity vector. This probe was used to 
measure the magnitude and direction of the mean velocity in a wing- 
body junction. The wing-body junction was formed by mounting a wing 
of constant thickness and semi-cylindrical leading edge on a flat plate. 
The measurements were obtained at two axial stations; namely, up- 
stream of the wing leading edge and in the separation region. Yaw an- 
gles greater than 9()0 were observed very close to the wall indicating 
flow reversal ahead of the wing. 

The strut-wall intersection losses were investigated by Barber 
[43]. The strut- wall junction was formed by mounting a symmetric air- 
foil section (65 series) normal to a flat plate. The losses were obtained 
by examining the wake region of the strut- wall junction. The effect of 
the strut thickness and incoming boundary layer thickness on the 
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strut-wall junction was also studied. It was found that the end-wall in- 
tersection losses were significantly lower for a thicker incoming boun- 
dary layer as compared to a thinner incoming boundary layer. It was 
also observed that the flow distortions extended up to one chord 
downstream of the intersection. 

The spatial extent of the horseshoe vortex legs from the wing- 
plate junction at an incidence to the main flow was studied by Rood 
[44]. It was observed that the larger of the two vortex legs produced by 
the wing at incidence dominated the flow downstream of the junction. 
He also noted that the amount of skew in the downstream region of 
the junction depended on the wing geometry and the incoming boun- 
dary layer thickness. It was concluded that the two-point velocity corr- 
elation measurements provided better evaluation of the spatial extents 
of the horseshoe vortex legs on either side of the wing than the mean 
velocity measurements alone. 

Kubendran and Harvey[45] showed the reduction in juncture 
drag with the use of leading-edge fillets. However, they noted a need 
for optimizing the fillet size to achieve overall improvement in the 
flow field. 

The turbulent flow measurements in wing-body junction using 
laser velocimetry were reported by Scheiman and Kubendran[46]. The 
wing-body junction was formed by mounting a wing, having curved 
surface (circular arc) on one side and a flat surface on the other side, 
normal to the flat plate. This construction of the wing was used to 
study the comer flow with and without the streamwise pressure gra- 
dient. The effect of sharp leading edge and comer fillets on the junc- 
tion flow field were also examined. In the case of juncture flow with a 
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sharp leading edge, no evidence of the presence of a horseshoe vortex 
was observed. It was noted, that the use of juncture fillets improved 
the flow characteristics near the wing trailing edge by removing re- 
gions of low momentum fluid. Recently, detailed measurements of 
mean and turbulent quantities using an inclined single sensor hot-wire 
probe in the blade end-wall comer region were obtained by Hazarika 
[47]. The blade end -wall comer region was constructed by mounting a 
blade (NACA 65-015 airfoil section) normal to a flat plate. Flow visua- 
lization using an ink-dot matrix method and oil film technique was 
also carried out. Based on total pressure, static pressure, mean velocity 
components, and the results of flow visualization, a horseshoe vortex 
and the comer eddy were identified upstream of the blade maximum 
thickness location[77]. Due to the action of inhomogeneous and aniso- 
tropic turbulence in the downstream direction the existence of a pair 
of contra rotating vortices were also reported. It was noted that the 
streamwise vorticity had a small influence on the static pressure field, 
whereas its effect was more pronounced on the total pressure, mean 
velocity, and the Bernoulli's surfaces. 

More recently, Bhargava, Raj, and Boldman[78] reported detailed 
measurements of the mean surface pressure and the wall shear stress 
vector on both the surfaces in the blade end-wall corner region. The 
measurements were obtained on the same model used by Hazarika[47], 
In the comer region, the magnitude of the wall shear stress was found 
to increase significantly compared to the value far upstream of the 
blade leading edge. The surface static pressure was observed to be sig- 
nificantly influenced on the end-wall, whereas very small changes 
were observed on the blade surface. 


45 


OR’GINAL PAGE IS 

OF POOR QUALITY 

2.3 Wall Shear Stress 

An important parameter in the analysis of near wall mean and 
turbulent flow properties is the wall shear stress. Besides providing 
the information about the friction drag force, it is also a useful 
parameter commonly used in the near wall similarity models. In the 
past, several techniques have been developed for the measurement of 
wall shear stress. Brown and Joubert[79] described various wall shear 
stress measurement methods while studying the effect of secondary 
forces on the floating element balance. A detailed review on the exist- 
ing wall shear stress measurement techniques with their limitations 
and advantages among each other was presented by Winter[80]. More 
recently, Hanratty and Campbell[81), while reviewing various wall 
shear stress measurement techniques, have discussed in detail the 
mass transfer analogy method. In pursuit of measuring the wall shear 
stress with accuracy and without disturbing the flow a nonintrusive 
laser interferometry technique has also been reported in the litera- 
ture[82-84]. Based on the available literature and following the work of 
Winter[80] a classification of various existing wall shear stress mea- 
surement techniques is shown in Fig. 2.4. A brief description of the 
different wall shear stress measurement methods is presented here. 
However, a very elaborate discussion on the Preston tube method, 
which was used in the present study, is included in this section. 

2.3.1 Momentum Balance Technique 

The momentum balance technique can be used for both devel- 
oping and fully-developed flows. In a fully-developed flow through a 
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Fig. 2.4 Classification of Wall Shear Stress Measurement Techniques 


duct a momentum balance provides an expression for the calculation 
of wall shear stress; 

T W = kl (dp/dx)f (2.22) 

where x w is the average value of the wall shear stress around the 
perimeter of the duct, ki is the ratio of duct cross-sectional area to 
the wetted perimeter, and (dp/dx)f represents frictional pressure 
drop. Ferriss[85] found that the superposition of additional pressure 
gradient due to contraction or diffusion in a duct passage caused large 
errors in the calculated value of the wall shear stress (reported in ref. 
79). Secondly, whereas no circumferential variation in x w was observed 

for the fully-developed region, significant circumferential changes 
were noticed in the developing region [86]. It is evident from the avail- 
able literature[87-88] that this technique has been widely used for 
calibrating various wall shear stress measuring devices such as the 
Preston tube, heated element probe, etc. In a developing two-dimen- 
sional flow such as the boundary layer under variable external pres- 
sure, a momentum integral equation (reference 89) may be used for 
evaluating the value of wall shear stress. This equation can be expre- 
ssed as 


Uoo 2 de/dx + (29 + 8*) Uoo dUoo/dx = t w /p (2.23) 

which, in dimensionless form, can be written as 

d0/dx + (2 + H) 0/Uoo (dUoo/dx) = Cf/2 (2.24) 

The value of free stream velocity gradient in the streamwise direction 
is related to the streamwise external pressure gradient as 

- Uoo dUoo/dx = l/p(dp/dx) (2.25) 
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and can be accurately evaluated. It Is to be noted that the momentum 
integral Eqn. 2.23 or 2.24 is valid only for two-dimensional flows, and 
a suitable momentum integral equation is required in a three-dimen- 
sional flow. To evaluate the wall shear stress using Eqn. 2.24 it is nec- 
essary to calculate the derivative of momentum thickness. This can be 
done by measuring boundary layer mean velocity profile at a number of 
streamwise stations. However, a careful velocity profile measurement 
does not necessarily give rise to an accurate evaluation of the wall 
shear stress, mainly because the momentum thickness changes slowly 
with the streamwise distance x. 

2.3.2 Direct Force Measurement Technique 

In this technique a floating element experiences a force caused 
by the shear which causes the element to deflect. This deflection is 
directly measured and, through calibration, provides the magnitude of 
the wall shear stress. This method of measurement does not employ 
any assumptions regarding the similarity in the flow near the wall and 
is, therefore, suitable in any flow situation. Kempfl90] was probably the 
first to use direct force measurement technique for the study of a flat 
plate boundary layer in water. Smith and Walker[91] used this method 
for an incompressible two-dimensional boundary layer with zero pres- 
sure gradient in air. Most of the earlier studies were limited to the 
zero pressure gradient case where the secondary effects due to pres- 
sure gradient, misalignment of the floating element, etc., were ne- 
glected. Brown and Joubert(791 studied the effect of secondary forces 
caused by the pressure gradient. The maximum secondary force was 
found to be approximately 15% of the wall shear force. The direct 
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force measurement technique has been shown to be useful in three- 
dimensional flows where both the magnitude and direction of the wall 
shear stress are important. Tennant, Pierce, and McAllister[92] devel- 
oped an omnidirectional floating element balance to measure the wall 
shear stress in three-dimensional flows. The basic difference between 
this device and all of the other two-dimensional floating element de- 
vices was its ability to simultaneously measure the magnitude and di- 
rection of the wall shear stress. McAllister, Pierce, and Tennant[93] 
pointed out that, because indirect devices sire calibrated only in two- 
dimensional flows there is some doubt regarding their obtaining re- 
sults in three-dimensional flows. The measurements obtained by 
McAllister, Pierce, and Tennant also showed that their technique was 
suitable for flows with varying wall shear stress. In addition, the results 
indicated good agreement for the wall shear stress direction and the 
limiting streamlines direction obtained by an oil film technique [93 J. 
The direct force measurement technique, however, is deficient in the 
following respects: (i) the size of the sensing element must be large to 
experience sufficient force for deflection. For this reason, the mea- 
sured wall shear stress is the average value over the sensing element 
area and, therefore, the floating element is difficult to use in flow situ- 
ations where the wall shear stress varies rapidly from point to point in 
a small region (for example, corner region). The diameter of the 
sensing element in the work of McAllister, Pierce, and Tennant was 
2.86 cm, whereas Brown and Joubert[79] used a smaller size (1.9 cm) 
sensing element, (ii) Traversing of the device is difficult because of 
device's stability requirements. Therefore, the test model must move 
while the device stays stationary and this is very cumbersome when 
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the model is large, (iii) The gap around the sensing element and mis- 
alignment of the sensing element are another problems associated 
with such a device which introduce error in the wall shear stress 
measurement. Recently Acharya, Bomstein, Escudier, and Vokurka[94] 
studied the effect of varying gap size and sensing element thickness. 
From this investigation it was concluded that the main source of un- 
certainty in determining the wall shear stress was the buoyancy force 
resulting from the non-uniformity of pressure within the gap sur- 
rounding the sensing element. Tennant, Pierce, and McAllister[92] 
reported good agreement between the wall shear stress measure- 
ments obtained with the floating element balance and the Preston tube 
in two-dimensional turbulent flows. A detailed study on the compari- 
son of the wall shear stress measurements, in three-dimensional tur- 
bulent flows, with different techniques was reported by Pierce and 
Krommenhoek[951. Different techniques compared were the Preston 
tube, the heated element, and the floating element balance. The 
Preston tube results were found to be within 5% of the values obtained 
with the floating element balance. 

2.3.3 Wall Similarity Technique 

The techniques developed for measuring the wall shear stress 
using the velocity profile similarity near a wall are also termed as indi- 
rect methods in the literature. Various devices (e.g., the Preston tube, 
heated element, etc.,) for the measurement of wall shear stress have 
been developed in the past using the principle of near wall flow simi- 
larity. This principle is based on the observation that in a two-dimen- 
sional turbulent flow near a wall, the motion of the fluid is influenced 
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by the wall variables only. These wall variables are the wall shear stress 
(x w ). kinematic viscosity (v), density of the fluid (p), and the distance 
from the wall (z). From the dimensional analysis it can be shown that 
[79], 


U/u* = f (zuVvl (2.26) 

where U is the local mean velocity at a distance z from the wall and u* 
= (tw/p) 0 - 5 is called the frictional velocity. Prandtl termed this near 
wall behavior the "Law of the Wall". It is quite evident from Eqn. 2.26 
that if such a functional relationship is known, then local wall shear 
stress can be calculated from the measurement of the other variables 
(for example, mean streamwise velocity, U). Studies have been made 
to establish the range of validity for the law of the wall. These studies 
indicated that in a two-dimensional boundary layer with adverse pres- 
sure gradient the following equations may be applicable[79] 

Region (1), the viscous sublayer, U/u* = zu*/v (2.27) 

Region (2), the logarithmic region, U/u* = 1/K ln(zu*/v) + B (2.28) 

Region (3), Half-power region, U/u* = k2(Apz/u*2)0-5 

+fl(A p v/u*3] (2.29) 

where K, k2 and B are universal constants and Ap = l/p(dp/dx). 

The above division of a turbulent boundary layer is based on the crite- 
rion that the deviation in velocity profile from the law of the wall is 
caused due to the adverse pressure gradient. It is obvious that the law 
of the wall is valid up to the edge of Region (2), and in Region (3) non- 
wall effects simulated through variable Ap dominate. The range of Re- 
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gion (2) was described by 30v/u* < z < 1.41 u*2/Ap and Apv/u*3 < 
0.05(79]. The methods and devices based on the above discussed 
principle and used in measuring the wall shear stress are briefly dis- 
cussed below. 

2.3.3. 1 Velocity Profile Method 

Clauser[96] proposed a method in which the mean velocity pro- 
file measurement is used to evaluate the magnitude of the wall shear 
stress. This method uses the logarithmic law (Eqn. 2.28) and can be 
rewritten as, 

U/Uoo = (l/K)(Cf/2)0 5 ln(zlWv) + (Cf/2)0-5[(l/K)ln(Cf/2)0-5+B) (2.30) 

In this method, a chart (called Clauser's chart) is constructed by plot- 
ting U/Uoo against logio(zUoo/v) for different values of Cf using Eqn. 

2.30. Then, the measured mean velocity profile is plotted on this chart 
and compared with the constant Cf curves. The plotted points closest 
to a constant Cf line gives an estimate of Cf. A modification to Clauser’s 

method was proposed by Bradshaw[97]. This method avoids plotting 
the velocity profile on U/Uoo, logio( z Uoo/v) axes. Instead, in this sim- 
plified method of Bradshaw, one suitable reference point (the value of 
zu*/v) is chosen on the inner law curve and plotted on the axes U/Uoo, 
z (which are the usual coordinates on which mean velocity profile is 
normally plotted) for different values of Uoo/u*. The resulting curve on 
the U/Uoo, z axes has a well defined intersection with the measured 
velocity profile in the region of U/Uoo ~ 0.6 for the more usual value of 
Uoo/u*. The value of U/Uoo corresponding to the point of intersection 
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and U/u* selected at the reference point then gives the value of Cf 
using, 

Cf = 2(uVUoo)2 = 2[(U/Uoo)/(U/u*) refl 2 
It was suggested that with little practice only two or three reference 
points need to be plotted, in the manner described above, to find the 
intersection and Cf value. 

2.3.3.2 Heat Transfer Similarity Method 

In this method wall shear stress is calculated from the mea- 
surement of heat transfer from a heated element provided the thick- 
ness of the thermal boundary layer over the heated element remains 
within the range of the distance of validity of the law of the wall. In 
this case, the heat transfer correlates only to the wall variables. Based 
on dimensional analysis it can be shown that [79], 

Q/lkA0 = Nu = f [ pcp/k, lu*/v] (2.31) 

where Q is the rate of heat convected from the heated element, k, (i 
and Cp are fluid properties, 1 is the characteristic length of the heated 
element and A0 is the temperature difference between the heated el- 
ement and the free stream. The temperature difference (A0) is as- 
sumed to be small so that the natural convection is absent. For a given 
fluid, the Prandtl number Pr (= (iCp/k) is constant and, therefore, the 

measurement of Q and A0 is required to calculate the wall shear stress 
(t w ) provided the functional relation, Eqn. 2.31 is known. Fage and 
Falkner[98] appear to be the first to use the heat transfer technique in 
a laminar boundary layer flow. Ludwieg[99], Liepmann and 
Skinner[100] and Bellhouse and Schultz[101] have used heated ele- 
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ment technique in turbulent flows. It was shown in the previously 
mentioned literature that the heat loss from the sensing element is 
proportional to the cube root of the wall shear stress! 1 00]. The thin- 
film probe developed by Bellhouse and Schultz! 101] was also useful in 
calculating the wall shear stress fluctuations. The condition to calcu- 
late the fluctuating value of the wall shear stress is that the thermal 
boundary layer should not only lie within the law of the wall but should 
remain entirely within the viscous sublayer. A similar restriction ap- 
plies if the same calibration is to be used for both the laminar and the 
turbulent flows. The range of application of the heated element probe 
is expressed by the inequality[50] 

Pr/Cf > Nu » 1 

The above mentioned inequality, expressed in terms of Reynolds num- 
ber based on effective sensor length, can be rewritten as[50] 

6.6/(Pr)0 5 < u * i e / v < 64 Pr 

If the working medium is air, the above mentioned inequality reduces 
to 


7.8 < u* l e /v < 46 

where le is the effective sensor length. In the above relation, the lower 
limit is the result of the boundary layer approximation. The upper 
limit is due to the fact that the thermal boundary layer of the film lies 
within the viscous sublayer of the turbulent boundary layer so that the 
linear velocity profile could be accurately assumed. 

McCroskey and Durbinfl02] have used a dual sensor probe to 
calculate both the magnitude and direction of the wall shear stress. 
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This is essential in the measurement in a three-dimensional turbulent 
flow. The dual sensor flush mounted probe was used in a three-dimen- 
sional turbulent flow by Higuchi and Peake[103]. It was observed that 
due to the small size of the sensor there will be a minimum thermal 
penetration into the skewed near wall flow and thus such a device can 
accurately predict the wall shear stress direction. The main advantage 
of such a method is that the miniature probes can be built and flush 
mounted and the point measurement of the wall shear stress vector 
can be accurately made without disturbing the flow field. 

Gessner and Jones[104] reported the directional characteristics 
of the local wall shear stress in a fully developed turbulent flow thro- 
ugh non-circular ducts. The measurements were obtained in square 
and rectangular channels. For both the channels the skewness of the 
local wall shear stress vector was found to increase in the close vici- 
nity of the comer. The reason suggested was the presence of large 
secondary flows in the corner region. The results also showed that 
skewness of the local wall shear stress vector was greater along the 
longer wall than along the shorter wall of the rectangular channel at 
points equidistant from the comer. This result was in agreement with 
the measured secondary flow profiles. These secondary flow profiles 
suggested that most of the flow away from the corner was directed 
outward along the longer wall. 

2.3.3.3 Similar Flows About Obstacles 

The method of wall shear stress measurement in this case is 
based on the principle that the velocity field about any small obstacle 
immersed entirely in the law of the wall region will be influenced by 
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the wall variables. It is important to note that the obstruction must be 
small enough to have a negligible effect upon the growth of the bound- 
ary layer in which it is placed. Therefore, any pressure difference (for 
example, in the case of the Preston tube, this is the pressure differ- 
ence between the Preston tube pressure and the local wall pressure) 
AP depends on the wall variables which can be written in a functional 
form as, 

AP = f [ t w . P> v, 1 ] (2.32) 

where 1 is the characteristic length parameter of the obstacle (e.g., di- 
ameter of the Preston tube). Using dimensional analysis, Eqn. 2.32 
becomes! 1 05] 

APl2/ pv 2= fi[ Xw i 2 / pv 2] (2.33) 

The form of the functional relationship can be obtained by calibrating 
the device in a known flow condition such as fully developed flow in a 
pipe. Several devices have been developed based on the above dis- 
cussed principle. Commonly used devices are the boundary layer 
fence, the Stanton tube, the razor blade technique and the Preston 
tube. The Preston tube, which is widely used, was proposed by 
Preston[105J and is discussed in detail here. The Preston tube has the 
advantage that geometrically similar tubes can be produced and thus a 
single calibration can be used for different diameter probes. Patel[87] 
produced widely used calibrations of the Preston tube covering a wide 
range of flow conditions and tube sizes. Patel also noted that the origi- 
nal calibration suggested by Preston! 105] was in error. Many other in- 
vestigators have proposed calibration equations for the Preston tube 
which are tabulated in Table 2.2. The calibration equations in Table 2.2 
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TABLE 2.2 Preston Tube Calibration Equations 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Definition of Variables: 

x* = logio (APd 2 /4pv 2 ), d + = du*/v, AP = Pp - Pg. Ps ' Surface Static Pressure 

y* = J oglO (T w d 2 /4pv 2 ), u* = (T w /p)0-5. Pp . Preston Tube d - Preston Tube Dia. 

Pressure. 






were obtained under different flow conditions and different standards. 
These conditions are described in Table 2.3. The information in 
Tables 2.2 and 2.3 was also presented by McAllister! 106]. The calibra- 
tion equations (see table 2.2) provided by Patel were obtained for a 
two-dimensional flow with pressure gradient. Patel also suggested the 
limiting values of pressure gradient, both favorable and adverse, within 
which his calibration equations might be valid. His proposed limits 
were based on the value of non-dimensional pressure gradient para- 
meter defined as 

Ap + = Apv/u*3 

The limits on pressure gradients are: 

(a) for adverse pressure gradients 

with maximum error 3%: 0 < Ap + < 0.01, du*/v < 200 

with maximum error 6%; 0 < Ap + < 0.015, du*/v < 250 

(b) for favorable Pressure gradients 

with maximum error 3%; 0 > Ap + > - 0.005, du*/v < 200 

with maximum error 6%; 0 > Ap + > - 0.007, du*/v < 200 

Patel also pointed out that these limits are a rough guide only. 

Brown and Joubert[79] observed that the Preston tube gave sat- 
isfactory results for much higher pressure gradient values than those 
suggested by Patel. The pressure gradient limit given by Brown and 
Joubert is Ap + = 0.023, whereas Patel[87] suggested a limit of Ap + < 
0.015 with 6 percent accuracy. Brown and Joubert also noted that the 
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Preston tube can be used with higher values of Ap + and without loss of 
accuracy provided du*/v (Reynolds number based on the Preston tube 
outer diameter) is reduced (in their work du*/v ranged from 30 to 
50). The main drawback associated with the calibration equations of 
Patel is that they do not match at the end of each region. Also the 
third equation is of transcendental nature and requires an iterative 
method to solve it. To overcome this problem Head and Ram[109] 
presented tabulated results for the calibration equations of Patel in 
terms of two new variables which in functioned form are related as 

AP/x w = F [APd 2 /pv 2 ] 

This expression was also plotted on the Up/Uoo, logio(dUoo/v) axes 
using Cf as the parameter (very much similar to the Clauser’s chart) 
where U p /Uoo = (AP/0.5pUoo 2 ) 0 - 5 . 

BertelrudlllO] suggested a single calibration equation (see Table 
2.2) for the Preston tube in terms of two variables AP/t w and du*/v. It 
was noted that the variables used in Patel's equations are the loga- 
rithms of the non-dimensional pressure and skin friction and, there- 
fore, are less suitable for evaluating the measurement accuracy. It was 
shown that suggested new parameters gave better accuracy compared 
to the usual variables x* and y*. 

It is important to note that all the techniques developed and 
mentioned under the category of flow about the obstacles (see Fig. 2.4) 
require the knowledge of local surface static pressure to calculate the 
magnitude of the wall shear stress. Conventionally, the surface pres- 
sure is measured by providing a tap on the surface. At many locations 
of interest the provision of surface static tap may be impracticable, for 
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example, near the trailing edge of a thin airfoil model. Moreover, a 
large number of surface taps on a surface in a three-dimensional flow, 
where the wall shear stress changes from point to point, need to be 
provided (for example, blade end-wall corner region). As a result, new 
designs of the Preston tube have been proposed either to eliminate 
the need of the measurement of surface static pressure or measure the 
surface pressure with a modified Preston tube! 1 1 1,1 12). A modified 
Preston tube proposed by Gupta[lll] consisted of two different dia- 
meter tubes attached together which rest on the surface side-by-side. 
The tube with smaller diameter was chamfered side way at 45® to 
increase the measured pressure difference. In another design two 
chamfered tubes of smaller diameters were placed symmetrically 
about the main tube. This second design, it was suggested, could be 
useful in three-dimensional flows. The modified designs of the Preston 
tube proposed by Gupta are shown in Fig. 2.5. The design shown in Fig. 
2.5(c) was suggested for three-dimensional flows in which pressure 
difference between the outer tubes is first used to align the probe in 
the local flow direction and then the pressure difference between 
central and outer tubes gives the value of the wall shear stress. The 
pressure difference measured for devices (a) and (b) (see Figs. 2.5(a) 
and 2.5(b)) was found to be a constant fraction of the Preston tube 
pressure. These values are also given in Fig. 2.5. It is important to 
mention that this new device was calibrated against the Preston tube 
and, therefore, depends on the calibration of the Preston tube for cal- 
culating the value of wall shear stress. Moreover, no results were pre- 
sented for device (c). It is also felt that it would be interesting to cali- 
brate these devices against direct force measuring device. 
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AP m /APn * 0.714 


APm/ APn = 0.798 


4S° 



45° 



Fig. 2.5 Modified Designs of Preston Tube by Gupta[l 1 1) 
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Bertelrud[l 12] in an effort to avoid the use of wall taps to mea- 
sure the surface pressure suggested a new design for the Preston tube. 
This probe consisted of a 1 mm outer diameter tube with a closed 
front end which was fixed to the wall inside a 3 mm outer diameter 
tube. A 0.5 mm hole was drilled at a distance of 5d (where d is the di- 
ameter of the outer tube) from the tip of the probe. Since the static 
pressure measured with this probe was not at the probe tip, a correc- 
tion to the pressure difference obtained with this probe was sug- 
gested. It was also noted that the difference in measured static pres- 
sure and true pressure was a function of L/d where L is the distance 
from the probe tip to the location of static tap. 

The available literature also shows that the Preston tube can be 
used in three-dimensional flows provided it is aligned in the flow di- 
rection! 1 13,1 14). Mojola and Young! 1151 and Shabaka[36] have made 
use of the Preston tube in three-dimensional flows (comer flows). 
They have not mentioned about the alignment of the tube in the flow 
direction. Prahlad!114) observed that consistent values of the local 
wall shear stresses were obtained provided the Preston tube was 
aligned in the local flow direction. Moreover, it was noted that an ac- 
curate alignment of the tube to the streamline direction was very 
cumbersome and, therefore, the Preston tube should also be calibrated 
for its yaw characteristics. The yaw calibration of three types of 
Preston tubes was presented by Prahlad. These three designs were the 
normal Preston tube, a single tube chamfered at 45°, and two tubes 
chamfered at an equal angle of 45° (termed the Conrad probe). For 
each design different size tubes were investigated. The yaw character- 
istics of the Preston tube were found to be represented by an equation 
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(Pp( 0 t) - P s )* = C a (Pp{a 0 ) - P s )* n a 

where (P p (ao) - Ps)* = (Pp(ao) - p s) d 2 /4pv 2 , P p (ao) is the pressure 

obtained from the Preston tube when aligned in the flow direction, 
and P s is the surface static pressure. (P p (a) - P s )* is defined similarly, 
P p (a) being the pressure recorded by the Preston tube when oriented 
at an angle a to the local flow direction. C a and n a are constants and 

can be calculated for a given Preston tube by calibration. The main ad- 
vantage of the chamfered tube was reported to be its greater sensitiv- 
ity to yaw angle near a = 0 as compared to a normal Preston tube. It 
was also suggested that for accurate determination of the wall shear 
stress vector with the Preston tube both the pressure gradient com- 
ponents in the wall shear direction and normal to it should be small. 
The method suggested by Prahlad[l 14], although simple in principle, 
requires an iterative method for evaluating the wall shear stress mag- 
nitude and direction. 

Braggl 116] reported measurements of the local wall shear stress 
for a turbulent flow in the comer region formed by two flat plates. The 
wall shear stress had an inflexion point in the close vicinity of the cor- 
ner. At any streamwise location, the maximum wall shear stress value 
in the vicinity of the comer was comparatively smaller than the values 
measured far away from the comer. Similar trends in the variation of 
the wall shear stress were observed for zero as well as adverse pres- 
sure gradient cases. Bragg also derived a momentum integral equation 
for the comer flows and noted that some of the terms in the equation 
were difficult to measure within the acceptable accuracy. Shabaka[36] 
measured wall shear stress using the Preston tube in the corner re- 
gion of the simplified wing body junction. His measurements showed 
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that, in the comer region, the value of wall shear stress on the body 
surface (wind tunnel wall) near the wing leading edge was substantially 
higher than the measured wall shear stress far upstream of the wing 
leading edge. At any axial station in the comer region Cf increased 
from its asymptotic value (Cf value far away from the corner line at any 
axial station) to a maximum value as the comer was approached. This 
trend was different than the one observed in the work of Bragg. In the 
work of Bragg, as discussed earlier, the wall shear stress decreased 
from the asymptotic value as the comer was approached. The probable 
reason for this difference in Bragg’s and Shabaka's results is due to the 

fact that a horseshoe vortex existed in Shabaka's case. Shabaka also 
noted that the maximum value of Cf t at different axial stations, de- 
creased on both the surfaces of the test model in the downstream di- 
rection along the comer. The location of the maximum value of Cf on 

both the surfaces moved away from the comer line in the downstream 
direction. The direction characteristic of the local wall shear stress 
was not considered either in Bragg's or Shabaka’s work. These charac- 
teristics generally would be influenced by the secondary flows of the 
"first kind" and the "second kind". 

In the case of rotating three-dimensional turbulent boundary 
layer flow, the use of the Preston tube was reported by Anand[117J. 
Patel’s calibration equations were used to evaluate the magnitude of 
the wall shear stress. The Preston tube was aligned in the direction of 
the limiting streamline angles obtained from the flow visualization. 
The value of the wall shear stress was also calculated using Clauser's 
method and good comparison was reported. 
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2.4 Wall Pressure Fluctuations and Wall Pressure-Velocity 

Correlations 

An extensive review on the wall pressure fluctuation measure- 
ments was presented by Willmarth[52]. The available literature also in- 
dicates that only a few attempts[49, 53-55] have been made to mea- 
sure the correlation of wall pressure fluctuation with the velocity fluc- 
tuations in a turbulent boundary layer. Moreover, most of these earlier 
studies were limited to two-dimensional turbulent flows. However, a 
considerable amount of information was obtained on the structure of a 
turbulent boundary layer through these measurements. 

The first measurement of the wall pressure fluctuation under a 
turbulent boundary layer was reported by Willmarth[118]. The mea- 
sured ratio of the RMS value of wall pressure fluctuation (p' w ) to the 
free stream dynamic head (qoJ was 0.0035 and was constant over a 
wide range of Mach and Reynolds numbers. 

Later extensive studies on the wall pressure fluctuation mea- 
surements were reported by Willmarth and Wooldridge] 11 9] and Bull 
[120]. Willmarth and Wooldridge! 119] reported measurements of the 
wall pressure fluctuation on the surface under a thick two-dimensional 
turbulent boundary layer. The ratio of the root-mean-square wall pres- 
sure to the wall shear stress was measured to be 2.19. It was reported 
that the power spectra of the wall pressure scaled with the free 
stream velocity and the boundary layer displacement thickness. From 
the measurements of space-time correlations of the wall pressure 
fluctuation information was obtained about the convection and decay 
characteristics of the pressure producing eddies. The convection 
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speed varied from 0.56 to 0.83 times the free stream velocity. The 
lower convection speeds were obtained when pressure fluctuations at 
high frequencies were correlated. On the other hand, higher convec- 
tion speeds were measured when the low frequency pressure fluctua- 
tions were correlated. It was concluded from these observations that 
the large and small scale pressure producing eddies decay after trav- 
elling a distance proportional to their scale. 

Bull[ 120) investigated the wall pressure fluctuations in a two-di- 
mensional turbulent boundary layer with zero pressure gradient. Space 
time correlation measurements were made in both the broad and 
narrow frequency bands. The objective was to gain information about 
the small scale components of the wall pressure fluctuation, which 
could not be obtained earlier, because the size of the pressure trans- 
ducer in relation to boundary layer thickness in previous studies was 
large. In this study di/8* varied from 0.15 to 0.51 (where di is the 

diameter of the pressure transducer and 8* is the boundary layer dis- 
placement thickness). In the work of Willmarth and Wooldridge! 1 19] 
the transducers used had a value of di/8* = 0.33. Bull's! 120] measure- 
ments showed that the ratio, p'w/^w increased from 2.11 to 2.80 as 
the value of Reg varied from 6400 to 33,800. Bull also noted that the 
observed trend for p' w /x w in his study was in contrast to the results 
obtained in fully developed flow through pipe as reported by Corcos 

1121] and Bakewell et al.[122]. The Corcos!121] and Bakewell, et al. 

1122] studies, however, showed a decreasing trend for p' w /Xw with 

increasing Reynolds number. The wall pressure fluctuation field was 
divided into two families of wave-numbers! 120]. The high wave-num- 
ber family was associated with the turbulent motion in a constant 
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stress layer. The other family was a low wave-number family associated 
with the large-scale eddy motion in a turbulent boundary layer. The 
broad band convection velocity results indicated that for small separa- 
tion distances there was a strong contribution from the constant 
stress layer and also possibly from the transition region. As the sep- 
aration distance increased to a few displacement thicknesses the 
broad band convection velocity was found to increase to its asymptotic 
value of 0.825 U«> which is a characteristic of large scale structure. 
This observation suggested that the contribution to the wall pressure 
correlation from the fine scale motion in the near wall region de- 
creased because of the rapid dissipation and dispersion of pressure 
sources due to the presence of high shear. 

The spatial resolution of a pressure field associated with a local 
turbulent flow is limited by the finite size of the pressure transducer 
sensing element! 123]. Consequently, a lack of resolution in space leads 
to an apparent inability to resolve in time. Corcos[123] presented an 
analytical method using spectral theory to evaluate the response func- 
tion of the transducer. In this approach it was necessary to know in 
advance cross-spectral density function r(co,C), where 0 ) is the fre- 
quency and C is the position vector. Based on experimental evidence, 
he represented r as 

if C = ( §.t|) 

r[co,tO] = <|>[G)] A[(oyu c ] exp(- io)£/U c ) 

r[o),0,ri] = <>[co] B[GJT|/U c ] 
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Using the above expressions, a relation for the attenuation (ratio of 
measured to actual frequency-spectral density) at a given frequency 
caused by the finite size of the transducer was obtained. The theory 
developed for the attenuation of the mean square value of the wall 
pressure fluctuation was applied to the measurements of Willmarth 
and Wooldridge! 119], and it was shown that the transducer used in 
the study of Willmarth and Wooldridge could not resolve a large frac- 
tion of the total pressure signal. 

The correction to the power spectra measured by a finite size 
transducer was also suggested by Willmarth and Roos[124]. A similarity 
model similar to that suggested by Corcos[123] was used to evaluate 
the correction to the measured parameters. It was suggested that the 
similarity model for the cross-spectral density function was not valid 
for spatial separations less than 0.78* and dimensionless frequencies 
greater than 3.0. The RMS value of the wall pressure fluctuation was 
found to be considerably attenuated due to the finite size of the pres- 
sure transducer. For a vanishingly small size of the pressure trans- 
ducer, p' w / x w was calculated to be 2.66. For R/5* = 0.221 the RMS 

value of the wall pressure fluctuation measured was 2.20 (R is the ra- 
dius of the pressure transducer). 

The influence of transverse curvature under a turbulent boundary 
layer on the wall pressure fluctuation and the structure of turbulence 
was studied by Willmarth and Yang! 125]. Measurements of the longitu- 
dinal space- time correlation of the wall pressure fluctuation showed 
that the streamwise convection speed was similar to that in the plane 
boundary layer. The explanation for the unchanged convection velocity 
was that the pressure producing eddies in the boundary layer with 
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transverse curvature might be smaller. This implied that the eddies 
were near to the wall where the mean velocity was lower. Two reasons 
were suggested for the reduction in size of the turbulent eddies due to 
transverse curvature. The first reason was the increase in fullness of 
the velocity profile. The second reason was that the transverse scale 
decreased more for larger eddies than for smaller eddies because of 
transverse curvature. It was also concluded that the effect of trans- 
verse curvature was negligible on the RMS value of the wall pressure 
fluctuation. 

The effect of roughness on the wall pressure fluctuation was in- 
vestigated by Blake[126], The measurements were made on the walls 
with different roughnesses using a pinhole microphone having a better 
high frequency resolution. The space-time decay rate was noted to be 
higher. The lower convection velocities were measured over the rough 
walls. The reason for the lower convection velocities was due to the 
retardation of the mean boundary layer motion by the surface rough- 
ness. It was also reported that, at moderately high frequencies where 
the pressure sources were nearer (z < 0.18) and the convection veloc- 
ity did not reach its asymptotic value, the spectra for all the walls be- 
haved as to -1 . This also required that the pressure spectrum level and 

wave-numbers would be determined by the wall shear stress and the 
source position, respectively. The value of p' w / T w was noted to be in- 
fluenced by the type of roughness. The RMS value of the wall pressure 
fluctuation for densely packed small scale and densely packed large 
scale roughness elements was smaller than for sparsely packed small 
scale roughness elements or smooth walls. On the other hand, for 
sparsely packed small scale roughness elements the RMS value of the 
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wall pressure fluctuation was larger than for the smooth walls. It was 
also noted that the roughness separation influenced the large scale 
turbulence structure, whereas the roughness height influenced the 
medium and very small scale turbulence structure. 

Schewe(127] studied the structure of turbulence and showed 
that the value of diu*/v = 19 (dimensionless diameter of the trans- 
ducer) was sufficient to resolve the essential structure of the turbulent 
wall pressure fluctuations. The power spectra measured with the sma- 
llest transducer showed the power law as - ci)”7/3. This law is similar to 
the one predicted for locally isotropic turbulence. The pressure struc- 
tures with high amplitude and the shape of short pulses were recog- 
nized using signal averaging in the time domain. These structures 
were found to have a characteristic frequency (o+) and longitudinal 
wave length (X+) of 0.52 and 145, respectively. The convection velo- 
city associated with such structures had a mean value of 0.53 Uoo. The 
RMS value of the wall pressure fluctuation (non-dimensionalized by the 

free stream dynamic head) increased considerably from 0.5% to 1% as 
the di + reduced to the smallest transducer used. Also the difference 

in the RMS value measured with the smallest transducer and the value 
corresponding to di + = 0 was only 4%. Based on this result it was 

concluded that the smallest transducer used in his study could resolve 
the pressure structure essential to turbulence. 

The wall pressure fluctuation measurements downstream of two 
rigid ribbons mounted in tandem and immersed in a turbulent bound- 
ary layer were reported by Beeler{128]. These ribbons are known as 
large eddy breakup devices (LEBUs). The intensity of wall pressure 
measured at a location (where the Cf value decreased by approximately 
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25%) downstream of the LEBUs showed a decrease of about 12.5% 
compared to the case when LEBUs were absent. This decrease in the 
wall pressure intensity also implies a reduction in noise generated due 
to the wall pressure fluctuations associated with a turbulent boundary 
layer. 

The wall pressure field for flow separation and re-attachment 
processes in a turbulent boundary layer was experimentally investi- 
gated by Farabee and Casarella[129]. Two types of flow were investi- 
gated: flow over a forward-facing step and flow over a backward-facing 
step. The experimental results showed that the separation-re-attach- 
ment process produced large amplitude, low frequency wall pressure 
fluctuations. The wall pressure fluctuations were found to be largest at 
re-attachment where the highly turbulent free-shear flow impinges on 
the surface. In particular, RMS values of the wall pressure fluctuation 
at re-attachment for the forward-facing step were 10 times larger than 
the pressures measured for an equilibrium flow. However, for the 
backward-facing step, with the same step height as in case of forward- 
facing step, RMS values of the wall pressure fluctuation were only 5 
times larger than the pressure measured for an equilibrium flow. This 
large increase in RMS values of the wall pressure fluctuation was at- 
tributed to the increase in the low frequency content of the pressure 
spectrum. 

Recently, measurements of RMS values of the wall pressure fluc- 
tuation and its frequency spectra in a wing-body junction were re- 
ported by Hasan, Casarella and Rood(130]. The wing-body junction was 
formed by mounting a wing, of semi-elliptic nose with constant thick- 
ness mid-body and a NACA 0020 tail profile, normal to the wind-tun- 
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nel side wall. In upstream region, between the saddle point and the 
wing leading edge and on the stagnation streamline, a significant in- 
crease in the wall pressure intensity (five times its far upstream value) 
was obtained. This observed increase was partially attributed to the 
impingement of the re-circulating flow within the comer region. The 
demarcation line of the horseshoe vortex around the wing was ob- 
served to coincide with the peak values of the intensity of wall pres- 
sure fluctuation at different streamwise locations in the comer region. 
Also, at any streamwise location in the comer region on the body sur- 
face, the frequency spectra of the wall pressure fluctuation showed an 
increase in spectral content towards the wing surface in the low fre- 
quency range (80-200 Hz). 

The effect of pressure gradient on the wall pressure fluctuation 
was also studied[49,131). Schloemer[131] noted that, in comparison 
to zero pressure gradient case, the ratio of RMS value of the wall pres- 
sure fluctuation to the free stream dynamic head was more for adverse 
pressure gradient and less for favorable pressure gradient. On the 
other hand, the normalized RMS value of the wall pressure fluctuation 
(non-dimensionalized by the wall shear stress) was also larger for the 
adverse pressure gradient compared to the zero pressure gradient 
case. The increased value of p' w /x w was higher compared to the in- 
crease observed for the ratio p’ w /c[°o. The spectral density was found 
to be influenced by the pressure gradient especially at low frequencies. 
At low frequencies the spectral density increased for both adverse and 
favorable pressure gradients compared to the zero pressure gradient 
case. At high frequencies the spectral density was approximately the 
same for the adverse and zero pressure gradient cases. 
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Bradshaw[49] noted that the low wave number pressure fluctua- 
tions are generated outside the inner layer. He suggested that this 
pressure fluctuation field consisted of all the spectral density at low 
wave number and contributed to the inactive motion in the inner re- 
gion. Bradshaw also suggested that the active inner layer motion pro- 
duced a high frequency contribution to the spectra which varied as 
x w 2 /k up to wave numbers of the order of the inverse of the sublayer 

thickness. The influence of the different flow conditions and the effect 
of the size of the pressure transducer on the wall pressure fluctuation 
is presented in Table 2.4. 

The correlation of fluctuating wall pressure and the streamwise 
component of velocity fluctuation was first reported by Seraflni[53J. 
The advantage of such measurement is to evaluate the contribution to 
the wall pressure fluctuation from different regions of the turbulent 
flow field. Two regions were identified where the streamwise compo- 
nent of the velocity fluctuation contributed to the wall pressure fluctu- 
ation. These two regions were the inner region near the wall and the 
outer region. It was reported that the contribution to the pressure 
fluctuation at the wall from the outer region was less than one-half of 
that from the inner region. 

The measurements of the wall pressure-velocity correlations 
were also reported by other investigators[49, 54.55J. The measure- 
ments of correlations of the wall pressure fluctuation with all the three 
components of velocity fluctuation were obtained and a qualitative 
model for the structure of turbulence near the wall was proposed 
[54,55]. It was proposed that the stretching of the vorticity produced 
by the viscous stresses in the sublayer leads to the fluctuation of wall 
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TABLE 2.4 Influence of Different Flow Conditions and Pressure Transducer Size 
on the Root-Mean-Square Value of Wall Pressure Fluctuation 
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pressure and velocity. It was also observed that the maximum magni- 
tude of the correlation coefficient (Rpu) was 0.3. 

Very recently, Kobashi and Ichijo[132] investigated the structure 
of coherent motions in a two-dimensional turbulent boundary layer. 
The measurements of the wall pressure fluctuation field and its corre- 
lation with streamwise and normal velocity components were ob- 
tained. From these measurements two types of coherent motions were 
identified. The first type included the large scale motions which pre- 
vailed through and outside the boundary layer. These large scale mo- 
tions were observed to be initiated by the instability of the mean flow 
and rotate in the direction of the mean flow shear. The second type of 
coherent motion was the small scale motion (also known as burst) 
which was found to be limited to a narrow region near the wall. The 
small scale motions were like vortices, which rotated against the mean 
flow shear. A vortex ring model of coherent structure was proposed 
and it was concluded that the burst was bounded by two pairs of 
counter rotating vortices which were inclined forward. The ejection 
and sweep which characterized the burst phenomena are the flows 
which are induced by the paired vortices in front and back of the 
burst. 

The influence of the pressure gradient on the wall pressure-ve- 
locity correlation was reported by Bradshaw[49]. The correlation coef- 
ficient Rpu was higher very close to the wall and lower away from the 
wall compared to the zero pressure gradient case. For a small distance 
z, correlation data was also presented in a narrow frequency band. At 
the lowest frequencies, the value of the correlation coefficient was 0.6. 
In the medium frequency range the value of the correlation coefficient 


77 


(Rpu) varied very little and was ~ 0.2 to 0.25. It was suggested that the 
variation from 0.2 to 0.25 corresponded to the wall pressure fluctua- 
tion generated in the inner layer. It was also noted that the v-compo- 
nent fluctuation in the inner region was less affected by the wall pres- 
sure fluctuations than the u-component. 

The discussion presented in sections 2. 1-2.4 suggests that a de- 
tailed experimental investigation of the influence of free stream tur- 
bulence on the development of a turbulent boundary layer downstream 
of distributed roughness, which is placed in the close vicinity of the 
leading edge of the flat plate, is needed. There is also a need to obtain 
a detailed measurement of the wall shear stress vector field in the 
blade end- wall corner region. The detailed knowledge of the wall 
pressure field and its correlation with the fluctuating velocity field is 
also desired to understand the near wall turbulence structure in three- 
dimensional flows. 


78 


CHAPTER m 


EQUIPMENT, INSTRUMENTATION, EXPERIMENTAL METHOD 

AND DATA PROCESSING 

The details of the equipment, instrumentation, methods of mea- 
surement, and the data processing techniques are described in this 
chapter. The present experimental study includes measurements of 
the mean velocity, wall shear stress, surface pressure, wall pressure 
fluctuation, and its correlation with the fluctuating velocity upstream 
and in the blade end-wall corner region. 

3.1 Equipment and Instrumentation 

3.1.1 Wind Tunnel 

A low subsonic open circuit wind tunnel, designed and cons- 
tructed at the Turbomachinery Laboratory of the City College of the 
City University of New York, was used for the experiments. A centri- 
fugal blower fitted with variable inlet guide vanes was driven by a 25 
H.P., 1770 rpm A.C. motor to generate the flow in the wind tunnel. Air 
velocity in the 46 cm x 46 cm test-section could be varied from about 
5 m/s to 35 m/s by adjusting the inlet guide vane setting. The non- 
uniformity in the free stream mean velocity in the test- section was 
within ± 0.5% at the maximum flow rate. The free stream turbulence 
intensity in the test-section was found to vary from 0.6% at a velocity 
of 6 m/s to 0.1% at a velocity of 20 m/s and higher. An inlet filter was 
used to separate dust from air passing through the wind tunnel. 


79 


9.1.2 Test-Model 

The blade end-wall comer test-model was constructed by moun- 
ting a blade with a symmetrical airfoil section on a flat plate. Two flat 
plates were designed to obtain the required measurements. The first 
flat plate (hereafter, will be referred to plate 'A') was constructed from 
two, 91 cm x 46 cm x 0.64 cm, 2024 T4 bare aluminum plates as 
shown in Fig. 3.1. Four 82.6 cm x 1.3 cm x 0.64 cm aluminum ribs 
were placed lengthwise between the plates and screwed together to 
produce flat surfaces. The holes were filled with epoxy and sanded to 
form smooth surfaces. The leading edge of the flat plate was made 
from a 1.9 cm diameter aluminum bar, whereas the trailing edge was 
constructed from a 2.5 cm x 1.9 cm bar. Six static pressure holes were 
provided on each side of the assembly with three holes in each row 
length-wise and two rows on each plate. Twelve polyethylene tubes 
connected to the wall static holes passed through the space between 
the plates and emerged from the side near the trailing edge of the flat 
plate assembly. These surface static taps were mainly provided so that 
the test-model could be aligned perfectly at zero incidence to the 
main flow. Flat plate A' was employed for the measurements of the 
mean velocity profiles, the wall shear stress, and the surface static 
pressure with the surface static tube. 

The second flat plate (hereafter, called plate 'B') was constructed 
from a single 2024 T4 bare aluminum plate as shown in Fig. 3.2(a). 
The leading edge of the flat plate was made circular in shape by using 
a standard 3/4 in. cutter. At four axial locations provisions were made 
to fix interchangeable circular disks shown in Figs. 3.2(b) and 3.2(c)). 
Two instrumentation and four dummy disks were constructed. In one 
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Fig. 3.2 (Contd.) Flat Plate Assembly- Plate B' 
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12.0 




of the instrumentation disks (hereafter, called Disk No. 1) five Kulite 
pressure transducers were installed (Fig. 3.3). The transducers were 
soft mounted using silicon rubber (RTV Dow Corning 3145). In the 
second instrumentation disk (hereafter, called Disk No. 2) nine sur- 
face pressure taps of 0.8 mm (1/32 in.) diameter were provided along 
its diameter (Fig. 3.4). The main dimensions of the flat plate 'B' were 
the same as those of the flat plate 'A' described earlier. A large number 
of surface static taps were also provided on the flat plate 'B' to obtain 
surface pressure distributions (Fig. 3.2(a)). 

Based on the compromise of obtaining a thick boundary layer on 
the blade and its stall characteristics, a NACA 65-015 base profile was 
chosen for the blade. The chord length of the blade was 25.4 cm. The 
leading and the trailing edges of the blade had diameters of 1.27 and 
0.25 cm, respectively (Fig. 3.5). The blade was made in two sections 
each with 21.9 cm span and a 0.95 cm hole in the spanwise direction 
at the position of maximum thickness on the chord line. A threaded 
rod passing through the hole in the flat plate and the blades held the 
entire assembly together. The distance between the leading edge of 
the flat plate and the blade was 62 cm. The initial seven cm from the 
leading edge of the flat plate was made rough by gluing a No. 4 sand 
paper to trigger the transition and artificially increase the boundary la- 
yer thickness on the flat plate (Fig. 3.6). The test-model assembly was 
bolted to the test section with the help of two threaded rods. The top 
and bottom test-section walls parallel to the flat plate were removed 
10 cm ahead of the blade leading edge to facilitate probe traverses. 
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Fig. 3.5 Airfoil for Test-Model 
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3.1.3 Traverse Mechanism and Rotary Device 

3. 1.3.1 Traverse Mechanism: 

The traverse mechanism was constructed with two compound 
slides held perpendicular to each other with an angle plate (Fig. 3.7). 
The assembly was bolted to a table anchored to the floor. A device to 
rotate the probe about its own axis was bolted to the connecting bar. 
The connecting bar was clamped to the compound slide with its axis 
perpendicular to the axis of the test-section. 

The least count of the compound slides were 0.0254 mm. The 
total traverse distance parallel to the axis of the test-section was 33 
cm, perpendicular to the axis. In the horizontal and vertical direc- 
tions, the traverse distance was 20 cm. The traverse distance of the 
probe in all the three directions could be extended by 5 cm by chang- 
ing the clamping position of the connecting bar on the vertical surface 
of the compound slide. Subsequent extension in traverse distance was 
possible by changing the position of the angle plate on the horizontal 
compound slide. 

3. 1.3.2 Rotaxy Device 

The rotary device was machined from brass stock in two prin- 
cipal components. The base, which was bolted to the supporting rods, 
carried a rotating part in the shape of a frustum of a cone with a hole 
through its axis and a chuck to hold the guide tube. Three adjustable 
spring loaded steel balls provided sufficient pressure between the base 
and the rotating part to hold them together snugly without impeding 
rotary motion. Two supporting rods bolted to the base held a guide 
block which had a hole aligned with the axis of the guide tube holder 
chuck. The protractor mounted on the rotating part and the vernier 
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Fig. 3.7 Photograph of Traverse Mechanism and Rotary Device for the 
Wall Shear Stress Measurements With the Preston Tube 



mounted on the base allowed the change in angle to be measured to 
the least count of 0.2 degrees. The Preston tube and guide tube as- 
sembly were held by the guide tube holder chuck. The rotary device 
was mounted on the traverse mechanism such that the axis of rotation 
made an angle of 90® with the axis of the test-section (Fig. 3.7). 

3.1.4 Grids 

A combination of two grids was used to create free stream tur- 
bulence level of 1.5% at 107 cm downstream of the grids (location of 
the blade leading edge) to avoid a separation bubble and increase the 
boundary layer thickness over the surfaces of the blade. The grids se- 
lected for this purpose were 1 mesh 2 mm diameter galvanized iron 
wire mesh followed by a 4 mesh 0.6 mm diameter galvanized iron wire 
mesh (Fig. 3.8). The wires of the grids were inclined 45® to each 
other. The grids were placed 45 cm ahead of the leading edge of the 
flat plate. The non-uniformity in the free stream mean velocity in the 
test section was less than 0.5% at a velocity of 27 m/s. 

3.1.5 Probes and Instrumentation 

The mean wall shear stress measurements were obtained with a 
specially constructed Preston tube. A Preston tube of 0.81 mm (0.032 
in.) outer diameter with an inner-to-outer-diameter ratio of 0.61 was 
used to measure the wall shear stress. The Preston tube was bent as 
shown in Fig. 3.9 so that the tip of the probe could be concentrically 
rotated about its axis of rotation. This was essential to align the probe 
tip in the flow direction at the measurement location. The stem of the 
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Fig. 3.9 Sketch of the Preston Tube Used in the Measurements 
of Wall Shear Stress 
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probe was soldered to increasing diameter tubes such that the probe 
could be mounted in the available fixtures. 

The mean surface static pressure was measured by two methods. 
In the first method, the wall static pressure was measured using the 
United Sensor Corporation Pitot-static probe of 1.59 mm (1/16 in.) 
diameter, resting on the surface. In the second method, a large num- 
ber of wall taps (also called a wall piezometer) provided on the surface 
of flat plate 'B' of the test-model and Disk No. 2 were used to measure 
the surface static pressure. 

A Kiel probe also made by the United Sensor Corporation was 
used to measure the reference total pressure which was maintained 
constant during the experiment. This was essential to maintain the 
free stream velocity almost the same at different times of data collec- 
tion. 

A total head probe of 0.81 mm (0.032 in.) outer diameter along 
with the pitot-static probe, described earlier, were used to obtain the 
mean velocity profile data. 

The different pressure measuring probes, discussed above, were 
connected to pre-calibrated variable reluctance type differential pres- 
sure transducers (Model DP 15) with a range of ± 0.5 PSID and CD 15 
carrier demodulator. The output was read on a digital voltmeter. 

The wall pressure fluctuation measurements were obtained with 
an integrated sensor pressure transducer (Model XCS-062-5G) manu- 
factured by Kulite Semiconductor. This transducer is a piezoresistive 
type consisting of an integrated silicon sensor (also called a diaph- 
ragm). A four arm wheatstone bridge is diffused in the surface of the 
silicon sensor using a photolithographic process. The transducer con- 
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sisted of a reference tube for gage pressure measurement. The pres- 
sure transducer was also provided with the temperature compensation 
module. The pressure transducers were the smallest available in size 
and had a diameter of the pressure sensitive area of 0.71 mm (0.028 
in.). Some of the important specifications of a typical piezoresistive 
type pressure transducer supplied by the manufacturer are given in 
Table 3.1. In Table 3.2 the characteristics of the five transducers used 
in the investigation are provided. 

A single normal sensor hot wire probe was used to obtain the 
streamwise fluctuating velocity component. The sensor was made from 
platinum coated tungsten wire with diameter and length of 4pm and 1 
mm, respectively. 

3.1.6 Calibration Equipment 

A Thermo System Inc. calibrator model 1125 was used to cali- 
brate the variable reluctance type pressure transducers against a mi- 
cromanometer with 0.0254 mm of water pressure resolution. A repre- 
sentative calibration curve for two differential pressure transducers is 
shown in Fig. 3.10. 

3.2 Methods of Measurement 

During the course of the flow visualization study which was con- 
ducted by Hazarika(47] on the same test-model, a separation bubble 
was observed on the blade surface in the blade-end-wall corner region. 

As a result, the entire experimental investigation was performed 
in two parts. In the first part, termed Case A, different flow variables of 
interest were measured without the presence of separation bubble. 
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TABLE 3. 1 Typical Specifications of a Piezoresistive 

Transducer (Model No. XCS- 062- 5G) as Supplied 
by Kulite Semiconductor Products Ltd. 


Parameter 
Bridge Type 

Rated Pressure 
Full Scale Output 
(nominal) 

Excitation voltage 
Input Impedance (min.) 
Output Impedance (nom) 
Combined Nonlinearity 
and Hysteresis 
Repeatability 
Natural Frequency 
Zero Balance 

Compensated Temperature 
Range 

Maximum Change of 

Sensitivity with Temperature 

Acceleration Sensitivity 

% FS/g 

Perpendicular 

Transverse 

Resolution 


Specification 

Fully active four arm Wheatstone 
bridge diffused into silicon 
diaphragm 
5 psig 
225 mv 

15 VDC 
1200 O 
2500 Cl 
0.5% BFSL 

0. 1% FS 
150 KHz 
3% FS 
80 F - 180 F 

+ 2% / 100 F 


0.005 

0.0005 

Infinite 
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TABLE 3.2 Specifications of the Pressure Transducers (Model XCS- 062- 5G) Used in 
the Investigation 
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Fig. 3.10 A Typical Calibration of Variable Reluctance Type Differential 
Pressure Transducers 


The flow variables measured included mean value of the surface static 
pressure and wall shear stress, mean velocity profiles, the fluctuating 
wall static pressure, and the streamwise fluctuating velocity compo- 
nent. In the second part, termed Case B, measurements were obtained 
in the presence of separation bubble. The flow parameters of interest 
in case B were mainly the surface static pressure and the wall shear 
stress. In this section, the methods of measuring the flow parameters 
discussed above along with their data processing procedures are de- 
scribed. 

3.2.1 Reference Flow Condition 

The reference total pressure was measured at the streamwise 
station which was located at 50.8 mm upstream of the blade leading 
edge and at 76 mm above the end-wall (flat plate) to determine the 
reference free stream velocity, Ue* (Fig. 3.11). At this location (i.e., x = 

- 50.8 mm, y = 0 mm, and z = 76 mm) the static pressure was found 
to be equal to the atmospheric pressure. The reference total pressure 
was maintained nearly constant for the entire experiment. The free 
stream velocity obtained from the measured reference total pressure 
was 27.3 ± 0.3 m/s during the ehtire experiment. The distances of dif- 
ferent streamwise stations with reference to the airfoil leading edge 
and the plate leading edge are shown in Table 3.3. Also included in 
Table 3.3 are corresponding values of x/c (c is the blade chord length). 
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TABLE 3.3 Location of Streamwise Stations 


Station 

No. 

Streamwise 
Distance xl 
(mm) 

Streamwise 
Distance x 
(mm) 

x/c 

(c = 254 mm) 

0 

366.0 

- 254.0 

- 1.0 

1 

518.4 

- 101.6 

- 0.4 

2 

543.8 

- 76.2 

- 0.3 

3 

569.2 

- 50.8 

- 0.2 

4 

594.6 

- 25.4 

- 0.1 

5 

645.4 

25.4 

0.1 

6 

696.2 

76.2 

0.3 

7 

747.0 

127.0 

0.5 

8 

772.4 

152.4 

0.6 

9 

823.2 

203.2 

0.8 

10 

871.0 

251.0 

0.99 


NOTE : 

(i) xl is the streamwise distance with reference to the plate leading 
edge. 

(ii) x is the streamwise distance with reference to the blade leading 
edge. 
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3.2.2 Mean Velocity Profile Measurement 
3.2.2. 1 Two-Dimensional Flow 

A detailed experimental investigation was undertaken to study 
the influence of free stream turbulence on the development of a tur- 
bulent boundary layer developing downstream of the distributed roug- 
hness strip. In order to investigate the above mentioned effect, the 
following four cases were studied: (i) the development of a turbulent 
boundary layer with natural transition and no free stream turbulence 
(Case 1), (ii) the influence of free stream turbulence on the develo- 
pment of a turbulent boundary layer with natural transition (Case 2), 
(iii) the development of a turbulent boundary layer downstream of an 
initially distributed roughness (Case 3), and (iv) the combined effect of 
initial roughness and the increased level of free stream turbulence 
(Case 4). 

For each case mean velocity profiles were measured at seven 
streamwise locations on the flat plate centerline. The two dimension- 
ality of the flow was checked by measuring the variation of the surface 
static pressure and the wall shear stress, at the last measurement lo- 
cation, over a distance of approximately 75 mm from the plate center- 
line in the spanwise direction. From the measurements of the above 
mentioned two parameters, at the last measurement location, it was 
found that the variation in surface static pressure was very negligible 
(about 0.12 mm of water column). However, the value of wall shear 
stress changed within ± 2%. 

The boundary layer thickness, 8, was obtained, from the mean 
velocity profile data, corresponding to U/Ue - 0.995, where U is the 
streamwise mean velocity within the boundary layer and Ue is the local 
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free stream velocity. The value of 5 was also calculated following Coles's 
work[ 133]. Using the combined law of the wall and wake to represent 
velocity profile completely, as proposed by Coles, a relation of 5 in 
terms of displacement thickness (8*) and the wake strength parame- 
ter (n) can be written as, 

8/8* = K (U e /u*)/(l + n) (3.1) 

where K is the Von-Karman constant and u* is the friction velocity. 
Once the value of 8*. u* and FI are known, the value of 8 can be calcu- 
lated from Eqn. 3.1. It is noted, that in Eqn. 3.1 the value of wake 
strength parameter, n, and u* are needed. Following Coles's[133] work 
the value of n was calculated using the relation, 

2n - ln(l + n) = K(U e /u*) - ln(8*U e /v) - KB - ln(K) (3.2) 

whereas, the value of u* was calculated from the measured wall shear 
stress. The values of the constants K and B used in Eqn. 3.2 were 0.41 
and 5.0, respectively. These values were suggested by Coles[134] in 
1968 Stanford conference and are generally considered standard in 
turbulent boundary layer work. 

The mean velocity profile data within the boundary layer, ob- 
tained using a pitot probe, was corrected for the pitot displacement 
effect as suggested by MacMillan[135]. 

The boundary layer integral parameters, viz., displacement 
thickness (8*), momentum thickness (0), Clauser's shape parameter 
(G), and Clauser's universal integral thickness (A) were calculated us- 
ing the standard definitions. The contribution to integrals for calcu- 
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lating various boundary layer integral parameters in the range 0 < 
zu*/v < 50 was evaluated following the work of Coles[134J. 

The various boundary layer parameters calculated were com- 
pared with the available empirical relations and also with the results of 
other investigators. 

3.2.2.2 Three-Dimensional Flow 

The measurements of mean velocity profile, using the pitot tube 
in combination with the pitot- static tube, were obtained in the flow 
field upstream of the blade leading edge. Five streamlines were se- 
lected upstream of the blade leading edge covering a distance of 0.4 
and 0.3 of blade chord length in the streamwise and lateral directions, 
respectively. The measurement locations for the mean velocity profile 
data on five streamlines are shown in Fig. 3.12. In Fig. 3.12 the symbol 
'S' corresponds to the streamline, whereas numeric subscripts corre- 
spond to the measurement location. The first numeric subscript cor- 
responds to the streamline number, whereas the second subscript 
corresponds to the number of the streamwise station. Traces of the 
above mentioned streamlines were obtained from the flow visualization 
studies conducted by Hazarika[47], 

The coordinates of the measurement locations on five stream- 
lines are tabulated in Table 3.4. The data of the limiting streamline di- 
rections, obtained from the flow visualization[47], at each location of 
measurement is also included in the Table 3.4. The static and total 
pressure probes were aligned with the flow direction at an angle at 
each measurement location as given in Table 3.4. For the flow angles < 
4^> the probes were aligned parallel to the test model (flat plate) cen- 
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Test Model Centerline 
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Fig. 3.12 Locations of Mean Velocity Profile Measurement Upstream of 
the Blade Leading Edge (For Coordinates of Measurement 
Locations, See Table 3.4) 


TABLE 3.4 Coordinates of Measurement Locations for Mean Velocity 



II 

CQ 

— 

o 
1 n 

CS 

« 

o o 
o o 

oo 
o m 

cs d 
*-< cs 



CS 

o q 

oo in 
co ^ 

— 

o o 

f— ( i— < 

CO 

in 
o n 

d oo 

CD 

co 

ti 

CQ 

O 

00 

6 

in 

o o 
o o 

m 

1 

00 
o in 

N 00 
CS 

o 
o in 

CD oo 
in 

oo 
o co 

in cb 

CD 

cs 

ll 

CQ 

0 
cs 

cb 

1 

o o 
o o 

o o 

<— i 

in 
O N 

in ^ 

CN 

<N 

o 

C 0 ID 
ID 

<o 

© ~ 

CO co 
CD 

r— < 

it 

CQ 

o 

CD 

o 

o o 
o o 

CO 

o oo 
n o 

m 
o cs 

cs 

CS 
o Tt 

p-H LO 

m 

(N 

o q 
o in 

CD 

O 

II 

CQ 

0 
in 

CS 

1 

o o 
o o 

1 1 

i i 

i i 

i t 


a 

a 

X 

<D 

a § 
fis 

a 

T 

_ u 

Is ^ 
1 3 
>> a 

T? 

OJ 

6 § 
is 

>. a 

<u 

¥ ™ 
C «■> 

3 3 
>> 8 

'z: *-• 

£ ttf) 

3 3 
>> 8 

9 

c n 


rH 

n 

< 

CS 

II 

< 

co 

it 

< 

ll 

< 

m 

ii 

< 


106 


Note: Sab represents the streamline in which letter "A" corresponds to the 

streamline number (A = 1 , 2 5) and "B" corresponds to streamwise 

station number (B = 0, 1 4) 

























terline because the pitot-static and the total head probes are known to 
be insensitive to flow angles less than 5® (based on calibration curves 
supplied by United Sensor Corporation). 

The velocity can be calculated from the total and static pressure 
measurements and is given by the relation 

U = {2Ci(V t - V s )/p}0.5 

where Vt and Vs are transducer outputs in volts of total and static 
pressures, respectively, and Ci is the calibration constant of the pres- 
sure transducer. 

3.2.3 Wall Shear Stress and Wall Static Pressure Measurements 
3.2.3. 1 Two-Dimensional Plow 

Two-dimensional wall shear stress data were collected on the 
flat plate 'A' with the blade removed. The traverses of the Preston tube 
and the static pressure probe were made at seven streamwise stations, 
station 1, and stations 5 to 10. The traverses in the Y-direction cov- 
ered a distance of nearly 76 mm from the test-model centerline. The 
grid showing the locations of the two-dimensional wall shear stress 
measurements is shown in Fig. 3.13. The Preston tube was aligned 
parallel to the test-model centerline because in this case flow stream- 
lines were parallel to the test-model centerline. The value of wall 
shear stress was also calculated using the mean velocity profile data 
with the help of Clauser's chart (see discussion on the method in sec- 
tion 2.3.3. 1). The measured values of the wall shear stress were com- 
pared with the well established empirical relations available in the lit- 
erature. 
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The yaw sensitivity of the Preston tube was examined in two- 
and three-dimensional flows. The yaw angle range covered was ± 30 
degrees. A detailed analysis was done to evaluate the error due to mis- 
alignment of the Preston tube with the local flow direction. 

The surface static pressure was measured with the static tube 
resting on the surface. The use of this technique (called surface tube 
method), to measure surface static pressure, was demonstrated by 
Sproston and Goksel[136) in a two-dimensional flow. They found that 
for a favorable pressure gradient the wall static pressure measured 
with the surface tube differed within 2% of the free stream dynamic 
head compared to the conventional method. 

It is well known that a static tube placed in the flow stream, 
away from the solid boundary, is insensitive to the yaw angle of ± 5®. 
However, beyond yaw angle of 5® the error in measured static pressure 
increases rapidly. The yaw characteristics of the surface tube in a two- 
dimensional flow without the presence of pressure gradient was ex- 
amined by Sproston and Goksel. The comparison between the static 
pressure measured with the surface tube and the wall tap showed 
small differences (< 0.5% of the free stream dynamic head) within the 
range of yaw angles (up to 10 degrees) investigated. Also the differ- 
ence in the two pressures was found to increase in the presence of 
pressure gradient. A detailed experimental investigation of the yaw 
characteristics of the surface tube in two and three-dimensional tur- 
bulent flows was performed in the present study. The surface static 
pressure measurements, on the flat plate centerline using plate 'B', 
were also obtained with the conventional method (i.e., with the wall 
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taps). The measurement locations of the surface taps in the flat plate 
centerline are shown in Fig. 3.2(a). 

3.2. 3.2 Three-Dimensional Flow 

The wall shear stress and the surface static pressure measure- 
ments were obtained upstream of the airfoil leading edge of the test- 
model at four streamwise stations viz., stations 1 to 4. The measure- 
ment points ahead of the airfoil leading edge are shown in Fig. 3.14 
(also see Fig.3. 12). It should be mentioned here, that all the mea- 
surements on the flat plate surface of the test model were made on 
the right half of the test model (looking in the downstream direction). 
On the flat plate the Preston tube was first aligned parallel to the test- 
model centerline with its axis of rotation parallel to Z-axis and then 
the probe was rotated for the maximum value of the pressure. Corre- 
sponding to this setting of the Preston tube, the reading of the pres- 
sure and the angular position of the Preston tube were recorded. This 
procedure was followed at each measuring grid point. Proper care was 
taken to make sure that the Preston tube tip rested on the flat plate at 
each measuring point. This was essential because the pressure mea- 
sured by the Preston tube at the wall is needed to evaluate the wall 
shear stress. The Preston tube tip contact with the flat plate surface 
was checked with an electrical circuit (Fig. 3.15) in addition to the vi- 
sual inspection. In the electrical circuit method the contact of the 
Preston tube tip was ensured once the non-zero voltage was recorded 
on the digital voltmeter. 

In the blade end-wall corner region, the Preston tube and the 
surface static pressure probe traverses were made at six streamwise 
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Fig. 3.14 Locations of Measurement Upstream of the Blade Leading Edge 


Test Model 
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Fig. 3.15 Sketch Showing Line Diagram of Electrical Circuit to Check 
the Preston Tube Tip Contact With the Surface 



stations, stations 5 to 10, The traverses In Y- and Z-directions ranged 
up to 76 mm from the test-model centerline. For the three-dimen- 
sional case the locations of measurement on the flat plate are shown in 
Fig. 3.16. In order to obtain adequate spatial resolution a large number 
of measuring points were selected near the comer line in the three- 
dimensional case as compared to the two-dimensional case. Figs. 
3.17(a) and 3.17(b) show the measurement locations on the blade 
surface. Also, for the three-dimensional case proper precautions were 
taken to make sure that the Preston tube tip rested on the flat plate 
and the airfoil surface at each measuring location. On the airfoil surface 
additional precautions were necessary because of the surface curva- 
ture. 

The Preston tube alignment procedure for the flat plate has al- 
ready been described. On the airfoil surface the Preston tube tip was 
first made parallel to the flat plate (axis of rotation parallel to Y-axis) 
and then it was made tangential to the airfoil surface at the point of 
measurement making sure the tip touched the airfoil surface. The 
Preston tube was then rotated for the maximum pressure and, for this 
setting of the Preston tube, data was recorded for the pressure and its 
angular position. 

In the blade end-wall comer region, detailed measurements of 
the wall static pressure were also obtained with the conventional tech- 
nique. The interchangeable instrument Disk No. 2 at three streamwise 
stations, stations 5, 7 and 9, was used to acquire the surface static 
pressure data. The disk was also rotated at each streamwise station so 
that the information about the surface static pressure over large surf- 
ace area could be obtained. 
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Test-Model in Corner Region 
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Test-Model in Corner Region 
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Fig. 3.17(b) Locations of Measurement on the Blade Surface of the 
Test-Model in Comer Region 


In the second case (Case B without the turbulence generating 
grids) wall shear stress data were collected on both the surfaces of the 
test-model in the blade end-wall comer region. The Preston tube and 
the pitot-static probe traverses were made at four axial stations, sta- 
tion 5 and stations 7 to 9. The measurement locations on both surfaces 
in this case are the same as in Case A (see Figs. 3.16 and 3.17). 

3.2.4 Calibration Equations of the Preston Tube and the Calculation 
Procedure 

The magnitude of the wall shear stress was calculated using 
Patel's[87] two-dimensional calibration equations. These equations 
were discussed in Chapter II. Patel used various sizes of Preston tubes 
ranging from an outer diameter of 0.60 mm to 12.65 mm. In the pre- 
sent study the outer diameter of the Preston tube was 0.81 mm and 
was within the range of Preston tubes used by Patel. The Preston tube 
and the static pressure probe data were converted into units of the 
wall shear stress using proper conversion factors. First the value of the 
parameter x* was calculated which is defined as, 

x* = logio (APd 2 /4pv 2 ) 

A suitable calibration equation was selected based on the magnitude of 
the parameter x*. and the value of y* was then calculated. Once the 
value of y* was known, the value of tw could be calculated from the 
relation, 

y* = log io (x w d 2 /4pv 2 ) 

The value of the free stream velocity (Uoo) was calculated from the 
recorded reference total pressure using the relation. 
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Uoo = [2(P t - Pa)ref/p) 0 - 5 

During the period of experimentation the atmospheric pressure and 
the temperature were recorded a number of times and their average 
values were used to calculate fluid properties (p and v). 

3.3 Wall Pressure Fluctuation and Wall Pressure-Velocity Correlation 

Measurements 

3.3.1 Wall Pressure Fluctuation Intensity 

The root-mean-square(RMS) value or the intensity of the wall 
pressure fluctuation was measured at four streamwise stations, stations 
1, 5, 7 and 9, using the flush-mounted piezoresistive type pressure 
transducers. At streamwise station 1, x/D = - 8, the wall pressure fluc- 
tuation intensity was obtained only on the stagnation streamline. At the 
other three streamwise stations, stations 5,7 and 9, the wall pressure 
fluctuation intensity was measured at selected locations (Fig. 3.18) on 
the flat plate surface of the test-model in the blade end-wall corner 
region. The output of the pressure transducer was directly fed to a 
Masscomp computer. The analog output was digitized at a 20 KHz 
sampling rate and 100,000 samples were taken. The excitation voltage 
to the pressure transducers was supplied with the differential DC am- 
plifier signal conditioning system (Model No. 7600A) supplied by Dy- 
namics Division of Waugh Controls. The excitation voltage was front 
panel controlled on the signal conditioning system. In the present 
work the excitation voltage was set at 15 V DC. The signal condition- 
ing system had a built-in 1 1 position filter switch. This switch had a 
capability of selecting a - 3 dB cutoff point from 1 Hz to a wideband 
position that made the - 3 dB frequency well above the highest desig- 
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All Dimensions sire in mm 
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Flat Plate Surface in Comer Region 


nated cutoff frequency. In the present study, the higher value of cutoff 
frequency used was 30 KHz. Table 3.5 presents specifications of the 
signal conditioning system as provided by the manufacturer. The am- 
plified data of the wall pressure fluctuation obtained from the signal 
conditioning system was directly fed to the A/D (analog to digital) unit 
of the Masscomp mini-computer. The A/D unit having 16 channels 
which, in combination with an S/H (sample and hold) unit, had a ca- 
pability of simultaneously sampling all the channels. The maximum 
throughput rate of the A/D with an S/H unit was 1 MHz. 

3.3.2 Power Spectra of the Fluctuating Wall Pressure 

The power spectral density function of the fluctuating wall pres- 
sure was obtained at different locations described earlier (see section 
3.3.1) using a Fast Fourier Transform (FFT) technique. This technique 
is commonly used because computational time and cost can be signifi- 
cantly reduced compared to the direct Fourier series method. The 
power spectral density function, <l>p w (f) at any frequency is related to 

the Fourier transform, Xpwff.T) by the relation! 137], 

<|> Pw (f) = (2/T) I Xpw(f.T) 1 2 (3.3) 

where T = N/f s is the finite time interval of the sample record x Pw (t) 

and N and fs represent the number of samples and sampling rate, re- 
spectively. The Fourier transform of the sample record, x p w W is given 

by- 

Xp w (f,T) = (l/f s ) I x p (n) exp[- i2xfn/f s ] (3.4) 
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TABLE 3.5 Performance Characteristics and Specifications of 
Differential DC Amplifier Signal Conditioning 
System as Supplied by Manufacturer (Waugh 
Controls Corporation) 


Parameter 
A Input 

1. Input Impedance 

2. Source Impedance 

3. Overscale Input 

4. Input Bias Current 


BL Tape Output 

1. Output Capability 

2. Capacitive Loading 

3. Output Impedance 

4. Output Level 

5. Output Zero Control 

(a) Adjustment Range 

(b) Resolution 

6. Zero Indicator 


Characteristics 

25 MQ shunted by 500 pF 
Specifications met with up to a 
source impedance of 1 KT2 (10 Kfl 
operation permitted) 

± 30 V DC, or peak AC without 
damage 

< 2 nA at 250C 

< 0.5 nA / OC 


± 10 V at 5 mA, limited to ± 

15 V at 25 mA. short-circuit 
protected 

up to 0. 1 |iF. no instability 
1.0 Q in series with 20 *iH 
Adjustable from 0 V to ± 10 V 

± 40 mv 
± 0.5 mv 

LED indicates the tape output zero 
status 

(a) Red if positive 

(b) Green if negative 

(c) off if within 1 mv of zero 
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TABLE 3.5 (Continued) 


c 


AC 


( All AC Specifications are 
independent of gain steps 
or Variable gain setting) 

1. Frequency Response 

2. Setting Time 

3. Overload Recovery Time 


± 1% (DC to 10 KHz) 

± 1 dB (10 KHz to 50 KHz) 

- 3 dB (above 100 KHz) 

< 25 ps to 0. 1% of final value 

< 50 ps to 5% of full scale for 


4. Slewing Rate 


5. Noise 


any overload up to 10 times full 
scale input not exceeding ± 20 V 
DC or peak AC 

3.7 V/ps, 20 V p-p output to 60 


KHz 





RTI 



Bandwidth 

5.0 pV 

rms 

0.1 

Hz 

to 

100 KHz 

3.0 pV 

rms 

0.1 

Hz 

to 

50 KHz 

2.0 pV 

rms 

0.1 

Hz 

to 

10 KHz 

4.0 pV 

P-P 

0.1 

Hz 

to 

100 KHZ 

1.0 pV 

p-p 

0.1 

Hz 

to 

10 KHz 


+ 200 pV rms RTO 


u dc 

1. Zero Drift (Constant 
Temperature) 

2. Temperature Coefficient 
of Drift 

3. Linearity J 

4. Gain Accuracy 


± 2 pV RTI, ± 200 pV RTO ± 200 
pV/vV(output offset) /8 hours after 
1/2 hr. warmup 

±0.5 nV/Oc RTI, ± 100 pV/Oc RTO 
±100 pV/V/Oc of output offset 
± 0.005% of full scale at DC 
± 0.1% in calibrate position of the 
CAL switch 
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TABLE 3.5 (Continued) 


5. Gain Stability ± 0.005% / 200 hours. 

± 0.005% / °C 

6. Zero Range ± 100 pV RTI of zero adjustment 

RTO set by output controls 


E. Common-Mode 

1. Rejection Ratio 

Gain 

1000 

100 

10 

1 

2. Operating Level 


3. Overscale 
4 Input Impedance 


Source Impedance 
0 ohm 350 Ohm 
126 dB 120 dB 

106 dB 100 dB 

86 dB 80 dB 

66 dB 66 dB 

± 50 V DC from DC to 60 Hz. 
Decrease at a rate of 6 dB/ Octave 
above 1 KHz. Limited to 1 V p-p up 
to 100 MHz 

± 75 V DC or peak AC without 
damage 

2000 MQ shunted by 1.5 pF 
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The Fourier components at discrete frequencies, fk are computed and 
related to Xp^f.T) as follows, 

Xpw(k) = Xp w (fj c ,T) fg = Z Xp^(n) cxp(- i27tkn/N] 

k = 0.1 N-l (3.5) 

where fk = kf s /N. The values of Fourier components are calculated at 
different frequencies, fk using Eqn. 3.5 in FFT technique. Further de- 
tails related to the FFT technique are given in reference! 137], In the 
present work, a standard FFT subroutine was used to calculate the val- 
ues of XpwfkHlSSl. Once the values of the Fourier components were 

known, the power spectral density functions at discrete frequencies 
were calculated using Eqn. 3.3 

3.3.3 Wall Pressure-Velocity Correlations 

The measurement of wall pressure-velocity correlations required 
information about the fluctuating velocity. The fluctuating velocity was 
measured with the single normal sensor hot wire probe. The correla- 
tion of the fluctuating wall pressure with the fluctuating velocity was 
measured at one streamwise station, station 1, upstream of the blade 
leading edge and three streamwise stations, stations 5,7 and 9, in the 
blade end-wall comer region. At each streamwise station in the comer 
region two locations were selected in the Y-direction where the pres- 
sure transducers were flush-mounted (Fig. 3.18). Corresponding to 
each Y-location the hot wire probe was traversed in the Z-direction 
covering a distance of approximately 30 mm. The wall pressure and 
velocity fluctuation data were simultaneously sampled at a 20 KHz 
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sampling rate. The number of samples recorded were 80,000 because 
of the space limitations of the floppy disk. 

3.3.4 Data Processing for Wall Pressure Fluctuation Intensity, Power 
Spectral Density Function, and Wall Pressure-Velocity 
Correlation 

The wall pressure and velocity fluctuation data were obtained in 
digitized form. The digitized data were first converted into units of 
volts using the conversion factor of the A/D unit. The A/D unit used 
had a 12 bit resolution. Therefore, for the voltage range selected, 
which was ± 5 volts, the conversion factor became 10/4096 volts/ 
digits. Once the digitized output was converted into volts further proc- 
essing depended on the variable of interest. The various parameters of 
interest which were calculated are briefly described in this section. 
3.3.4. 1 Wall pressure Fluctuation Intensity 

The digitized wall pressure output was converted into pressure 
units using an appropriate calibration constant (see Table 3.2) for the 
pressure transducer. The wall pressure fluctuation intensity was cal- 
culated from the relation, 

(p'w) 2 - I £ (Pwl 2 )/N - IpWl 2 (3.6) 

where p w = ( £ Pw)/N and Pw represent the mean and instantaneous 

t 

values of the wall pressure, respectively. In Eqn. 3.6, N represents the 
total number of samples. The value of p' w then gives the required wall 

pressure fluctuation intensity. The wall pressure fluctuation intensity 
was also calculated from the power spectrum. If 0p w (f) represents 

power spectral density function of wall pressure, then the total area 
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under the plot of power spectral density function versus frequency 
represents square of the wall pressure fluctuation intensity. 
Mathematically, 

(p'w) 2 * l “ fp w (f) df (3.7) 

The calculation of p' w using Eqn. 3.7 gave a check on its value ob- 
tained from Eqn. 3.6. 

3.3.4.2 Power Spectra of Wall Pressure Fluctuation 

A general procedure to obtain the power spectral density func- 
tion was described in section 3.3.2. The digitized wall pressure data 
were first converted into pressure units. The standard FFT subroutine 
was then used to obtain discrete Fourier components. Once the Four- 
ier components were known, then the power spectral density function 
at discrete frequencies was calculated using Eqn. 3.3. 

3.3.4.3 Wall Pressure-Velocity Correlation 

The digitized wall pressure data were converted into pressure 
units as described earlier. The digitized output in volts from the hot 
wire was converted into the instantaneous velocity using King’s law, 

E b 2 = Ai + BiU n (3.8) 

where E b and U are the instantaneous voltage output of hot wire and 
velocity, respectively. In Eqn. 3.8, Ai, Bi and n are constants which 
were obtained by calibration of the hot wire. The calibration procedure 
for the hot wire is given in appendix A. The mean component of ve- 
locity was calculated using the relation, 

U = ( IU)/N (3.9) 
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Once the mean value of wall pressure and velocity were known, the 
wall pressure-velocity correlation was calculated from the relation, 

Pw u = ( X Pw U)/N - Pw u (3.10) 

The correlation coefficient was calculated using the relation, 

Rpu = Pw u/p'w u’ (3. 1 1) 

where u' is the RMS value of the velocity fluctuation which was calcu- 
lated like p’ w described earlier. Computer programs were written to 

calculate the various parameters described above and are presented in 
appendix B. 


127 


CHAPTER IV 


RESULTS AND DISCUSSION ON COMBINED EFFECTS OF INITIAL 
ROUGHNESS AND FREE STREAM TURBULENCE ON THE 
TURBULENT BOUNDARY LAYER DEVELOPMENT 


The mean velocity profiles, the surface static pressure and the 
wall shear stress measurements were obtained for four cases de- 
scribed in chapter III (section 3.2) at seven streamwise stations on a 
flat plate. The data obtained for these four cases was then used to ana- 
lyze the influence of free stream turbulence on a turbulent boundary 
layer developing downstream of the distributed roughness strip in the 
vicinity of the plate leading edge. The results of the above mentioned 
investigation are presented and discussed here. 

4.1 Boundary Layer Development in the Absence of a Transition 

Promoting Device and Free Stream Turbulence (Case 1) 

A turbulent boundary layer developing with a natural transition 
and low level of free stream turbulence (usually present in wind-tun- 
nels) is considered ideal and is well documented in the literature. 
Therefore, this case was investigated first so that the other cases 
could be compared against this ideal case. For such an ideal case, sim- 
ilarity in mean velocity profiles at different streamwise locations is at- 
tained if the velocity defect non-dimensionalized by the local friction 
velocity, (U e -U)/u*, is plotted against z/8[56, 139), where U e and u* 

are local free stream velocity and friction velocity, respectively, and 8 
is the local boundary layer thickness. The mean velocity profiles at four 
streamwise locations (viz., at xl = 366.0 mm, 518.4 mm, 620.0 mm, 
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and 823.2 mm) on (U e -U)/u* and z/8 axes are shown in Fig. 4.1(a). 
Four streamwise locations sire chosen mainly to minimize overlapping 
of the data points. The velocity profiles are clearly observed to be self- 
similar except at x = 366 mm. This deviation from self- similarity can 
possibly be because the boundary layer has not yet attained its fully-de- 
veloped state. The mean velocity profiles in the inner law variables, at 
seven streamwise locations, are shown in Fig. 4.1(b). The values of fric- 
tion velocity (u*) used in Fig. 4.1, described above, were calculated 
from the local wall shear stress measured with the Preston tube. The 
solid lines in Fig. 4.1(b) represent the logarithmic law of the wall 
written as 

u+ = U/u* * 1/K In z+ + B (4. 1) 

where z + = zu*/v, and K and B are constants with the values as de- 
scribed earlier (section 3.2). The law of the wall, in general, is ob- 
served to exist in a region z + < 250 and results agree with Eqn. 4.1 
with a small change in the value of B. A larger shift from the solid line 
at xl = 366 mm, compared to that at the other measurement loca- 
tions, again suggests some deviation from its fully-developed state as 
was observed in Fig. 4.1(a). 

The skin friction coefficient (Cf) measured with the Preston tube 
and its comparison with the Cf values calculated from Clauser’s[96] 
technique, using the mean velocity profile data in the inner region, is 
shown in Fig. 4.2. The skin friction coefficient values were also calcu- 
lated by using empirical relations proposed by other investigators. 
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(a) Velocity-Defect Plot (b) Logarithmic-Law Plot 

(Note: Ordinate Scale Refers to Lowest Profile. 
Also Ordinate Scale Shifted by 2 Units 
for Second Profile Onwards) 
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Fig. 4.2 Variation of Cf With Ree for Case 1 



These empirical relations are 

Ludwieg-TillmanI 140): Cf = 0.246 Ree -0 - 268 10“0-678H ( 4 . 2 ) 

Coles[141]: Cf = 0.3 exp(-1.33H)/(logio Ree) 1 -? 4 + 0.31H (4.3) 

Robertson and Holt[63]: Cf = (4.4 + 3.8 log 10 Ree) -2 (4.4) 

The values of Cf calculated from Eqns. 4.2-4.4 are also shown in Fig. 
4.2. The measured Cf values compare well (maximum difference < 

1.5%) with the values obtained from Eqn. 4.2 and Eqn. 4.4. The values 
of Cf calculated using Eqn. 4.3 are consistently smaller (= 4%). The Cf 

values obtained with Clauser’s method are consistently larger (=* 5%). 

The variation of the boundary layer integral parameters 8, 5*/8 
and 5/0 with streamwise distance is shown in Fig. 4.3. The values of 8 

calculated from Eqn.3.1 are found to be in good agreement with the 
values obtained from the velocity profiles corresponding to U/Ue = 
0.995. Clauser[139] noted, that the usual boundary layer thickness (8) 
was not a well defined thickness. Therefore, the boundary layer thick- 
ness should be defined such that it could be determined mainly from 
the parameters of the universal plot. He further noted, that the inte- 
gral expression defined as 

A = 5f° (U e - U)/u* d(z/8) (4.5) 

0 

showed the required property. The integral parameter A is known as 
Clauser's integral thickness. For a zero pressure gradient turbulent 
boundary layer, the value of A/5 was reported to be 3.6(139]. In the 
present study the average calculated value of A/8 is 3.8. In Fig. 4.3 
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other cases are also shown for comparison and will be discussed later. 
The boundary layer data for four cases investigated is presented in 
Table 4.1. 

Fig. 4.4 shows the variation of Clauser's shape parameter (G), 
shape factor (H = 8*/0), and Coles's wake profile parameter (II). The 
Clauser's shape parameter is defined as 

G = H(U - U e )/u*] 2 d(z/A) = (1 - 1/H) (2/Cf )0-5 (4.6) 

O 

The average value of G = 7.15 (see Fig. 4.4(a)) is obtained in the pre- 
sent investigation. Clauser(139] reported a value of G = 6.8 for a con- 
stant pressure turbulent boundary layer. In Fig. 4.4 other cases in- 
vestigated are also shown for the comparison and will be discussed 
later. 

The variation of the boundary layer integral parameters for the 
four cases, with Reynolds number based on momentum thickness, is 
shown in Fig. 4.5. In Fig. 4.5(c) values of 5/0 obtained from an empiri- 
cal relation reported by Robertson and Holt(63], written as 

5/0 = 8.22 + 0.475 log 10 Ree (4.7) 

are also included. The values of 8/0 obtained in the present study were 

consistently smaller (by maximum of 10%) than the values predicted 
from Eqn. 4.7. The variation with Reg of Clauser’s shape parameter 

(G), shape factor (H), and Coles’s wake profile parameter (TT) for the 
four cases is shown in Fig. 4.6. The value of H was also calculated from 
an empirical relation reported by Robertson and Holt[63) 

H = 1.09 + 0.40/(0.72 logio R e0 - 1.0) (4.8) 


134 



TABLE 4. 1 Details of Mean Flow Parameters for Four Cases 


Case 1 : Natural Transition, No Grids 


xl 

(mm) 

8 

(mm) 

e 

(mm) 

Cfx 10 3 

H 

366.0 

■gm 

1.243 

3.536 

1.396 

518.4 


1.502 

3.339 

1.410 

569.2 

■tn 

1.569 

3.301 

1.412 

620.0 

14.6 

1.621 

3.254 

1.414 

670.8 

15.0 

1.703 

3.191 

1.414 

747.0 

16.2 

1.811 

3.173 

1.407 

823.2 

17.5 

1.893 

3.173 

1.398 


Case 2 : Natural Transition With Grids 


xl 

(mm) 

TUoo 

8 

(mm) 

e 

(mm) 

Cfx 103 

H 

366.0 

0.0172 

13.2 

1.368 

3.626 

MMM 

518.4 

0.0155 

18.4 

1.621 

3.473 

i ififl 

569.2 

0.0153 

19.0 

1.745 

3.443 

WmSm 

620.0 

0.0149 

19.9 

1.788 

3.436 

■fff 

670.8 

0.0144 

21.0 

1.888 

3.362 


747.0 

0.0139 

21.3 

1.947 

3.340 

mtm 

823.2 

0.0136 

22.2 

2.038 

3.323 

m 
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TABLE 4. 1 (Continued) 


Case 3 : Initial Roughness, No Grids 


xl 

(mm) 

6 

(mm) 

e 

(mm) 

Cfx 10^ 

H 

366.0 

31.8 

3.026 

3.369 

1.285 

518.4 

29.0 

2.616 

3.347 

1.299 

569.2 

29.8 

2.693 

3.323 

1.301 

620.0 

31.4 

2.773 

3.298 

1.300 

670.8 

31.7 

2.824 

3.271 

1.300 

747.0 

32.2 

2.895 

3.221 

1.303 

823.2 

33.8 

2.944 

3.269 

1.303 


Case 4 : Initial Roughness and Grids 


xl 

(mm) 

Tuoo 

S 

(mm) 

e 

(mm) 

Cfx 103 

H 

366.0 

0.0172 

26.7 

2.887 

3.347 

1.316 

518.4 

0.0155 

28.2 

2.817 

3.207 

1.313 

569.2 

0.0153 

28.4 

2.869 

3.183 

1.315 

620.0 

0.0149 

28.9 

2.951 

3.184 

1.315 

670.8 

0.0144 

29.5 

2.984 

3.169 

1.314 

747.0 

0.0139 

30.9 

3.049 

3.140 

1.314 

823.2 

0.0136 

32.0 

3.077 

3.196 

1.311 
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The values of H obtained from Eqn. 4.8 are shown as a solid line In Fig. 
4.6(b). In general, Eqn. 4.8 predicts lower values (= 3%) compared to 
the values obtained in the present investigation. In Fig. 4.6(c) the solid 
line corresponds to Coles's empirical relation 

n = 0.55[1 - exp(- 0.243^0 5 . 0.298^)1 (4.9) 

where ^ = Ree/425 - 1.0 

The values of n obtained in the present study are generally higher 
compared to those predicted from Eqn. 4.9. 

4.2 Turbulent Boundary Layer Development With Natural Transition 
and the Turbulence Generating Grids (Case 2) 

Fig. 4.7(a) shows mean velocity profiles at different stream wise 

stations, in the presence of increased free stream turbulence, on (U e - 

U)/u* and z/8, axes. The velocity profiles are observed to be self-simi- 
lar except at xl = 366 mm. In this case also, as in case 1, the boundary 
layer appears to have not reached its fully-developed state at xl = 366 
mm. The mean velocity profiles, in inner law variables at different 
streamwise locations, are shown in Fig. 4.7(b). The region in which 
law of the wall exists has not changed with the increased free stream 
turbulence level compared to case 1. 

The comparison of mean velocity profiles non-dimensionalized 
by local free stream velocity (Ue). for four cases at four streamwise 
locations, is shown in Fig. 4.8. The other cases ( viz., cases 3 and 4) 
are also shown in Fig. 4.8 for comparison purposes and will be dis- 
cussed later (see sections 4.3 & 4.4). The mean velocity profiles are 
observed to be fuller compared to case 1 at each streamwise location. 
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Fig. 4.7 Development of Mean Velocity Profiles for Case 2 
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The fullness of velocity profile implies larger velocity gradient at the 
wall and, therefore, increased wall shear stress. This trend of incre- 
ased wall shear stress due to the fullness of velocity profiles is in good 
agreement with the measured values of wall shear stress (see Fig. 
4.12). 

The comparison of velocity-defect profiles for four cases at four 
streamwise stations is shown in Fig. 4.9. Cases 3 and 4 will be dis- 
cussed in sections 4.3 and 4.4, respectively. The influence of in- 
creased free stream turbulence is to decrease the velocity-defect. Sim- 
ilar observations i.e., fullness of velocity profiles and decrease in vel- 
ocity-defect because of the increase in free stream turbulence, were 
reported earlier by Evans[60]. 

Fig. 4.10 shows the comparison of velocity profiles, for four cases 
at four streamwise locations, in inner law variables. In case 2, the wake 
component is observed to decrease as compared to case 1 which is in 
agreement with the calculated value of the wake profile parameter, n 
(see Fig. 4.4(c)). Earlier studies by other investigators for example, 
Evans[60] and Raghunathan & McAdam[68) have reported similar in- 
fluences of increased level of free stream turbulence. The decrease in 
wake component and velocity-defect can be attributed to the incre- 
ased mixing in the outer region due to the increased level of the free 
stream turbulence. 

The comparison of Cf values at various streamwise locations by 
different methods is shown in Fig.4. 11. The Ludwieg-Tillman relation 
(Eqn. 4.2) predicts closer values (within 3%), whereas the differences 
with Clauser’s method are within 4% to 6%. The comparison of Cf val- 
ues at the same streamwise locations, measured with the Preston tube 
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Fig. 4.9 (Contd.) Comparison of Mean Velocity Profiles on 
Velocity- Defect Plot 
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Fig. 4. 1 1 Variation of Cf With Reg for Case 2 


for four cases, is shown in Fig.4. 12(a). The comparison of Cf with Reg 

for the four cases is shown in Fig. 4.12(b). The increase in the value of 
Cf because of the increased free stream turbulence level as compared 

to case 1 is *• 5%. Bradshaw[142) proposed an empirical relation (rep- 
orted by Raghunathan & McAdam[68]) to evaluate increase in Cf value 

due to the increase in free stream turbulence as 

ACf/CfO = 3.2 Tuoo (4.10) 

where Cfo is the value corresponding to the zero turbulence level case. 
Cf values obtained in the present study are found to be in good agree- 
ment with the values predicted by Eqn. 4.10. The recent studies of 
Raghunathan & McAdam[68] and Hancock & Bradshaw(69) have shown 
a non-linear effect on Cf at low levels of free stream turbulence. For a 

free stream turbulence in the range 0 < Tu<» < 0.04, Raghunathan & 
McAdam[68] proposed an empirical relation (Eqn. 2.16) for evaluating 
change in Cf. 

The Cf values predicted by Eqn. (2.16) are smaller than the 
measured values. No significant changes in ACf/Cfo at different stream- 

wise locations are observed implying an independence of the change 
in Cf on Reg with the increase in free stream turbulence level. Similar 

observations were reported by other investigators[68,69]. 

At the same streamwise location, the increase in 8 for case 2 
compared to case 1 in the present study is approximately 35% (see 
Fig. 4.3(a)). The increase in 8 observed at the same streamwise loca- 
tion can be due to the shift in the origin of the boundary layer. It is 
well established in the literature that the increase in free stream tur- 
bulence level leads to early transition and an upstream shift of the ori- 
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(a) Variation of Cf as a Function of Streamwise Distance 
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Fig. 4.12 Variation and Comparison of Skin Friction Coefficient 
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gin of the turbulent boundary layer. Therefore, it would be more ap- 
propriate to compare the variation of 8 with Reg in the two cases. It is 

observed that at the same value of Re0 the value of 6 is higher in case 2 

compared to its value in case 1 (Fig. 4.5(a)). At the same streamwise 
distance, the non-dimensional value of the displacement thickness 
(5*/ 5) is observed to decrease with the introduction of free stream 
turbulence (Fig. 4.3(b)). The decrease in 5*/8 on the average is * 22%. 
The main reason for this decrease in 8*/ 8 is because 8 increases much 
more than 8* as the turbulence level is increased. The value of 8/0 in- 
creased in case 2 compared to case 1 (Fig. 4.3(c)). The average in- 
crease in 8/0 is found to be » 22%. The comparison of results obtained 
in the present investigation with the values of 8/0 obtained from em- 
pirical relation, Eqn. 2.12, proposed by Raghunathan and McAdam[68] 
suggests that Eqn. 2.12 underestimates the changes in 8/0 because of 
the increase in the free stream turbulence level. 

The value of G is observed to decrease due to the presence of 
free stream turbulence, the average decrease being = 13% (Fig. 4.4(a)). 
The results obtained in the present study compare well with the val- 
ues obtained from an empirical relation reported by Raghunathan and 
McAdam[68] 

AG/G 0 = -1/3 (u'oo/u*) (4.11) 

where Gq corresponds to the zero free stream turbulence case. The 

average value of G = 6.2 (less than the value of 7.15 of case 1) again 
suggests that the mean velocity profiles are far from the self- similar 
state of case 1. 
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The values of H decreased compared to case 1 . The decrease In 
the value of H is * 5%. The value of H was also calculated from an em- 
pirical relation (Eqn. 2.14) proposed by Raghunathan and McAdam. 
The values obtained in the present measurements are little higher 
compared to those predicted by Eqn. (2.14). 

Fig. 4.4(c) shows the effect of increased free stream turbulence 
on n at the same streamwise locations. The decrease in the value of n 
compared to case 1 is in agreement with the decrease in wake com- 
ponent observed in Fig. 4.10. The value of II was also calculated using 
empirical relation (Eqn. 2.15) proposed by Raghunathan and McAdam. 
Present measurements suggest that Eqn. 2.15 underestimates the ef- 
fect of Tuoo on changes in n. 

To summarize, the effect of increased free stream turbulence on 
various mean flow parameters, it can be stated that the trends in the 
variation of different mean flow boundary layer parameters is similar to 
those reported by other investigators (see Table 2.1). 

4.3. Influence of Initial Roughness on the Turbulent Boundary Layer 

Development Without the Turbulence Generating Grids (Case 3) 

In this section the results of a boundary layer developing down- 
stream an initially distributed roughness strip and its comparison with 
case 1 are discussed. 

Fig. 4.13(a) shows velocity profiles at four streamwise locations 
on a velocity-defect plot. The self-similarity of velocity profiles is quite 
evident except at the first station (located = 300 mm from the end of 
roughness strip). The velocity profiles in inner wall variables are 
shown in Fig. 4.13(b) for seven streamwise locations. The law of the 
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wall region is observed to extend over a larger distance (z + < 600) 
compared to case 1. The velocity profiles in the presence of initial 
roughness are fuller compared to the velocity profiles in case 1 (see 
Fig. 4.8). The fullness of the velocity profile implies a higher wall shear 
stress which is in good agreement with the measured Cf values except 
at xl = 366 mm (see Fig. 4.12(a)). At the first station the measured Cf 
value is smaller compared to the Cf value in case 1 and, therefore, is in 
contradiction with the observed fuller profile. A comparison of Cf val- 
ues for case 1 and case 3 shows lower Cf values close to the roughness 
strip. Further downstream Cf values are slightly higher compared to 
case 1 (see Fig. 4.12(a)). It is observed that a distance of = 40 8 q (5o is 

the boundary layer thickness at xl = 366 mm for case 1) from the 
roughness strip is needed for Cf values to reach the values obtained in 
easel. A comparison of Cf values by different techniques is shown in 
Fig. 4.14. The Cf values calculated from Ludwieg-Tillman relation differ 
by < 2% from the measured values. Smaller differences {= 3%) in Cf 
values, obtained by Clauser's technique and measured values, are ob- 
served compared to those observed in case 1 . 

The variation of Cf with Ree for case 1 and case 3 is shown in 
Fig. 4.12(b). Although, similar values of Re 0 could not be achieved in 
the two cases, but it was quite clear that at the same Reg, the values of 
Cf in case 3 are larger than for the case 1. The data obtained by 
Purtell, Klebanoff and Buckley! 143] and Huffman, Zimmerman and 
Bennet[62] are also included in Fig. 4.12(b). These two data sets were 
selected as their study also used a roughness strip to artificially 
thicken the boundary layer. The values of Cf in the present investiga- 
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Measured 



Fig. 4.14 Variation of Cf With Re© for Case 3 





tion are higher compared to these two data sets at similar values of 
Ree. 

A comparison of the velocity profiles for cases 1 and 3, on the 
velocity-defect plot shows a decrease in the velocity-defect up to the 
last measurement location (see Fig. 4.9). Although, the self-similarity is 
observed at different streamwise locations (see Fig. 4.13(a)), the mean 
flow has not reached the self-similarity distribution of case 1 up to the 
last measurement location which is located at a distance of - 65 So 

from the end of roughness strip. A comparison of the velocity profiles 
in the presence of initial roughness with case 1 on a log-law plot (see 
Fig. 4.10) shows a decrease in the wake component Implying a lower 
wake profile parameter, n, which agrees with the calculated value of n 
(see Fig. 4.4(c)). 

The boundary layer thickness for case 3 increases considerably 
(* 100%) compared to the value for case 1 (see Fig. 4.3(a)). A higher 
value of 8 is observed at the first station. A similar trend in the varia- 
tion of 8 was obtained from Coles' method (Eqn. 3.1). 

The value of 5*/8 for case 3 is lower (= 25%) than the value in 
case 1 (Fig. 4.3(b)). This is mainly because 8 for case 3 is much higher 
than the value of 8 in case 1. Moreover, the decrease in 8*/ 8 is ob- 
served to be independent of Re©. On the contrary, 8/0 is observed to 

increase (« 23%) compared to the value of case 1 (see Fig. 4.3(c)). The 
increase in 8/0 is also observed to be independent of Ree- 

A considerable decrease (= 20%) in the value of G for case 3 
compared to the value for case 1 is observed (Fig. 4.4(a)). The average 
value of G is = 5.7, which is different from the case 1 value (= 7.15), 
also suggests that the flow has not reached the self-similar state of 
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case 1. The values of H for case 3 decrease (~ 8%) in the presence of 
initial roughness compared to case 1 values (Fig. 4.4(b)). 

4.4 Combined Influence of Initial Roughness and Free Stream 
Turbulence (Case 4) 

The self- similarity of mean velocity profiles, on a velocity- defect 
plot, at different streamwise locations is clearly observed in Fig. 
4.15(a). At the first streamwise station in this case also, as in previ- 
ously discussed cases, some deviation from similarity law is observed. 
The presence of initial roughness in combination with free stream 
turbulence seems to have a little influence on the law of the wall (Fig. 
4.15(b)). 

In comparison to case 1, the mean velocity profiles are observed 
to be fuller at all streamwise locations implying larger wall shear stress 
(Fig. 4.8). On the contrary, the measured Cf values show lower values 

compared to case 1 except at the last measurement location (see Fig. 
4.12(a)). A comparison of Cf values for two cases, i.e., case 4 with case 

1, suggests the need for a larger distance downstream of roughness 
strip = 65 5o, compared to case 3, to reach a smooth plate value of wall 

shear stress. The above mentioned discrepancy of fuller velocity pro- 
files but decreased Cf values can be resolved if the comparison of case 

4 is made with case 3. In comparison to case 3, mean velocity profiles 
for case 4 are observed to be less full implying a lower wall shear 
stress which is in agreement with the measured Cf values (see Figs. 4.8 

and Fig. 4.12(a)). 

A comparison with case 1 of mean velocity profiles on velocity- 
deflect plot shows a decrease in velocity-defect and the decrease is 
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observed to be smaller compared to case 3 (Fig. 4.9). The observed 
change in the velocity-defect law (Fig. 4.9) in combination with the 
existence of self similarity (Fig. 4.15) implies that also in this case 
(case 4) the mean flow is far from achieving the self-similar state of 
case 1. 

At the same streamwise station, in comparison to case 1, the 
wake component is observed to be smaller and agrees with the calcu- 
lated value of n (see Fig. 4.10 and Fig. 4.4(c)). However, this decrease 
in wake component is smaller compared to the decrease observed in 
case 3. This observation is also confirmed by the variation of fl (Fig. 
4.4(c)). Based on the above mentioned discussion on mean velocity 
profiles it appears that the influence of increased free stream turbu- 
lence level on a turbulent boundary layer with and without the pres- 
ence of initial roughness have opposite effects. This observation will 
find more support once the changes in boundary layer integral param- 
eters are discussed. 

The value of 8 for case 4, even though, increased compared to 
case 1 but a decrease in 8 is observed compared to the value for case 3 
(Fig. 4.3(a)). In comparison to case 1, 8*/8 values for case 4 are smaller 
by = 16% (Fig. 4.3(b)). This may be attributed to the increase in the 
value of 8 compared to the value for case 1. However, an increase (= 
10-16%) in 8*/8 values are noticed compared to the values for case 3. 
The observed increase in the values of 8*/8 is because of the decrease 
in the values of 8 and the increase in the values of 8* compared to the 
values for case 3. The values of 8/0 increased compared to the values 
for case 1, with the increase being smaller than that observed for case 
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3 (Fig. 4.3(c)). A decrease of « 9% in 8/0 values was observed com- 
pared to the values of case 3. 

A comparison of the boundary layer thickness, shape factor and 
skin friction coefficient for the four cases at the last streamwise sta- 
tion, xl = 823.2 mm, is presented in Table 4.2. It is clearly observed 
that in the presence of initial roughness, the increase in free stream 
turbulence level has an opposing effect on the above mentioned pa- 
rameters compared to the case without the initial roughness. 

The average value of G calculated was = 6.0 and is far from the 
value obtained for case 1 (Fig. 4.4(a)). Moreover, the value of G is larger 
(by = 5%) compared to its value of case 3. The shape factor, H, varia- 
tion shows a decrease of = 7% compared to the values for case l(Fig. 
4.4(b)). However, a smaller increase in the values of H was obtained 
compared to the values in case 3. A comparison of measured Cf values 

for case 4 with other techniques is shown in Fig. 4.16. The differences 
observed are larger with Clauser's method than with the Ludwieg- 
Tillman law. 

A comparison of the influence of free stream turbulence on diff- 
erent boundary layer parameters, with and without an initial roughness 
is presented in Table 4.3. This comparison clearly shows that in the 
presence of initial roughness, the influence of free stream turbulence 
on a turbulent boundary layer has an opposite effect compared to the 
case without the initial roughness. The numerical values given in Table 
4.3 are the average values taken over all the streamwise stations inves- 
tigated. One of the plausible reasons can be the nature in which local 
wall shear stress changes. This is because in a turbulent boundary layer 
various boundary layer parameters are dependent on Cf. 
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TABLE 4.3 Influence of Free Stream Tiubulence on Mean Flow 
Boundary Layer Parameters With and Without 
the Initial Roughness 



Without Initial Roughness 

With Initial Roughness 

Parameter 

Qualitative 

% variation 
in parameter 

Qualitative 

% variation 
in parameter 

Mean velocity 
profile fullness 

More full 


Less full 


Defect 

velocity 

Decreases 


Increases 


Wake 

component 

Decreases 


Increases 


6 

Increases 

35.0 

Decreases 

- (3.0 - 8.0) 

8*/5 

Decreases 

- 22.0 

Increases 

10.0 - 16.0 

5/6 

Increases 

22.0 

i 

Decreases 

- 9.0 

G 

Decreases 

-13.0 

Increases 

3.0 - 6.0 

H 

Decreases 

- 5.0 

Increases 

< 2.0 

Cf 

Increases 

5.0 

Decreases 

- (2.0 - 4.0) 

n 

Decreases 

- (30.0 - 40.0) 

Increases 

20.0 - 70.0 
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CHAPTER V 


RESULTS AND DISCUSSION ON THE SURFACE STATIC PRESSURE, 
THE WALL SHEAR STRESS, AND THE MEAN VELOCITY PROFILES IN 
A BLADE END-WALL CORNER REGION 


Experimental results on the surface static pressure and the local 
wall shear stress vector in an upstream region and in the blade end- 
wall corner region are presented and discussed in this chapter. The 
measurements of mean velocity profiles obtained only in an upstream 
region of the blade leading edge are also described here. Yaw sensiti- 
vities of the static tube resting on the surface and the Preston tube 
were experimentally investigated in two- and three-dimensional flows. 
Measurement of various flow parameters, as mentioned above, were 
also made on the flat plate centerline without the presence of the 
blade (i.e., two-dimensional flow). The results of measurements for two 
-dimensional flow are also discussed in this chapter. As described in 
Chapter III, two flow conditions were investigated: the first case (Case 
A) is without the presence of separation bubble, and in the second 
case (Case B) a separation bubble existed in the blade end-wall corner 
region. Therefore, the results of two cases, Case A and Case B, are pre- 
sented in section A and section B, respectively. 


SECTION A 

5.1 Two-Dimensional Flow (Case A) 

There were two main reasons, as mentioned earlier in Chapter 
III, to obtain measurements of various flow parameters of interest in a 
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two-dimensional flow. The first reason was to check the proper 
working of the different probes used in the investigation. This was 
achieved by comparing the results obtained in the present investiga- 
tion for a two-dimensional flow case with similar results available in 
the literature. Secondly, these results are also used to compare stress 
values obtained in the blade end-wall corner region. 

5.1.1 Surface Static Pressure and Mean Velocity Profiles 

The surface static pressure coefficient, Cp, variation on the flat 
plate centerline in the streamwise direction is shown in Fig. 5.1. In 
the streamwise direction, a small decrease (approximately 1.5% of the 
free stream dynamic pressure) in Cp value is observed implying an ex- 
istence of favorable pressure gradient. The existence of favorable pres- 
sure gradient in streamwise direction is also evident from the free 
stream mean velocity measured at different streamwise stations (Fig. 
5.2). This observed decrease in the surface static pressure or increase 
in the free stream mean velocity in the streamwise direction can be 
attributed to the boundary layer growth on the sidewalls of the wind 
tunnel test-section. In Fig. 5.1, the surface static pressure measured 
using the conventional (wall tap) technique is also included. The Cp 
values obtained with the two techniques agree within one percent of 
the free stream dynamic pressure. 

It needs to be mentioned here, that the flow conditions for the 
case under discussion were similar to Case 4 discussed in detail in 
Chapter IV. The detailed description on the mean velocity profiles are 
presented in section 4.4. 
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5.1.2 Wall Shear Stress 

The wall shear stress measured with the Preston tube at differ- 
ent streamwise stations on the flat plate centerline (two-dimensional 
case) is presented in Table 5.1. The values of the wall shear stress at 
different streamwise stations were also obtained from the measured 
mean velocity profiles using Clauser’s[961 method and the empirical 
relation proposed by Ludwieg and Tillman[140J. The wall shear stress 
values obtained with the Preston tube, at different streamwise stations, 
are observed to agree within 1.5% as compared to Ludwieg-Tillman 
relation. However, Cf values obtained by Clauser's technique are, in 

general, found to be higher (within 5%) compared to the values ob- 
tained with the Preston tube. These larger differences, as compared to 
the latter technique, are probably due to changes in the values of con- 
stants used in the logarithmic law of the wall. 

5.2 Three-Dimensional Flow (Case A) 

The results of the yaw sensitivities of the static tube resting on 
the surface and the Preston tube are discussed in detail in appendix C. 
However, briefly the main observations are presented here. The re- 
sults of the surface static pressure and the wall shear stress mea- 
surements upstream and in the blade end-wall comer region are pre- 
sented in this section. 

5.2.1 Taw Sensitivities of the Surface Static Tube and the Preston 

Tube 

The yaw characteristics of the surface static tube showed that 
the smaller value (approximately 20%) of the surface static pressure 
would be obtained if the static tube was not aligned in the local flow 
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TABLE 5.1 Comparison of Cf Values by Different Methods 
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direction (see Fig. C2 of appendix C). Also the surface static tube was 
found to be insensitive to the yaw angle within ± 5^. The experimental 
results on the yaw characteristics of the Preston tube showed that it 
was insensitive to the yaw angle of ± 5®. It was also observed that if the 
Preston tube was not aligned in the local flow direction, it would have 
recorded approximately 15% less total pressure (see Fig. C3 of ap- 
pendix C). For the wall shear stress measurement, this observed error 
in the Preston tube output due to misalignment would lead to ap- 
proximately 10% error in the calculated value of the wall shear stress. 

5.2.2 Surface Static Pressure in Upstream Region of the Blade Leading 

Edge 

The surface static pressure coefficient variation in the Y-direc- 
tion upstream of the blade leading edge at four streamwise stations is 
shown in Fig. 5.3. All linear distances are non-dimensionalized by the 
blade leading edge diameter. The available literature, as presented in 
Chapter II, clearly shows that the blade leading edge shape signifi- 
cantly influences the flow characteristics in an upstream region and in 
the blade end-wall corner region. The characteristic length scale, 
therefore, can be the maximum thickness of the blade leading edge. In 
the present study, because the blade leading edge shape is circular, 
the diameter of the blade leading edge is considered a suitable choice 
for non-dimensionalizing various linear distances. In an upstream re- 
gion of the blade leading edge, the surface static pressure coefficient, 
Cp, decreases in the Y-direction and the rate of decrease is observed 
to be maximum at x/D = - 2.0 in the vicinity of the blade leading edge. 
At a distance of y/D = 6, the value of Cp is found to approach approxi- 
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Fig. 5.3 Surface Static Pressure Distribution Upstream of the 
Blade Leading Edge 
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♦ Conventional 


Fig. 5.4 Surface Static Pressure Variation on Stagnation Streamline- A 
Comparison Between Surface Static Tube and Conventional 
Techniques 






mately (within 2% of the free stream dynamic pressure) to its value far 
upstream, i.e., at x/D = - 20, of the blade leading edge. 

The variation of Cp on the stagnation streamline, in an upstream 
region of the blade leading edge, is shown in Fig. 5.4. The measured 
value of Cp obtained with conventional wall taps is also included in Fig. 
5.4. The maximum difference in measured Cp values by two tech- 
niques is found to be less than 2% of the free stream dynamic pres- 
sure. 

5.2.3 Wall Shear Stress in Upstream Region of the Blade Leading Edge 

The variation of the local wall shear stress, normalized by the far 
upstream value (at x/D = - 20, y/D = 0), at four streamwise stations 
upstream of the blade leading edge is shown in Fig. 5.5. This 
presentation of the local wall shear stress rather than the skin friction 
coefficient, is chosen for two main reasons. Firstly, the local free 
stream velocity is expected to change at different streamwise stations. 
Moreover, in the blade end-wall corner region, it is difficult to define 
the local free stream velocity because of the interacting nature of the 
shear layers. Secondly, this presentation clearly shows how the local 
wall shear stress, in the blade end-wall corner region, varies relative 
to the far upstream value at x/D = - 20. It is important to mention 
here, that the local wall shear stress at x/D = - 20 is approximately 
equal (within 3%) to the two-dimensional value at the same stream- 
wise location. The normalized local wall shear stress decreases as the 
blade leading edge is approached. A maximum decrease in wall shear 
stress of * 40%, from the value far upstream of the blade leading edge, 
occurs on the stagnation streamline where the adverse pressure gra- 
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Fig. 5.6 Variation of Surface Static Pressure Along Different Streamlines 
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dient Is also found to be maximum (see Fig. 5.3). The wall shear stress 
values, at different streamwise stations, reach the two-dimensional 
value within a distance of - 6D in the Y-direction. The wall shear 
stress values obtained with the floating element balance by 
McAllister! 1061 in an upstream region of his test-model at approxi- 
mately the same x/D distance are also presented in Fig. 5.5. The test- 
model used in McAllister's study consisted of a cylindrical body 
mounted on a flat plate. The body had a large circular leading edge 
(12.7 cm diameter) and was tapered in streamwise direction. It should 
be mentioned that the blade leading edge diameter of the present 
study is 1/10 that of the McAllister's study. It is quite evident that the 
trend in the variation of the local wall shear stress in Y-direction ob- 
tained in McAllister's measurements is similar to that obtained in the 
present study. However, the magnitude of the normalized wall shear 
stress, at different y/D distances, is relatively higher in McAllister’s 
study than in the present study. The possible reason for these ob- 
served differences of the order of 25% in the two cases can be because 
of the pressure gradient effects. 

5.2.4 Mean Velocity Profiles Upstream of the Blade End-Wall Comer 

Region 

The mean velocity profiles were measured on the flat plate sur- 
face upstream of the blade leading edge on five streamlines (see Fig. 
3.12 for the measurement locations). The traces of these five stream- 
lines were obtained from flow visualization using an oil dot matrix 
method. The local wall shear stress and the surface static pressure 
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were also measured at different measurement locations. The results of 
the above mentioned measurements are discussed in this section. 

The variation of surface static pressure coefficient along five 
streamlines is shown in Fig. 5.6. It is clearly observed that the maxi- 
mum adverse pressure gradient in the streamwise direction exists on 
the stagnation streamline (streamline no.l) where the value of Cp 
changed from 0.05 to 0.24. Moreover, away from the stagnation 
streamline i.e., as the distance y increases, the adverse pressure gra- 
dient in the streamwise direction decreases. In fact, the value of Cp on 
streamline no. 5, does not change significantly except between 
streamwise stations 3 and 4, where the value of Cp changed by 1% of 
the free stream dynamic pressure. 

Fig. 5.7(a) shows the comparison of non-dimensional mean ve- 
locity profiles at x/D = - 20 for two- and three-dimensioned flow cases. 
The mean velocity profiles on the logarithmic inner law of the wall 
plot for two- and three-dimensioned flows at x/D = - 20 is shown in 
Fig. 5.7(b). The solid line in Fig. 5.7(b) represents 

u+ sU/u* = (1/K) In z+ +B (5.1) 

where the values of the universal constants, K eind B, used are the 
same as described in Chapter IV. The value of u* used in Fig. 5.7(b) 
was calculated from the local wall shear stress measured with the 
Preston tube. The mean velocity profiles are observed to be self-similar 
in the outer-layer variables (Ue and 8*). However, on the logarithmic 
law plot, deviation from the solid line for the three-dimensional flow 
case is observed. 
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Fig. 5.7 Variation of Mean Velocity Profiles at Streamwise Location of x/D = - 20 



Figs. 5.8 - 5.12 show mean velocity profiles in the outer-layer 
and inner-layer variables on 5 streamlines. Along each streamline, at 
different streamwise stations, the mean velocity profiles show self- 
similarity (Figs. 5.8(a) - 5.12(a)). However, the inner region of the 
boundary layer on various streamlines is influenced differently. If the 
mean velocity profile is represented as a combination of the law of wall 
and the law of the wake, as proposed by Coles! 133], and is written as 

u+ a U/u * = (1/K) In z+ + B + (n/K) W[z/8] (5.2) 

where K and B are universal constants with the values given in Chapter 
IV. n is the wake strength parameter, and W[z/5] is the wake function 
with the condition that at z = 5, W[z/8] = 2. Along streamline no.l, the 
values of constants K and B change. Also, the value of wake strength 
parameter increases as the blade leading edge is approached. This 
trend of an increase in the value of n in the streamwise direction is 
observed up to streamline no. 3. However, on streamlines no. 4 and 5, 
the variation of n in streamwise direction changes. In fact on stream- 
line no. 5, the value of n decreases as the blade leading edge is ap- 
proached. Hsing and Teng[35] during their investigation of the up- 
stream region of a wing-body junction also reported a variation in the 
values of K, B, and n. They noted that the observed variation in the 
values of K, B, and n could be attributed to the changes in adverse 
pressure gradient and three-dimensional curvature of the streamlines. 

The variation of constants K and B and the wake strength pa- 
rameter is shown in Fig. 5.13. The value of constant 1/K increases in 
the streamwise direction up to streamline no. 3. However, on stream- 
lines no. 4 and 5, the value of 1/K decreases. 
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Fig. 5.8 Variation of Mean Velocity Profiles Along Streamline No. 
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Fig. 5. 1 1 Variation of Mean Velocity Profiles Along Streamline No. 
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Fig. 5. 13 Variation of 1/K, B, and n in Combined Law of the Wall and 
the Wake in Upstream Region of the Blade Leading Edge 
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The trend in the variation of B is just opposite that for 1/K. On 
the first three streamlines, the value of B decreases as the blade lead- 
ing edge is approached. However, on streamlines no. 4 and 5 the value 
of B increases in the streamwise direction. 

The variation of the wake strength parameter, n, along different 
streamlines is shown in Fig. 5.13(c). The variation in the value of II in 
the streamwise direction is found to be consistent with the observa- 
tions made in the logarithmic law plot (see Figs. 5.8(b) - 5.12(b)). 

Fig. 5.14 shows the variation of non-dimensional boundary layer 
integral parameters (8*/5, 8/0, H) along different streamlines. The 
value of 8*/8 increases as the blade leading edge is approached on all 
the streamlines, except on the streamline no. 5 (Fig. 5.14(a)). On 
streamline no. 5 the value of 8*/8 decreases in the streamwise direc- 
tion. 

The value of 8/0 decreases as the blade leading edge is ap- 
proached on all the streamlines, except on streamline no. 5 (Fig. 
5.14(b)). On streamline no. 5 the value of 8/0 increases slightly. 

The value of shape factor, H, increases in the streamwise direc- 
tion along different streamlines except the last two streamlines (Fig. 
5.14(c)). 

The variation of Clauser's shape parameter, G, and the skin fric- 
tion coefficient, Cf, is shown in Fig. 5.15. The value of G increases in 

the streamwise direction along all the streamlines except the last 
streamline (Fig. 5.15(a)). Along streamline no. 5, the value of G de- 
creases in the streamwise direction. 

Since there exists an adverse pressure gradient in the stream- 
wise direction and a favorable pressure gradient in the Y-direction, it 
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Fig. 5.14 Variation of Boundary Layer Integral Parameters Upstream 
of the Blade Leading Edge 
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Fig. 5.15 Variation of Clauser's Shape Parameter and Skin Friction 
Coefficient in Upstream Region of the Blade Leading Edge 
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is of interest to analyze the development of mean velocity profiles at 
different streamwise stations. 

Figs. 5.16 - 5.19 show the variation in the mean velocity profiles 
in the Y-direction. These profiles are presented in terms of the inner- 
and outer-layer variables at four different streamwise stations. At 
streamwise station 4 (x/D = - 2) the non-dimensional mean velocity 
profiles, in outer-layer variables, show lack of self-similarity in the Y- 
direction (Fig. 5.16(a)). This is possibly due to the presence of a high 
degree of skewness and strong favorable pressure gradient. Moreover, 
on the logarithmic law plot, a change in the values of K, B, and n are 
evident. The wake component of the boundary layer decreases as the 
distance y/D increases. This is also evident from the calculation of the 
wake strength parameter, n, (Fig. 5.13(c)). 

At x/D = - 4 (streamwise station 3) self-similarity in the mean 
velocity profiles, in outer-layer variables, at different y/D distances is 
observed (Fig. 5.17(a)). Also, at this streamwise station the wake com- 
ponent of the boundary layer decreases as y/D increases (Fig. 5.17(b)). 

At x/D = - 6 and - 8, a similaj* trend, as at streamwise stations 3 
and 4, in the variation of the wake component is evident(Figs. 5.18(b) 
and 5.19(b)). Furthermore, it is observed that the strength of the wake 
component at the last two streamwise stations, stations 1 and 2, is 
relatively smaller than that observed at stations 3 and 4. 

It is also interesting to note that various non-dimensional boun- 
dary layer integral parameters change in the Y-direction in an opposite 
manner to that observed for the streamwise direction (Fig. 5.14). For 
example, the value of 5*/5 decreases in the Y-direction, whereas it in- 
creases in the X-direction (Fig. 5.14(a)). It is quite evident from the 
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Fig. 5.17 Variation of Mean Velocity Profiles in Y-Direction at x/D = 
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Fig. 5.18 Variation of Mean Velocity Profiles in Y-Direction at x/D * 
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Fig. 5.19 Variation of Mean Velocity Profiles in Y-Direction at x/D = 


foregoing discussion that the observed differences in the variation of 
various boundary layer integral parameters in the X- and Y-directions 
can be attributed to the existence of varying adverse and favorable 
pressure gradients, respectively. 

Based on the measured mean velocity profiles, surface static 
pressure, and the local wall shear stress in an upstream region of the 
blade end-wall comer, it can be stated that the extent of influence of 
the blade leading edge in the upstream region is * 8D ( where D is the 
blade leading edge diameter). However, in the lateral direction (Y-di- 
rection) the flow is influenced by the blade leading edge up to a dis- 
tance of = 6D. 

5.2.5 Surface Static Pressure Variation in Blade End-Wall Comer 

Region 

Fig. 5.20(a) shows the variation of Cp on the flat plate surface 
(end-wall) in the blade end-wall comer region. At each streamwise 
station up to x/D =12, except at x/D = 2, the surface static pressure 
decreases as the comerline (the locus of point of intersection of the 
blade surface and the flat plate) is approached relative to the value far 
away from the cornerline. This variation occurs because the stream- 
lines are more curved near the comerline compared to those far away 
from the blade surface. However, at x/D = 2, the value of Cp first de- 
creases as the blade surface is approached and then increases by a 
small amount in the close vicinity of the comerline. Moreover, the 
trend in the Cp variation reverses downstream of x/D = 12. Also, the 
rate at which surface static pressure decreases at any streamwise sta- 
tion (except at x/D = 16 and 19.75) as the cornerline is approached is 
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larger upstream compared to that downstream of the maximum blade 
thickness location. In the streamwise direction, i.e., along the comer, 
the surface static pressure first decreases up to x/D = 10 and, there- 
after, it increases as the flow approaches the blade trailing edge. This 
trend in the variation of Cp along the comer is mainly because of the 
presence of the streamwise curvature. The flow is expected to ac- 
celerate from the blade leading edge up to its maximum thickness lo- 
cation and then decelerate as the flow approaches towards the blade 
trailing edge. 

On the blade surface, the static pressure variation at different 
streamwise stations is smaller in comparison to that observed on the 
end-wall surface (Fig. 5.20(b)). At each streamwise station up to x/D = 
12, the surface static pressure increases relative to values far away 
from the cornerline as the comerline (end-wall surface) is ap- 
proached. However, at the last two streamwise stations (x/D = 16 and 
19.75), the surface static pressure decreases as the end-wall surface is 
approached. In general, at any streamwise station, the observed in- 
crease or decrease in Cp towards the comerline is within 5% of the 
free stream dynamic pressure. Also, in the streamwise direction the 
surface static pressure variation is as expected on the blade surface 
i.e., it first decreases up to the blade maximum thickness location and, 
thereafter, it increases. Shabaka[36] reported that his preliminary 
measurements of surface static pressure on the wing surface of a sim- 
plified wing-body junction, at a given streamwise station, showed very 
little change along the span. Sepri[42] also measured the surface static 
pressure on the blade surface (airfoil section), at different streamwise 
stations, in a blade end-wall corner region and found very small 
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changes in Cp values along its span. The spanwise surface static pres- 
sure variation at different streamwise stations on the blade surface in a 
blade end-wall comer region observed in the present investigation is 
in good agreement with the results reported by other investigators 
[36. 42]. 

Fig. 5.21 shows the comparison of measured surface static pres- 
sure on the blade surface, at different streamwise stations and span- 
wise positions, with the values obtained from potential flow calcula- 
tions. The potential flow calculation is based on the vortex panel 
method described by Kuethe and Chow[144]. In this numerical tech- 
nique, the blade surface is replaced by vortex panels on which the cir- 
culation density varies linearly. A computer program developed also by 
Kuethe and Chow was used. This program required information about 
the blade geometry as input data. Therefore, the available computer 
program was modified to incorporate the generation of the blade ge- 
ometry. At each streamwise station the measured values of Cp at a dis- 
tance of z = 6D, i.e., far away from the end-wall surface, are generally 
higher (within 10% of the free stream dynamic pressure) compared to 
the values obtained from the potential flow calculation. Moreover, the 
differences in the values of Cp by the two methods are observed to be 
larger before the maximum thickness section of the blade than at 
streamwise stations following it. These observed differences can pos- 
sibly be because the potential flow calculation does not take into con- 
sideration the displacement thickness effect. However, the measured 
Cp distribution clearly shows that the trend in the variation of Cp is 
similar to that obtained by the potential flow calculation. 
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Fig. 5.2 1 Surface Static Pressure Distribution on Blade Surface- A Comparison 
Between Measured and Potential Flow Calculation (Case A) 



The measured distribution of Cp, at various streamwise stations 
on the blade surface, in the close vicinity of the end-wall surface (at 
z/D = 0.16) is also shown in Fig. 5.21. This distribution of Cp com- 
pared to its value far away from the end-wall surface (at z/D = 6) sug- 
gests that the surface static pressure on the blade surface, in blade 
end-wall corner region, is not influenced by the corner effects. The 
variation of the measured Cp on the end-wall surface upstream, in the 
comer region at a distance of y = 6D, and on the flat plate centerline 
(two-dimensional case) is also included in Fig. 5.21. The comparison of 
the Cp distribution for the above mentioned cases shows that Cp values 
in the three-dimensional case, upstream of the blade leading edge, are 
higher than the values for the two-dimensional flow. However, in the 
corner region, Cp values at different streamwise stations are consis- 
tently smaller than the values for the two-dimensional flow. The larger 
differences in Cp in the comer region compared to the values for two- 
dimensional flow implies the existence of corner effects up to dis- 
tances of y/D = 6. 

As mentioned in chapter III, the distribution of surface static 
pressure, at three streamwise stations, in the blade end-wall corner 
region was also obtained by using the conventional technique (i.e., wall 
taps). Fig. 5.22 shows the comparison of Cp obtained by the surface 
static tube with the values of Cp measured with the wall taps. The 
three streamwise locations were chosen such that the accuracy of the 
surface static tube could be examined in strong to mild favorable and 
adverse pressure gradients and also in the presence of varying skew- 
ness of the turbulent boundary layer. At all streamwise locations, ex- 
cept at x/D = 2 where the boundary layer is strongly skewed, Cp values 
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Fig. 5.22 Surface Static Pressure Variation In Blade End-Wall Comer 
Region- A Comparison Between Surface Static Tube and 
Conventional Methods (Case A) 


obtained from the static tube are In agreement with the values ob- 
tained from the wall taps within ± 2% (maximum) of the free stream 
dynamic pressure. However, at x/D = 2, the differences in Cp obtained 
with the two techniques are observed to be higher i.e., a maximum of 
7% of the free stream dynamic pressure. These large differences in 
Cp, at x/D = 2, resulting from the two different measurement tech- 
niques, probably can be attributed to the presence of large skewness 
and strong secondary flows. 

It is also observed that the absolute magnitude of the surface 
static pressure measured by the two techniques differ as much as 50% 
if expressed as a fraction of one another. Generally, if the measured 
surface static pressure is used for the calculation of velocity, the error 
in surface static pressure expressed as a fraction of free stream dy- 
namic pressure is commonly used. However, in calculating parameters 
such as the wall shear stress, the absolute magnitude of the surface 
static pressure is required. The calculations of the wall shear stress 
using the observed errors in surface static pressure were made. The 
results of the analysis showed that large errors (of the order of = 50%) 
in the surface static pressure measurements would cause very small 
errors in the calculated value of the wall shear stress. The maximum 
error in the wall shear stress was observed to be within 7%. 

5.2.6 Wall Shear Stress Variation in Blade End-Wall Comer Region 

The variation of normalized local wall shear stress on the end- 
wall surface at different streamwise stations in the comer region is 
shown in Fig. 5.23(a). At different streamwise locations up to x/D = 12, 
the wall shear stress increases as the blade surface is approached and 
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very close to the comerline its value decreases. However, at the last 
two streamwise locations, the value of wall shear stress first decreases 
as the cornerllne is approached and then increases. In the close 
vicinity of the blade surface it is found again to decrease. It is quite 
evident, that the maximum value of the normalized wall shear stress 

and its location from the comerline changes at different streamwise 
stations (Figs. 5.24 and 5.25). The maximum value of (x w /'tw e Jmax oc- 
curs at x/D » 6 and it is almost 2.4 times the value of wall shear stress 
far upstream (at x/D = - 20, y/D = 0) of the comer region as shown in 
Fig. 5.24. 

At a given streamwise location, on the blade surface, the local 
wall shear stress also attains a maximum value in a region closer to the 
end-wall surface and then decreases at far distances acquiring an 
asymptotic value (the wall shear stress values at z/D = 6) as shown in 
Fig. 5.23(b). Also, at streamwise stations downstream of the maximum 
blade thickness location, the trend in the variation of the wall shear 
stress in the Z-direction is different relative to the trend upstream of 
the maximum blade thickness section (Fig. 5.23(b)). Moreover, at each 
streamwise station, the maximum value of the wall shear stress and its 

location from the end-wall surface changes (Figs. 5.24 and 5.25). How- 
ever, the trend in the variation of the location of (%/x w Jmax on the 

blade surface differs from that observed on the end-wall surface (Fig. 
5.25). 

It is also interesting to note that, at different streamwise loca- 
tions in the comer region the maximum values of the normalized wall 
shear stress are consistently higher on the blade surface than the val- 
ues obtained on the end-wall surface of the test-model(Fig. 5.24). 
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Fig. 5.24 Variation of Maximum and Asymptotic Values of Wall Shear 
Stress in Blade End- Wall Comer Region for Case A 
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Fig. 5.25 Variation of the Location of Maximum Wall Shear Stress 
in Blade End- Wall Comer Region for Case A 
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These observed differences might possibly have been caused by the 
blade curvature. The asymptotic value (value of the wall shear stress at 
y/D = 6) of the wall shear stress at different streamwise locations on 
the end-wall surface in the comer region is approximately 35% higher 
than the value far upstream of the comer (Fig. 5.24). It is also impor- 
tant to mention that the local wall shear stress, at different stream- 
wise stations and at a distance of y/D = 6 is found to be higher 
(maximum 30%) than the value at corresponding locations in an ab- 
sence of the blade (two-dimensional case). Therefore, the comer ef- 
fects are dominant even up to distance of y = 6 times the blade leading 
edge diameter. 

On the blade surface preceding its maximum thickness section, 

especially at x/D = 2 and 6, the asymptotic values of the wall shear 
stress are of the same order as (Tw/'*w e Jmax- However, beyond x/D = 

6, asymptotic value of the wall shear stress decreases in the stream- 
wise direction (Fig. 5.24). 

The variation of the location of (Tw/tWoJmax a ^ on g the comer on 

both the surfaces of the test-model is shown in Fig. 5.24. On the end- 
wall surface, the location of the maximum value of the normalized wall 

shear stress moves from a distance of 0.3D to 0.9D from the blade 
surface. However, on the blade surface, the location of ('Cw/'tw 00 )max 

moves closer to the end-wall surface and then moves away from it at 
the last streamwise location, to a distance of = 3.6 times the blade 
leading edge diameter. 

The possible reasons for an increase in the wall shear stress val- 
ues in the vicinity of the comerline and changes in the maximum value 
and the location of the maximum value with respect to the cornerline 
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are discussed here. The plausible explanation is presented in combi- 
nation with the flow visualization and three-dimensional mean velocity 
profile measurements, on the same test- model, by Hazarika, Raj, and 
Boldman[77], A detailed discussion of the flow visualization and mean 
velocity profile measurements in the blade end-wall comer region is 
given in reference[47]. It is quite evident from the earlier investi- 
gations[36, 47, 77] in a blade end-wall comer region that a horseshoe 
shaped vortex is formed due to the deflection of a two-dimensional 
vortex sheet as it encounters the blade leading edge. The two legs of 
the horseshoe vortex wrap around the blade. As this horseshoe vortex 
moves along the comer it stretches or compresses depending on 
whether it encounters acceleration or deceleration in the flow. 
Moreover, at a given streamwise location in the comer region, the 
presence of the horseshoe vortex causes high momentum fluid to 
move from the outer region towards the cornerline on both the sur- 
faces. As a result, the mean velocity gradient is expected to increase 
near the surface consequently giving rise to an increase in the value of 
wall shear stress. Considering the end-wall surface of the test-model, 
the maximum value of the wall shear stress at any streamwise location 

increases up to the blade maximum thickness section (Fig. 5.24). This 
increase in (xw/xwjmax can be associated with the stretching of the 

horseshoe vortex because of acceleration of the flow in the streamwise 
direction. However, downstream of the blade maximum thickness lo- 
cation the horseshoe vortex gets attenuated because of deceleration of 
the flow and, as a result, the value of (Tw/tWoJmax decreases. The 

variation in the location of the maximum value of the wall shear stress 
along the comer on both surfaces of the test-model suggests that the 
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horseshoe vortex gets distorted as it moves in the downstream direc- 
tion. 

5.2.7 Wall Shear Stress Vectors in Blade End-Wall Comer Region 

The wall shear stress vectors obtained with the Preston tube at 
different streamwise locations on the end-wall surface are shown in 
Fig. 5.26(a). The dashed lines correspond to the limiting streamline 
directions obtained from the flow visualization study conducted by 
Hazarika[47] on the same test-model using an oil dot matrix method. 
The solid lines correspond to the wall shear stress directions obtained 
with the Preston tube. The magnitude of both solid and dashed line 
vectors is based on measurements with the Preston tube. A separation 
line, obtained from flow visualization, is also shown in Fig. 5.26(a). A 
variation in the wall shear stress direction, at a fixed streamwise sta- 
tion in the Y-direction, is clearly observed. The maximum value of the 
wall shear stress direction is observed at streamwise station 5 (x/D = 
2) located near the blade leading edge. At streamwise stations 5 and 6, 
the magnitude of the wall shear stress direction increases as the blade 
surface is approached. The flow visualization studies revealed that, up 
to the point of maximum thickness of the blade, the flow turned away 
from the comerline. Subsequently, in the region downstream, the flow 
was directed towards the cornerline. Also, the maximum magnitude 
and direction of the wall shear stress occurs between the blade surface 
and the separation line. 

The wall shear stress vectors at different streamwise stations on 
the blade surface are shown in Fig. 5.26(b). At x/D= 2, the magnitude 
of the wall shear stress direction is larger (maximum value 11 de- 
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Free End 


(b) Airfoil Surface 


Fig. 5.26 Wall Shear Stress Vectors in Blade End-Wall Comer Region 
for Case A (for Measurement Locations See Figs. 3.14, 3.16 
and 3.17) 
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grees) than the values at other streamwise locations. Also, a very small 
change in the wall shear stress direction In the Z- direction is ob- 
served. A similar trend in the variation of the limiting streamline di- 
rection is evident from the flow visualization study (see Tables 5.2 and 
5.3). 

The comparison of wall shear stress directions obtained by the 
two methods shows that the value of the angles obtained with the 
Preston tube generally differ by 6 to 8 degrees compared to those ob- 
tained from the flow visualization study. One of the reasons for this 
difference can be because the Preston tube measures the flow direc- 
tion at a distance equal to the effective center of the Preston tube 
opening. In a three-dimensional flow, the local flow direction changes 
in a direction normal to the surface. McAllister, Pierce, and Tennant 
[93] have also reported the existence of differences in the wall shear 
stress direction obtained from the velocity probe resting on the surf- 
ace and the limiting streamline directions obtained using the oil str- 
eak pattern. Another reason is the yaw sensitivity of the Preston tube. 
As discussed in section 5.2.1, the Preston tube used in the present 
study was found to be insensitive to the yaw angle of ± 5®. This is one 
reason why the wall shear stress directions obtained with the Preston 
tube beyond x/D = 6 on the end-wall surface were zero. The limiting 
streamline directions at streamwise stations beyond x/D = 6 show 
small angles within the range in which the Preston tube is insensitive. 
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Table 5.2 Measured Wall Shear Stress Directions With the Preston Tube on End-Wall 

Surface in Corner Region for Case A (Limiting Streamline Directions Obtained 
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SECTION B 


5.3 Three-Dimensional Plow (Case B) 

The measured mean velocity profiles at a stream wise location of 
x/D = - 20, the surface static pressures, and the wall shear stress vec- 
tors in the presence of the separation bubble in blade end-wall comer 
region are discussed in this section. 

5.3.1 Mean Velocity Profile Upstream of the Blade Leading Edge 

(Case B) 

The non-dimensional mean velocity profile in terms of outer- 
layer variables at a streamwise location of x/D = - 20 is shown in Fig. 
5.27(a). The non-dimensional mean velocity profile at the same 
streamwise location in an absence of the blade (two-dimensional flow 
case) is also included in Fig. 5.27(a). The mean velocity profiles for the 
two flow conditions (two- and three-dimensional flows) are clearly ob- 
served to be self-similar. 

Fig. 5.27(b) shows the mean velocity profiles in terms of inner 
law variables at x/D = - 20 for the two- and three-dimensional flows. 
The solid line in Fig. 5.27(b) corresponds to the logarithmic law of the 
wall with the universal constants mentioned earlier. For both the flow 
conditions the inner law of the wall exists but with different values of 
the universal constants in the logarithmic law of the wall. 

5.3.2 Surface Static Pressure Distribution in Blade End-Wall Comer 

Region (Case B) 

Fig. 5.28 shows the distribution of the surface static pressure at 
five streamwise locations on both surfaces in the blade end-wall comer 
region. On the end-wall surface at streamwise location up to x/D =12 
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the surface static pressure decreases as the blade surface is ap- 
proached (Fig. 5.28(a)). Also the rate in which the surface static pres- 
sure decreases changes at different streamwise locations up to x/D = 
12. At the streamwise location nearer to the blade leading edge (x/D = 
2) the values of Cp first decrease as the blade surface is approached. 
However, in the vicinity of the blade surface Cp increases. At a 
streamwise location of x/D = 16, the surface static pressure increases 
as the blade surface is approached. 

On the blade surface, the surface static pressure variation at a 
fixed streamwise station is smaller than on the end-wall surface (Fig. 
5.28(b)). A maximum change in Cp in the Z-direction of approximately 
7% of the free stream dynamic pressure is observed at a streamwise 
location of x/D = 2. 

The trend in the variation of the surface static pressure on both 
the surfaces of the test-model for Case B is similar to the trend for 
Case A. A detailed discussion of the comparison of different measured 
parameters for Case A and Case B is presented in section 5.4. 

5.3.3 Wall Shear Stress Variation in Blade End-Wall Comer Region 

(Case B) 

The variation of normalized wall shear stress at four streamwise 
stations on both surfaces of the test-model in the blade end-wall cor- 
ner region is shown in Fig. 5.29. On the end-wall surface at each 
streamwise location the magnitude of the wall shear stress increases 
as the blade surface is approached. However, in the vicinity of the 

blade surface the wall shear stress decreases (Fig. 5.29(a)). The maxi- 
mum value of tw/xwoo and its location with respect to the cornerline 
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(a) Plate Surface (Case B) (b) Blade Surface (Case B) 
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Fig. 5.29 Wall Shear Stress Variation in Blade End-Wall Comer Region for Case B 


(i.e., the blade surface) at each streamwlse station is different (Figs. 

5.30 and 5.31). Among the four streamwise stations considered, the 
maximum value of (Tw/%Jmax occurs at x/D = 2, and it is approxi- 
mately 2.2 times the far upstream value (at x/D = -20. y/D = 0) of wall 
shear stress as shown in Fig. 5.30. The asymptotic values of Xw/xw,*, at 

different streamwise locations are higher (by approximately 30%, ma- 
ximum) compared to the wall shear stress at x/D = - 20 (Fig. 5.30). 
These higher asymptotic values of Xw/xwoo. at y/D = 6. relative to far 

upstream value of wall shear stress, imply the existence of corner ef- 
fects on the wall shear stress up to distances of y = 6 times the blade 
leading edge diameter. 

On the blade surface at x/D = 2, the normalized value of wall 
shear stress increases rapidly in the vicinity of the end-wall surface 
(Fig. 5.29(b)). Furthermore, once the local wall shear stress attains its 
maximum value at x/D = 2, it does not decrease as is observed at other 
streamwise stations. A similar variation of the wall shear stress is also 
observed at streamwise location of x/D = 6 (see Fig. 5.37). It should be 
mentioned that at x/D = 2 and 6, a similar trend in the variation of lo- 
cal wall shear stress was observed for Case A (see Figs. 5.23(b) and 
5.37). There exists a distinctive maximum value of x w /xw« at each 

streamwise location. The location of the maximum value with respect 

to the end-wall surface varies in the downstream direction as shown in 
Fig. 5.31. Also, the maximum value of (tw/xwjmax occurs at x/D = 2 

where its value is 2.68, indicating an increase of 168% relative to the 
far upstream value of the wall shear stress (Fig. 5.30). 

At x/D = 12 on the blade surface, a significant decrease in the 
wall shear stress value is observed. The value of Xw/xw^ starts de- 
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Fig. 5.31 Variation of the Location of Maximum Wall Shear Stress 
for Case B 
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creasing sharply at z/D » 1.19 and continues to decrease up to z/D = 
3.8. In a distance of approximately 2.6D the value of TwAwoo decreased 

by * 80% (Fig. 5.29(b)). This large decrease in the wall shear stress 
value is directly related to the presence of the separation bubble. It 
can also be stated that in the Z-direction the separation bubble ex- 
tends from z/D * 1.2 to 2.6. Also, at x/D = 16 on the blade surface, the 
variation of the wall shear stress shows the usual trend i.e., in the 
vicinity of the end-wall surface, the wall shear stress increases. There- 
after, it decreases and achieves an asymptotic value. This variation at 
x/D =16 implies the flow has reattached to the blade surface again. 
The above mentioned discussion also suggests that the separation bub- 
ble formed somewhere after x/D = 10, and then the flow reattached 
before the streamwise location of the x/D = 16. These results are 
found to be in good agreement with the flow visualization studies (Fig. 
5.32). Fig. 5.32 shows a photographic view of the limiting streamlines 
obtained with the oil dot matrix technique. The separation bubble can 
be clearly seen in Fig. 5.32. In the Z-direction the separation bubble 
extends from a distance of z = 2D to 5.6D, whereas in the streamwise 
direction it starts after x = 10D and ends between x = 14D to 16D (Fig. 
5.32). 

The variation of the maximum and asymptotic values of TwAw*, 

on both surfaces of the test-model in the blade end-wall comer region 

is shown in Fig. 5.30. On both the surfaces the variation of 
(iwAwoJmax is similar to that observed in Case A (see Fig. 5.24). 

Therefore, the reasons given earlier are applicable for this case. A 
sudden drop in the asymptotic value observed at x/D = 12 is mainly 
attributed to the presence of the separation bubble. The maximum val- 
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Obtained From Flow Visualization for Case B 





ues of xw/xwoo at x/D =10 and 12 are smaller on the blade surface rel- 
ative to the values on the end-wall surface. This is also possibly caused 
by the presence of the separation bubble on the blade surface. 

Fig. 5.31 shows the variation of the location of the maximum val- 
ues of Xw/xwo,, on both the surfaces in the blade end-wall comer re- 
gion. On the end-wall surface, the location of (tw/xwjmax is observed 

to change from a distance of y = 0.32D to 0.53D. On the blade surface, 
the location of the maximum value of the wall shear stress first moves 
closer to the comerline (i.e., nearer to the end-wall surface) up to 
streamwise location of x/D =12 and, thereafter, it moves away from 
the end-wall surface as the flow progresses closer to the blade trailing 
edge. 

5.3.4 Wall Shear Stress Vectors in Blade End-Wall Comer Region 

(Case B) 

5.33(a) shows wall shear stress vectors at different streamwise 
locations on the end-wall surface of the test-model. For comparison 
purposes, the wall shear stress vectors for case A are also included. 
The wall shear stress vectors represented by solid and dashed lines 
correspond to Case A and Case B, respectively. Also, only at streamwise 
location of x/D = 2, the wall shear stress direction was obtained by 
using the Preston tube. However, at other streamwise locations, the 
directional information was obtained from the flow visualization study 
because the wall shear stress directions were found too small for the 
Preston tube to sense. On the end-wall surface, the wall shear stress at 
various measurement locations is higher for Case A than for Case B, 
whereas the magnitude of the wall shear stress directions for the two 
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Fig. 5.33(a) Comparison of Wall Shear Stress Vectors on the Flat Plate Surface 
of the Test-Model For Case A and Case B 


cases did not differ appreciably. One of the main reasons for the ob- 
served increase in the magnitude of the wall shear stress for Case A 
relative to Case B is the presence of an increased level of free stream 
turbulence in Case A. As mentioned in Chapter III, the turbulence 
generating grids were employed in Case A to avoid the formation of a 
separation bubble. Also, in Case B the free stream velocity was 23 m/s 
compared to 27.3 m/s for Case A. 

A comparison of the wall shear stress directions on the end-wall 
surface obtained from the flow visualization for Case A and Case B (see 
Tables 5.2 and 5.4) indicates small differences in the magnitude of the 
wall shear stress directions at different streamwise locations, except 
at x/D = 2. At a streamwise location of x/D = 2, the maximum differ- 
ence in singles for the two cases is approximately 5 to 6 degrees. 

On the blade surface the wall shear stress vectors at different 
streamwise locations are shown in Fig. 5.33(b). In this figure the vec- 
tors represented by dashed lines correspond to the limiting stream- 
line direction obtained from the flow visualization. The separation bub- 
ble is also traced by using the flow visualization picture. The wall shear 
stress directions obtained with the Preston tube, at different stream- 
wise locations, vary only slightly in the Z-direction (see also Table 5.5). 
The variation of the wall shear stress direction obtained from the flow 
visualization technique is similar to the variation obtained with the 
Preston tube. 
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Table 5.4 Measured Wall Shear Stress Directions With the Preston Tube on 
End-Wall Surface in Corner Region for Case B (Limiting Streamline 
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Table 5.5 Measured Wall Shear Stress Directions With the Preston Tube on Blade 
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5.4 Comparison of Surface Static Pressure and Wall Shear Stress for 

Case A and Case B 

As mentioned in Chapter III, the free stream velocity for Case A 
and Case B was different. Therefore, a quantitative comparison of the 
results can not be made. However, the results for the two cases are 
compared qualitatively. Furthermore, the effects of the presence of 
separation bubble on the measured surface static pressure and the wall 
shear stress are qualitatively analyzed. 

5.4.1 Comparison of Surface Static Pressure 

Fig. 5.34 shows a comparison of the surface static pressure at 
five streamwise locations on the end-wall surface of the test-model for 
Case A and Case B. A similar trend in the variation of the surface static 
pressure at different streamwise locations is observed for the two 
cases. No significant influence of the separation bubble is seen in the 
measured surface static pressure. 

The comparison of surface static pressure at five streamwise lo- 
cations on the blade surface for the two cases is shown in Fig. 5.35. At 
each of the streamwise locations, the variation in surface static pres- 
sure in the Z-direction is small (maximum of 7% of the free stream 
dynamic pressure). There is no influence of the separation bubble, 
which existed on the blade surface, on the surface static pressure dis- 
tribution. 

5.4.2 Comparison of Wall Shear Stress 

Fig. 5.36 shows the comparison of the normalized wall shear 
stress at four streamwise locations on the end-wall surface for case A 
and case B. The variation of normalized wall shear stress for the two 
cases at different streamwise locations clearly shows a similar trend. 
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Fig. 5.35 Surface Static Pressure Variation on the Blade Surface of the 
Test-Model- A Comparison Between Case A and Case B 
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Fig. 5.36 Wall Shear Stress Variation on the Flat Plate Surface of the 
Test-Model- A Comparison Between Case A and Case B 


The magnitude of Xw/xw,*, at eac h streamwise locations is in general 

higher for Case A than for Case B. The higher value of the wall shear 
stress for Case A is caused by the increased level of free stream turbu- 
lence intensity and increased magnitude of the free stream velocity. 

The comparison of normalized wall shear stress variation on the 
blade surface at five streamwise locations for the two cases is shown in 
Fig. 5.37. The variation of the wall shear stress at different streamwise 
stations (except at station 8 (x/D = 12)) is same for both the cases. At 
x/D =12, the variation of the wall shear stress is different for two 
cases. For Case B the value of normalized wall shear stress starts de- 
creasing rapidly at z = 1.2D and continues to decrease up to z = 3.8D. 
Thereafter, it changes only by a small amount (within 1%). This sharp 
decrease in the value of normalized wall shear stress is mainly due to 
the presence of the separation bubble in Case B. 

Fig. 5.38(a) shows a comparison of the maximum and the 
asymptotic values of the normalized wall shear stress on the end-wall 

surface at different streamwise stations for the two cases. A similar 
trend in the variation of maximum and asymptotic values of xw/xw*, for 

the two cases is observed. 

A comparison of the maximum and the asymptotic values of 
xw/xwoo on the blade surface for the two cases is shown in Fig. 5.38(b). 

A similar variation of the maximum and the asymptotic values of the 

normalized wall shear stress for the two cases is evident in the figure. 
Also, both the maximum and the asymptotic values of xw/xwoo f° r Case 

A are higher relative to the values for case B. The decrease in the 
maximum and the asymptotic values of Xw/xw^ for Case B can be at- 
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(a) Blade Surface (Case A) (b) Blade Surface (Case B) 
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Fig. 5.37 Wall Shear Stress Variation on the Blade Surface of the 
Test-Model* A Comparison Between Case A and Case B 




Fig. 5.38 Variation of Maximum and Asymptotic Values of Wall Shear 
Stress- A Comparison Between Case A and Case B 
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tributed to a lower level of free stream turbulence intensity and lower 
free stream velocity in case B. 

Fig. 5.39(a) shows a comparison of the location of the maximum 
value of tw/Twoo on the end-wall surface for the two cases. A similar 

variation of the location of (iw/ %Jmax for two cases is evident in Fig. 

5.39(a). Also, the location of the maximum value of the normalized wall 
shear stress at different streamwise stations is almost identical. 

A comparison of the location of the maximum values of the nor- 
malized wall shear stress on the blade surface for the two cases is 

shown in Fig. 5.39(b). A similar variation of the location of the maxi- 
mum value of iw/xwoo is noticed for the two cases. However, in Case B 

the location of (xw/xwjmax moves closer to the end-wall surface. 
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CHAPTER VI 


RESULTS AND DISCUSSION ON THE WALL PRESSURE 
FLUCTUATION AND ITS CORRELATION WITH THE VELOCITY 
FLUCTUATION IN BLADE END-WALL CORNER REGION 


Experimental results of the wall pressure fluctuation and its cor- 
relation with the velocity fluctuation in the blade end-wall comer re- 
gion are presented in this chapter. Measurements of the wall pressure 
fluctuation and its correlation with the velocity fluctuation in the ab- 
sence of the blade i.e., for two-dimensional flow, were also obtained 
and compared with the blade end-wall comer flow results. 

6.1 Two-Dimensional Flow 

6.1.1 Wall Pressure Fluctuation Intensity 

The RMS value or the intensity of the wall pressure fluctuation, 
p’ w , normalized by the free stream dynamic pressure, qoo obtained in 
the present investigation is 2.79 xl0‘2. The value of p'w/floo obtained 
in the present study is higher by a factor of approximately 3 compared 
to the values reported in the literature (see Table 2.4). The non-di- 
mensional diameter, di + (= diu*/v), of the pressure transducer has 

been reported to influence the value of the wall pressure fluctuation 
intensity[52,1271. The value of p’w/floo increased with a decrease in 
the value of di + mainly because of the improved spatial resolution[52J. 
The value of dj + = 60, in the present study, lies in the lower range of 

the values reported by other investigators (Table 2.4). The other pos- 
sible reason for the observed high value of p’w/ Too in the present study 

appears to be the presence of the initial roughness or the combined 
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effects of the initial roughness and the increased level of free stream 
turbulence. The results presented in Chapter IV indicated that the 
mean flow boundary layer parameters and the inner region of the 
boundary layer were influenced by the presence of the initial rough- 
ness. Therefore, the wall pressure fluctuation intensity is expected to 
be influenced by the presence of the initial roughness. 

The value of p’ w was also measured without the turbulence gen- 
erating grids at the same streamwise location to estimate the effect of 
increased free stream turbulence level on the value of p' w . The value of 
p’ w increased by only approximately 1.4% with the increased level of 

free stream turbulence. This result suggests that the high value of 
P'w/ floo obtained in the present investigation is probably caused by the 

initial roughness. As a consequence of high wall pressure fluctuation 
intensity, the value of p'w/ x w is also relatively larger than the usual 
value of * 3[52]. In the present study, the value of p' w / x w obtained is 
8.4. 

6.1.2 Frequency Spectra of the Wall Pressure Fluctuation 

The frequency spectral distribution with and without the tur- 
bulence generating grids is shown in Fig. 6.1. No significant change is 
observed for the two cases at frequencies below 3500 Hz. However, 
the spectral content increases at frequencies above 3500 Hz. This im- 
plies that the increase in free stream turbulence level increases the 
spectral content of the wall pressure fluctuation. This is in agreement 
with the observed increase in the overall wall pressure fluctuation in- 
tensity because of the increased free stream turbulence level. 
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6.1.3 Wall Pressure-Velocity Correlation 

The correlation of the wall pressure fluctuation with the 
streamwise fluctuating velocity component across the boundary layer is 
shown in Fig. 6.2. The results for the case in which turbulence gener- 
ating grids were not used are also included. The trend in the variation 
of Rpu for z/8 > 0.15 is observed to be the same for both the cases. For 
0.15 < z/8 < 0.6, the value of Rp U first increases up to a distance of z/8 
~ 0.4 and, thereafter, slightly decreases. The maximum value of Rpu is 
approximately 0.05 for both of the cases. The turbulent boundary layer 
can be divided into two regions, namely, the inner and the outer re- 
gions. In the inner region, 0.15 < z/8 < 0.4, the value of Rpu increases 
and attains its maximum value. This implies that the fluctuating veloc- 
ity field in this region contributes more to the wall pressure fluctua- 
tion field than the outer region. In the outer region, z/8 > 0.4, the 
value of Rpu decreases compared to its value in the inner region which 
is the indication of a smaller contribution from the outer region. Also, 
close to the boundary layer edge the value of Rpu increases again. The 
increase in the value of Rp U in the inner region was also reported by 
Serafini[53J. Moreover, Serafini's study revealed that the intermittent 
region of the turbulent boundary layer contributed approximately half 
to the wall pressure fluctuation field than that of the inner region. 
However, in the present study, the values of Rpu close to the boundary 
layer edge are found to be of the same order as the maximum value of 
Rpu in the inner region. 
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6.2 Three-Dimensional Flow- Upstream of the Blade Leading Edge 

6.2.1 Wall Pressure Fluctuation Intensity 

The intensity of the wall pressure fluctuation measured up- 
stream of the blade leading edge at x/D = - 8 is the same as obtained 
for the current two-dimensional case. Therefore, from the point of 
view of the wall pressure fluctuation field, the flow in the upstream 
region is not influenced by the presence of the blade at an axial dis- 
tance of 8D. 

6.2.2 Frequency Spectra of the Wall Pressure Fluctuation 

Fig. 6.3 shows the frequency spectral distribution at x/D = - 8. 
For comparison, the frequency spectral distribution at the same strea- 
mwise location in absence of the blade is also included in Fig. 6.3. It is 
evident that at streamwise location, x/D = - 8, the spectral content at 
different frequencies with and without the blade is the same. This 
again implies, as mentioned earlier, that in the upstream region the 
wall pressure fluctuation field is not influenced by the presence of the 
blade at x/D = - 8. 

6.2.3 Wall Pressure-Velocity Correlation 

The variation of correlation coefficient, Rpu in the Z-direction 
upstream of the blade leading edge at x/D = - 8 is shown in Fig. 6.4. At 
the same streamwise location, the results obtained in the absence of 
the blade i.e., for the two-dimensional case are also included in Fig. 
6.4. A similar trend in the variation of Rpu for the two cases is evident. 
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6.3 Three-Dimensional Flow- Blade End-Wall Comer Region 
6.3.1 Wall Pressure Fluctuation Intensity 

Fig. 6.5(a) shows variation of the wall pressure fluctuation inten- 
sity at three streamwise stations in the blade end-wall comer region 
normalized by the far upstream value (p’wO) measured at a streamwise 

location of x/D = - 8. At all of the streamwise stations in the blade end- 
wall comer region the value of p' w is less than the far upstream value. 

At x/D = 2, the normalized wall pressure fluctuation intensity is ob- 
served to increase in a small distance close to the blade surface and, 
thereafter, it decreases with some scatter. The maximum value of the 
normalized wall pressure fluctuation intensity is found to be 0.75. Also, 
the maximum value of the normalized wall pressure fluctuation inten- 
sity is located close to the separation line. The traces of the separation 
line obtained from the flow visualization study of Hazarika[47] are 

shown in Fig. 5.26. The variation of normalized wall pressure fluctua- 
tion intensity in the Y-direction at x/D = 2 also shows that p’ w /P’wO 

increases with a decrease in the mean surface pressure, whereas, its 
value decreases as the mean surface pressure increases (see Fig. 
5.20(a) for mean surface static pressure distributions). 

At x/D = 10, the normalized wall pressure fluctuation intensity 

has a maximum value of 0.83 which occurs very close to the blade sur- 
face as shown in Fig. 6.5(a). The values of p' w /P’wO at x/D =10 are 

consistently higher than the values at x/D = 2. In the streamwise di- 
rection the mean surface static pressure, as discussed in Chapter V 
(Fig. 5.20(a)), decreased from x/D = 2 to x/D =10 and, thereafter, its 
value increased. In the Y-direction at x/D = 10, the mean surface static 
pressure increases, whereas, the wall pressure fluctuation intensity 
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Fig. 6.5 Variation of Normalized Wall Pressure Fluctuation Intensity 

in Blade End-Wall Comer Region: (a) Normalized by the Value 
Far Upstream of the Blade Leading Edge; (b) Normalized by 
the Local Wall Shear Stress 
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decreases with some scatter (see Fig. 5.20(a) and 6.5(a)). In the 
streamwise direction, as in the Y-direction, the wall pressure fluctua- 
tion intensity increases with the decrease in surface static pressure. 

At x/D =16, very close to the blade surface the normalized wall 
pressure fluctuation intensity has a maximum value of 0.68 (Fig. 
6.5(a)). In the Y-direction the normalized wall pressure fluctuation in- 
tensity first decreases and attains a second maximum value of ap- 
proximately 0.67 and, thereafter, it slowly decreases with some scat- 
ter. The mean surface static pressure in the Y-direction at x/D =16 
slightly decreases up to y/D = 1.5, whereas, the value of wall pressure 
fluctuation intensity increases, except very close to the blade surface 
(see Figs. 5.20(a) and 6.5(a)). At distances of y/D > 1.5, the mean sur- 
face static pressure and the wall pressure fluctuation intensity values 
show small changes. In the streamwise direction the adverse pressure 
gradient exists beyond x/D = 10 (see Fig. 5.20(a)). However, the values 
of wall pressure fluctuation intensity at x/D =16 are smaller than the 
values at x/D = 2 and x/D = 10 as shown in Fig. 6.5(a). 

The above discussion on the variation of the wall pressure 
fluctuation intensity at different streamwise stations in the blade end- 
wall comer region clearly indicates that the wall pressure fluctuation 
intensity increases in the presence of a favorable pressure gradient, 
whereas, its value decreases in the presence of adverse pressure gra- 
dient. The aforementioned observation has been observed both in the 
streamwise and lateral directions. However, the available measure- 
ments in a two-dimensional turbulent boundary layer showed a higher 
value of wall pressure fluctuation intensity in the presence of adverse 
pressure gradient and lower value in the presence of favorable pres- 
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sure gradient compared to the zero pressure gradient case[131). At all 
the three streamwise stations examined in the blade end-wall comer 
region, the maximum normalized wall pressure fluctuation intensity 
occurred between the blade surface and the separation line. 

The presence of the horseshoe vortex in the blade end-wall cor- 
ner region and the changes in its strength appear to significantly in- 
fluence the variation of the wall pressure fluctuation intensity. The ob- 
served increase in values of wall pressure fluctuation intensity at x/D = 
10 compared to the values at x/D = 2 is probably associated with the 
stretching of the horseshoe vortex due to the existence of a favorable 
pressure gradient in this region. Moreover, beyond x/D = 10 the 
horseshoe vortex becomes attenuated because of the adverse pressure 
gradient. Therefore, the observed decrease in the values of the wall 
pressure fluctuation intensity at x/D =16 can possibly be attributed to 
the attenuation of the horseshoe vortex. 

The measurements of wall pressure fluctuation intensity, in the 
simplified wing-body junction, were recently also reported by Hasan, 
Casarella and Rood[130J. The wing used in their study had a semi-el- 
liptic leading edge followed by a constant thickness section for some 
length and then a tapered tail section. The investigation of Hasan et al. 
revealed that the normalized values (normalized by the value far up- 
stream of the wing leading edge) of wall pressure fluctuation intensity 
in the blade end-wall comer region were greater than unity. However, 
the normalized values of wall pressure fluctuation intensity are found 
to be less than unity in the present investigation. The value of p' w 0> 

used for normalizing the wall pressure fluctuation intensity values in 
the corner region in the present investigation, is relatively high corn- 
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pared to the values for two-dimensional flows available in the litera- 
ture. The high value of p'wo obtained in the present study might have 

caused the observed difference with the study of Hasan et al. 

Fig. 6.5(b) shows the variation in the Y-direction of the wall 
pressure fluctuation intensity normalized by the local time averaged 
wall shear stress. The value of p'w/^w varies from approximately 3-6. 

6.3.2 Frequency Spectra of the Wall Pressure Fluctuation 

Fig. 6.6(a) shows the variation of the frequency spectral distribu- 
tion in the Y-direction at x/D = 2. The spectral content at different 
frequencies increases, except at y/D = 1.65, as the blade surface is ap- 
proached. However, at y/D = 1.65 the spectral content for f < 2500 Hz 
is higher and for f > 2500 Hz is smaller compared to the value at y/D = 
1.02. The observed increase in spectral content at y/D = 1.65 is con- 
sistent with the maximum value of wall pressure fluctuation intensity 
obtained at this streamwise station (see Fig. 6.5(a)). 

The variation of frequency spectral distribution in the Y-direc- 
tion at x/D = 10 is shown in Fig. 6.6(b). Although the change in spectral 
content for different y locations is small at this streamwise station the 
spectral content is observed to increase as the blade surface is ap- 
proached. 

The variation of frequency spectral distribution in the Y-direc- 
tion at x/D = 16 is shown in Fig.6.6(c). At this streamwise station the 
variation of the spectral distribution with different y-locations is simi- 
lar to that observed at x/D = 2. However, the changes in spectral con- 
tent at the different y locations are not significant. 
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6.3.3 Wall Pressure -Velocity Correlation 

The variation of correlation coefficient normal to the end-wall at 
three streamwise stations in the blade end-wall comer region is 
shown in Fig.6.7. Firstly, it is observed that a large scatter in the values 
of correlation coefficient exists. Also, at each streamwise station the 
variation of Rpu at the two y-locations is similar within the 
experimental scatter. In general, it is noticed that for z/8 < 0.6 the 
magnitude of correlation coefficient increases with the increase in 
distance in the Y-direction. However, this trend appears to reverse 
near the shear layer edge. 

The values of Rpu beyond z/8 > 0.2 at x/D =10 are smaller than 
the values at x/D = 2. These lower values of Rp U are due to the exis- 
tence of a favorable pressure gradient between these two streamwise 
stations. However, at x/D = 16, the values of Rp U beyond z/8 > 0.2 are 
higher than the values at the other two streamwise stations with some 
scatter. The existence of a strong adverse pressure gradient(Fig. 
5.20(a)) beyond the streamwise location of x/D =10 appears to be the 
reason for the observed increase in the values of Rp U at x/D = 16. In 
the two-dimensional turbulent boundary layer flow, the increase in the 
value of Rpu in the presence of the adverse pressure gradient relative 
to the zero pressure gradient case was reported by Bradshaw[49]. 
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Fig. 6.7 Variation of Wall Pressure-Velocity Correlation Coefficient 
at Different Streamwise Locations in Blade End-Wall 
Comer Region 





CHAPTER VH 


SUMMARY OF RESULTS 

The combined effects of initial roughness and the increased 
level of free stream turbulence on the development of a two-dimen- 
sional turbulent boundary layer were investigated. The wall shear 
stress and the mean surface static pressure measurements, upstream 
of the blade leading edge and in the blade end- wall corner region, 
were obtained with and without a separation bubble on the blade. Mea- 
surements of wall pressure fluctuation and its correlation with the ve- 
locity fluctuation were made upstream of the blade and in the blade 
end-wall comer region. The results of the above mentioned measure- 
ments are presented and discussed in Chapters 4-6. Based on the re- 
sults of the various measurements, the following main conclusions 
were drawn: 

(1) The surface static pressure on the end-wall varied significantly 
in the streamwise and lateral (Y) directions upstream of the 
blade leading edge (Fig. 5.3). An adverse pressure gradient ex- 
isted in the streamwise direction because the flow decelerated 
as it approached the blade leading edge. However, the surface 
static pressure, at different streamwise stations, decreased in 
the Y-direction. The rate of decrease of the surface static pres- 
sure in Y-direction was observed to be a maximum at x = - 2D 
where the streamlines had large curvature due to the presence 
of the blade leading edge. The values of Cp, at different stream- 
wise stations, were observed to approach the far upstream value 
(value at x/D = - 20) within 2% of the free stream dynamic pres- 
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sure at a distance of y = 6D. This result implied the influence of 
the blade on the surface static pressure diminished within a dis- 
tance of approximately 6D in the Y-direction. 

(2) The wall shear stress decreased in the upstream region as the 
blade leading edge is approached. The maximum decrease in the 
value of wall shear stress was approximately 40% relative to the 
value far upstream of the blade leading edge. This decrease in 
wall shear stress occurred on the stagnation streamline where 
the maximum adverse pressure gradient was observed (Fig. 5.5). 
The wall shear stress and surface static pressure in the up- 
stream region at different streamwise stations approached the 
fax upstream values within a distance of approximately 6D in the 
Y-direction. 

(3) The mean velocity profiles obtained along different streamlines 
in the upstream region of the blade leading edge indicated vari- 
ation in the values of the constants ( K, B and n) in Eqn. 5.2 (Fig. 
5.13). The values of 1/K (the coefficient of the term In z+ in 
Eqn. 5.2) and n increased, whereas, the value of B decreased in 
the streamwise direction. This trend prevailed along the stream- 
lines nearest to the stagnation streamline. However, the trend in 
the variation of these constants reversed along the streamlines 
far away from the stagnation streamline. The observed variation 
in the values of these constants was attributed to the changes in 
adverse pressure gradient and curvature of the streamlines. The 
mean flow boundary layer parameters in the upstream region 
were influenced by the changes in pressure gradient (Figs. 5.14 
and 5.15). The various boundary layer parameters showed an 
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opposite trend in the streamwise and lateral directions because 
of the presence of adverse and favorable pressure gradients, 
respectively. 

(4) At different streamwise stations, in the blade end -wall corner 
region, the surface static pressure varied significantly on the 
end- wall surface. However, on the blade surface at any stream- 
wise station a small variation (< 5% of the free stream dynamic 
pressure) in the surface static pressure existed along its span 
(Fig. 5.20). 

(5) At different streamwise stations on the end- wall surface, except 
at first station in the comer region, the surface static pressure 
measurements obtained by the surface static tube and the wall 
taps were in good agreement (within ± 2% of the free stream 
dynamic pressure (Fig. 5.22)). However, at the first streamwise 
station (x/D = 2), the differences in Cp obtained by the two 
techniques were higher (maximum 7% of the free stream dy- 
namic pressure). These observed large differences in Cp at x/D = 
2 resulting from the two methods of measurements were at- 
tributed to the presence of large skewness in the flow. 

(6) At different streamwise stations on the end-wall surface, except 
at the last two stations, the normalized wall shear stress values 
increased in the corner region as the blade surface was ap- 
proached (Fig. 5.23(a)). The maximum value of the normalized 
wall shear stress and also its location from the blade surface 
were observed to change at different streamwise stations (Figs. 
5.24 and 5.25). 
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(7) The variation of the wall shear stress on the blade surface was 
not the same as on the end-wall surface. However, on both the 
blade surface and the end-wall surface, the maximum value of 
the normalized wall shear stress and its location from the end- 
wall surface changed at different streamwise stations (Fig. 5.24 
and 5.25). Moreover, the maximum values of the normalized wall 
shear stress at different streamwise stations in the corner re- 
gion were found to be consistently higher on the blade surface 
compared to the values on the end-wall surface (Fig. 5.24). 

(8) The trend in the variation of the wall shear stress directions ob- 
tained from the Preston tube measurements, on both the sur- 
faces of the test-model, was in agreement with the values ob- 
tained from the flow visualization method. The differences in the 
wall shear stress directions obtained by the two methods were 
within the ranges of yaw angles in which the Preston tube was 
found to be insensitive. 

(9) In the presence of a separation bubble (Case B), the wall shear 
stress variation at different streamwise stations on the end-wall 
surface was similar to that observed without the separation bub- 
ble, Case A (Fig. 5.36). On the blade surface, a similar trend in 
the variation of wall shear stress at different streamwise stations, 
except at x/D = 12, for the two cases was observed (Fig. 5.37). At 
x/D = 12 on the blade surface, for Case B, the wall shear stress 
decreased sharply at z/D ■ 1.2 and continued to decrease up to 
z/D * 3.8. This decrease in the wall shear stress was mainly 
caused by the presence of the separation bubble . 
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( 1 0) The values of wall pressure fluctuation intensity on the end-wall 
surface, at different streamwise stations in the blade end-wall 
comer region, were significantly influenced by the changes in 
the strength of the horseshoe vortex as it moved along the cor- 
ner (Fig. 6.5(a)). At any streamwise station in the blade end-wall 
corner region, the spectral contents at different frequencies 
were observed to increase as the blade surface was approached 
(Fig. 6.6). 

(11) The values of correlation coefficient, in the streamwise direc- 
tion, were observed to decrease and increase because of the 
presence of favorable and adverse pressure gradients, respec- 
tively, in the blade end-wall comer region. Also, at distances be- 
yond z/5 = 0.2, the values of R pu were higher in the inner region 

compared to the outer region of the shear layer. 

(12) The development of turbulent boundary layer in the presence of 
natural transition (Case 1) showed the existence of self-similarity 
at all streamwise stations, except at the first streamwise station 
(Fig. 4.1(a)). The lack of self-similarity at the first streamwise 
station suggested that the boundary layer had not attained its 
fully developed state. Similar results on the existence of self- 
similarity at different streamwise stations were evident either if 
initial roughness was introduced. Case 3 (Fig. 4.13(a)) or when 
the free stream turbulence level was increased with or without 
the presence of initial roughness, viz.. Cases 4 and 2, respec- 
tively (see Figs. 4.15(a) and 4.7(a)). 

(13) The velocity-defect was observed to decrease at a given stream- 
wise station with an increase in the free stream turbulence in 
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the absence of initial roughness (Fig. 4.9). A fullness of the mean 
velocity profile indicates an increase in the value of wall shear 
stress. The measured wall shear stress values in Case 2 were 
found to be higher than the values in Case 1 and, therefore, sup- 
port the above observation (Fig. 4.12). However, the velocity-de- 
fect increased and the mean velocity profiles became less full 
with an increased level of free stream turbulence in the pres- 
ence of initial roughness (see Figs. 4.8 and 4.9). 

( 1 4) The comparison of Cf values at the same streamwise stations 

with and without initial roughness i.e., Case 3 with Case 1, indi- 
cated lower values of Cf closer to the roughness strip. Further 
downstream Cf increased relative to Case 1. It was noted that a 
distance of « 40 8o (8o being the boundary layer thickness at the 

first streamwise station in Case 1) from the roughness strip was 
required for Cf to reach the values obtained in Case 1 (Fig. 
4.12(a)). However, the comparison of Cf results for Case 4 with 
Case 1 indicated that a larger distance « 65 8o downstream of 
the roughness strip was needed to reach the Cf values of Case 1 . 

(15) A comparison of mean flow boundary layer parameters for four 
cases revealed that in the presence of initial roughness, the in- 
creased level of free stream turbulence had an effect opposite in 
nature compared to the case when initial roughness was absent 
(see Tables 4.2 and 4.3 and Figs. 4.3-4. 6). The variation of wall 
shear stress appeared to be one of the possible reasons for this 
observed trend. 
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APPENDIX A 


NORMAL SENSOR HOT-WIRE PROBE CALIBRATION 

The normal sensor hot-wire probe was calibrated against the 
pitot- static tube. The anemometer output E 5 and the velocity U are 

related by King’s law written as, 

Eb 2 = Ai+BiU n (A. 1) 

In the above equation Ai, Bi and n are calibration constants which are 
obtained by calibrating the hot-wire probe. The hot-wire and the pitot- 
static probes were placed in the wind tunnel test-section where the 
uniform flow existed. Corresponding to the known velocity, which was 
measured with the pitot-static probe, anemometer output was recor- 
ded. The outputs from both the probes were simultaneously digitized 
using the MASSCOMP mini computer system. The number of data sets, 
containing hot-wire and pitot-static probes outputs, were obtained by 

varying the free stream velocity. For different values of exponent n, 
constants Aj and Bj in Eqn. A.1 were calculated in a least square 

sense. The optimum values of Aj, Bi and n were obtained such that the 

normalized standard deviation for calculated velocity is minimized. 
The normalized standard deviation for velocity was calculated using 
the relation, 

e u = [ (1/N1) 1(1- U c /U)2 1 0.5 (A. 2) 

where U is the measured velocity, U c is the value of velocity calculated 

by using Eqn. A.1 and N1 is the number of measured data sets. The 
above mentioned procedure for calculating the optimum values of Aj, 
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Bi and n was proposed by Bruun and Tropea[145]. A typical calibration 
curve is shown in Fig. A.l. Fig. A.2 shows the variation of normalized 
standard deviation for the calculated velocity against the value of expo- 
nent n. It is quite evident that the optimum value of n exists at which 
the error in calculated velocity is minimum (Fig. A.2). A listing of the 
computer program used for calibrating the hot-wire probe using the 
digitized data is also included here. 


Eb 2 (Volt) 





nnnnnnnnnnunnnnnnnnnrin 


INCLUDE '/USR/INCLUDE/MR.F • 


PROGRAM TO CALCULATE OPTIMUM VALUES OF CALIBRATION 
CONSTANTS IN KING'S LAW IN LEAST SQUARE SENSE FOR 
SINGLE NORMAL SENSOR MOT-WIRE PROBE 
PROGRAM REAOS DIGITIZED DATA STORED ON OISC FILES 
UNDER NAME 'DAT .1 JK ' 

PROGRAM WRITTEN IN USER INTERACTIVE MODE 
CALOATA - DIGITIZED SAMPLES 
NCH - TOTAL NO. CF DIGITIZED CHANNELS 
TNSAMPL - TOTAL NO. OF DIGITIZED SAMPLES COLLECTED 
RFRQ - SAMPLING RATE 

BRFRQ - THROUGHPUT RATE OF A/D UNIT * 1MHZ 
M - NO. OF DATA SETS 

CC * CALIBRATION CONST. FOR VARIABLE RELUCTANCE TYPE 
DIFFERENTIAL PRESSURE TRANSDUCER IN IN. /VOLT 
CF1 * CONVERSION FACTOR FROM 1IN. OF HG TO N/M**2 
CF2 * CONVERSION FACTOR FROM 1IN. OF WATER TO N/M**2 
RAIR * GAS CONST. FOR AIR IN J/KG.K 
AVMNl(I) AND AVMN2U) ARE USED AS OUTPUTS OF CH. 1 
AND CH.2, RESPECTIVELY 

AVMN1 IS USED TO CALCULATE MEASURED VELOCITY. UM(I) 

INTEGERS CALDATAOOOGOC) 

DIMENSION AVMN1 (30) ,AVMN2( 30 ) , UM (30) 

INTEGER NCH, TNSAMPL, M.NTCH 
REAL RFRO, BRFRQ 
CHARACTER*20 FLNM 
CHARACTER*20 OUTFILE 

PRINT *, 'ENTER FILE NAME WHERE DATA ARE STORED' 

READ 1001, FLNM 
1001 FCRMAT(A) 

QPEN( 9, FORM* 'UNFORMATTED '.STATUS* 'OLD' ,FILE*FLNM) 
REWIND 9 

READ(9) NCH, TNSAMPL, RFRO. BRFRQ. M 
PRINT +, NCH, TNSAMPL, RFRQ, BRFRO.M 
NTCH * TNSAMPL 
NSAMPL * TNSAMPL/NCH 

CALL DISCREADKALDATA, NTCH, NCH, M, RFRO, AVMN1.AVMN2) 
PRINT +, 'ENTER FILE NAME FOR OUTPUT AS CALOUT.IJ' 
READ 1001, OUTFILE 
PRINT 1001, OUTFILE 

OPEN (8, FORM* 'FORMATTED' »FILE*OUTFILE ) 

REWIND 8 

WRITE (8 ,101) TEXT 
101 FORMAT ( 2X ,A) 

DC 10 I * 1, M 

10 PRINT *, 'DATASET £• , I , AVMN1 ( I ) , AVMN2 ( I ) 

PRINT ♦, ' ENTER TRANSDUCER SENSITIVITY IN IN. /VOLT' 
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READ ♦ , CC 

PRINT +, • ENTER ATM. PRESSURE IN IN. OF HG = PATH' 

READ ♦ , PATM 

PRINT *» * ENTER ROOM TEMP. IN DEG. F = TATM' 

READ +. TATM 
CF1 = 3386.6 
CF2 = 249.09 
RAIR * 287.1 

TABS * 5.0*(TATM - 32.0J/9.0 ♦ 273.15 
RHO = PATM * CF 1/ (RA IR 4 TABS) 

WRITE C 8, 102) RFRQ, NSAMPl, PATM, TABS, CC ,RHO 

102 FCRMAT(/2X, 'SAMPLING RATE *• ,F7 . 1 , 1 OX SAMPLE S/CHANNEL = ', 

1 I6//2X » ' ATM . PRESS. (IN. HG) * ' ,F7 .3 »5X • 'ATM. TEMP (K ) = ', 

2 F8.3//2X, 'TRANSDUCER SENSITIVITY =',F7.4//2X, 

3 'AIR DENSITY =',F7.4/I 
WRITE (3,103) 

103 FORMAT (/5X» ' I ' »2X » 'UREF (VOL TS ) ' »2X»'BRIDGE OUTPUT (VOLT S )'/ ) 
DO 11 I = 1 » M 

11 WRI TE (8 » 104 ) I .AVMNKI ) .AVMN2(I) 

104 FORMAT (3X,I2,2(4X,F8.5)/) 

DO 20 I = 1 »M 

UM ( I ) * SORT (2 .0+AVMN1 ( I)*CC*CF2/RH01 
20 PRINT *, UM ( I ) « AVMN2 ( I ) 

CALL CALBRT (UM , AVMN2 ,M , ACP T,BOPT»OPTN,ERMIN) 

WRITE (8,107) AOPT.BCPT.GPTN 

107 F0RMAT(//5X, 'OPTIMUM A = • . F8 .5//5X , • OPT I MUM B =',F8.5// 

1 5X, 'OPTIMUM N *• »F8.5/J 

CLOSE (8) 

STOP 

END 


SUBROUTINE DISCREAD 

READS DIGITIZED DATA AND CALCULATES MEAN VALUES OF 
DATA (IN YOLTS) FROM TWO CHANNELS 

SUBROUTINE D I SCREAD ( CALDAT A »NTCH ,NCH ,M , RFRQ , VMN1 ,YMN2) 
INTEGERS CALDATA(M*NTCH) 

DIMENSION VMN1 (M) * YMN2 ( M) 

CHARACTERS TEXT 
READ(9) TEXT 
PRINT 1001, TEXT 
1001 FORMAT (A ) 

READ(9) CALDATA 
REWIND 9 

PRINT *, 'DATA READING FINISHED* 

CLOSE (9) 

NSAMPL * NTCH/NCH 
J = 1 
K = NTCH 
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DO 44 II * I,M 

sum * o.o 

SUM2 * 0.0 

DO 33 I = J,K,NCH 

0ATA1 * FLCAT(CALDATA(I ))*5. 0/2043.0 
DATA2 * FLOAKCALDATAdn) )*5. 0/2048.0 
SUM1 = SUMI ♦ CATA1 

SUM2 * SUM2 + DATA2 

33 CONTINUE 

VMNl(II) * SUMI /FLO AT (N SAM PL ) 

VMN2UI) * SUM 2 /FLO AT (N SAUPL ) 

J * J ♦ NTCH 
44 K = K ♦ NTCH 

PRINT *, 'SAMPLES SEPARATEO ANO AVERAGING DONE' 

RETURN 

END 


SUBROUTINE CALBRT 

PROGRAM TO CALCULATE OPTIMUM VALUE OF CONSTANTS IN 
CALIBRATION EQUATION OF HOT-HIRE 
CAIBRATION EQUATIONS E +*2 * A ♦ B U**N 
INPUTS 
X * U 
FX * E 

M * NO OF DATA SETS 
OUTPUTS 

OPTIMUM VALUES OF A* 8» N 

ERMIN = MINIMUM ERROR IN CALCULATED VALUE OF X 

SUBROUTINE CALBRT (X ,FX .M.AOPT ,BOPT ,OPTN , ERMIN J 
DIMENSION X(M) »FX(M ) 

EXPNT = 0.3 
DEXP = 0.01 
J = 1 

HR I TE (8 »105 ) 

50 CALL LSTSQtX ,FX .EXPNT.M .A. B ) 

SUM = 0.0 
DO 1 I = 1 »M 

XC * ( (FXC I)**2 - A)/BM+(1. O/EXPNT ) 

FACT * (1.0 - XC/X(I))**2 
1 SUM * SUM ♦ FACT 

ERX * SCRT (SUM/FLOAT (M) ) 

IF (J .EC. 1) GO TO 10 
IF (ERX .LT. ERMIN) GO TO 10 
GO TO 30 
10 ERMIN = ERX 
A OPT = A 
BQPT = 8 
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OPTN « EXP NT 

33 IF (EXPNT .LE. 0.6) THEN 

105 FCRMAT (8X»'A'»10X»*B'»1QX« *N • ,7X . »ERR (U) • /) 
WRITE (8 *106) A.B.EXPNT.ZRX 

106 FORHA T(2X»3(3X»F8.5)»2X»E12.5) 

EXPNT » EXPNT ♦ DEXP 

J = J ♦ 1 
GC TO 50 
EN DIF 
RETURN 
ENC 


SUBROUTINE LSTSQ 

SUBROUTINE TO CALCULATE CONSTANTS (A ANC 6) IN 
EQUATION: FX**2 « A ♦ B X**EXPNT IN LEAST 

SQUARE SENSE FCR A GIVEN VALUE OF EXPNT 

SUBROUTINE LSTSQ (X »FX » EXPNT ,M»A ,B) 

CIMENSION X(H) ,FX(M) 

SMX * 0.0 
SMFX * 0.0 
SMXFX « 0.0 
SMX2 * 0.0 
DO 1 I * 1 »M 
XI » X ( I ) **EXPNT 
FX1 * FX ( I )+*2 
SMX * SMX ♦ XI 
SMFX « SMFX + FX1 
XFX • XI * FX1 
SMXFX » SMXFX ♦ XFX 
X2 « XI ♦ XI 
1 SMX2 • SMX2 ♦ X2 

DET « SMX2*FL0AT(M) - SMX*SMX 
A « CSMFX*SMX2 - SMX*SMXFX ) /CET 
B * CSMXFX+F LOAT (M) - SHFX*SMX)/DET 
RETURN 
END 
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APPENDIX B 


COMPUTER PROGRAMS FOR CALCULATING WALL PRESSURE- 
VELOCITY CORRELATION AND POWER SPECTRA OF FLUCTUATING 

WALL PRESSURE 


A brief description on the computer programs, written for the 
calculation of wall pressure-velocity correlation and the power spectra 
of the fluctuating wall pressure, and their FORTRAN listings are pre- 
sented in this appendix. 

B.l Computer Program for the Calculation of Wall Pressure-Velocity 

Correlation 

A computer program was written to calculate the wall pressure- 
velocity correlation incorporating the steps outlined in Chapter III 
(Section 3.3.4). The digitized data in the program was processed by 
dividing the complete record length of the signal in to number of 
blocks containing equal number of digitized samples. The various re- 
quired parameters were calculated first for each block and then aver- 
aged over the number of blocks selected. The temperature correction 
for the hot-wire voltage output was included in the computer program. 
In an absence of the fluctuating velocity data, a provision was made in 
the program to analyze only the fluctuating wall pressure data for the 
calculation of mean and RMS values of the wall pressure fluctuation. A 
FORTRAN listing (PROGRAM 1) of the computer program is included 
here. 
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PROGRAM 1 


INCLUDE • /USR/INCLUDE/HR.F* 

PROGRAM TO CALCULATE CORRELATION OF WALL PRESSURE 
FLUCTUATION WITH THE VELOCITY FLUCTUATIONS USING 
DIGITIZED DATA 

PROGRAM READS DIGITIZED DATA STORED ON DISK FILES UNDER 
NAME 'DAT . I JK * 

PROGRAM WRITTEN IN USER INTERACTIVE NODE 
NCH- TOTAL NO. OF DIGITIZED CHANNELS 
TNSAMPL- TOTAL NO. OF SAMPLES COLLECTED 
RFRQ- SAMPLING RATE 

BRFRQ- THROUGHPUT RATE OF A/C UNIT * 1MKZ 
MDP- NO. OF DATA SETS 
CALDATA- DIGITIZED SAMPLES 

INTEGERS CALDATA (lOOOOOO ) 

INTEGER NCK, TNSAMPL, MDP 
REAL RFRQ, BRFRQ 
CHARACTER*20 FLNH 

PRINT *, 'ENTER FILE NAME WHERE DATA ARE STORED* 

READ 101, FLNM 
101 FORMAT ( A J 

OPEN (9, FORM* 'UNFORMATTED', ST ATUS= 'OLD • ,F I LE=FLNM) 

REWIND 9 

READ (9) NCH, TNSAMPL, RFRQ, BRFRQ, MDP 
PRINT +, NCH, TNSAMPL, RFRQ, BRFRQ, MDP 
NSAMPL * TNSAMPL/NCH 

CALL OISCREADtCALDATA, NSAMPL , NCH , MDP , RFRQ ) 

STOP 

END 

SUBROUTINE DISCREAD 

SUBROUTINE DI SCRE AD (CALDATA .NSAMPL , NCH .MDP .RFRQ ) 

READS DIGITIZED DATA AND PROCESS IT TO OBTAIN 
CORRELATION OF HALL PRESSURE FLUCTUATION WITH 
THE VELOCITY FLUCTUATION 

CALCULATIONS DONE IN BLOCKS EACH HAVING 4096 
SAMPLES MAXIMUM 

A, B, AN- CALIBRATION CONSTANTS IN KING'S LAW 
CELT- TEMPERATURE DIFFERENCE IN DEG. F FOR 
CORRECTING HOT-WIRE OUTPUT 

CC * CALIBRATION CONST. FOR VARIABLE RELUCTANCE TYPE 
DIFFERENTIAL PRESSURE TRANSDUCER IN IN. /VOLT 
CF1 = CONVERSION FACTOR FOR 1 IN. OF HG TO N/M**2 

CF2 * CONVERSION FACTOR FOR 1 IN. OF WATER TO N/M**2 

CF3 = CONVERSION FACTOR FOR 1 PSI TO N/K**2 

RAIR * GAS CONSTANT FOR AIR IN J/KG.K 
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CH. 1 USED FOR REF. TOT. PRESSURE 

CH. 2 USED FOR HALL PRESSURE FLUCTUATION 

CH. 3 USED FOR VELOCITY FLUCTUATION 


INTECER*2 CALOATA (NCH.NSAMPL ) 

DIMENSION PI (4096), UI (4096 ) 

CHARACTER*80 TEXT 
CHARACTER*20 OUTFILE 
READ(9) TEXT 
PRINT 102. TEXT 
102 FORMAT C2X, A) 

READI9) CALDATA 
REWINO 9 

PRINT *, ’DATA READING FINISHED* 

CLOSE (9) 

PRINT *, 'ENTER FILE NAME FOR OUTPUT DATA AS DATOUT.IJK' 
READ 101, OUTFILE 
101 FORMAT (A) 

PRINT 102, OUTFILE 

OPEN (8, FORM* 'FORMATTED ' ,FILE*OUTFILE } 

REWIND 6 

WRITE (8,200) TEXT 

PRINT *, 'ENTER BLOCK SIZE OF SAMPLES- NBLKSZ ' 

READ *, NBLKSZ 

1 PRINT *, 'ENTER NO. OF BLOCK- NBLK ' 

READ *. NBLK 

MBLK * INT (NSAMPL/N8LKSZ) 

IFCNBLK .GT. MBLK) THEN 

PRINT *. 'REDUCE NO. OF BLOCKS < OR * TO', MBLK 

GO TO 1 

ENOIF 

PRINT *. 'ENTER ATM. PRESS. IN IN. OF HG = PATH* 

READ *, PATM 

PRINT *, 'ENTER ROOM TEMP. IN DEG. F * TATM* 

READ *, TATM 

PRINT +, 'ENTER CALIBRATION CONST. IN IN. /VOLT' 

READ *, CC 

PRINT *, 'ENTER AMPLIFICATION GAIN FACTOR' 

READ *. GF 

PRINT ♦, 'ENTER TRANSDUCER SENSITIVITY IN MY/PSI* 

READ *, S 
CF1 * 3386.6 
CF2 * 249.09 
RAIR * 287.1 
CF3 * 6E94 .76 

TABS * 5.0+tTATM - 32.0J/9.0 ♦ 273.15 
RHO * PATM+CF 1/ (RA IR+TA2S) 

WRITE (8 ,201) RFRQ.NSAMPL 
WRITE (8 ,202) NBLKSZ ,NBLK 
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WRITE (8. 203) PATH. TABS, CF.S.CC.RHO 
C BRANCHING 

PRINT *. 'ENTER I FCR RMS CR 2 FOR PU-CORREL AT ION* 
READ ♦, INUM 

IF (INUM .EQ. 1) GO TO 500 
PRINT *, 'ENTER A. B, N :• 

READ *» ' A « 8, AN 
PRINT • A » B. N:'» A, 8, AN 

WRITE (8*204 ) A, B. AN 
PRINT *, 'ENTER OELT ' 

READ ♦, DELT 
OELC * DELT*5. 0/9.0 
ALPHA * 0.0042 
R * 1.5 

FAC1 « ALPHA/I 2 .0* ( R - I.0>) 

FAC 2 * - FAC1*0ELC 

UREF * 0.0 

UM * 0.0 

UI2H « C.O 

PI2M « 0.0 

PM « 0.0 

PUIM * 0.0 

DO 11 II « l.NBLK 

SUREF * 0.0 

SPI « 0.0 

SPI2 * 0.0 

SUI * 0.0 

SUI2 * 0.0 

SPUI * 0.0 

CC 12 I * 1 »NBLKSZ 

K * I ♦ (II - 1 )*NBLKSZ 

VDAT * FL0AT(CALDATA(1,K))*5. 0/2048.0 

PCAT * FLOAT (C ALDATA (2 »K) ) 45 .0/2048 .0 

UCAT * FLOAT (C ALDATA (3«K))45 .0/2048.0 

UDAT » UDAT*(1 - FAC2) 

PI (I ) * PDAT*1000.0*CF3/(GF*S> 

UI(I) * ((UDAT4U0AT - A )/B > ♦♦ ( 1 .O/AN) 

L'I2 « UI(I)*UI(I) 

PI2 * PI(I )*PI (I) 

SPI « SPI ♦ PI C I } 

SUI « SUI ♦ UI(I) 

SUI2 « SUI2 ♦ UI2 

SPUI * SPUI 4 P I ( I ) 4UI (I) 

SPI2 * SPI 2 4 PI2 
12 SUREF * SUREF 4 YOAT 

SUREF » SUREF/FLOAT (NBLKSZ ) 

UREF * SQRT(2.0*SUREF4CC*CF2/RH0) ♦ UREF 
PM « SP 1/FLOAT (NBLKSZ) 4 PM 
PI2M = SPI2/FL0AT (NBLKSZ) ♦ PI2M 
UM * 3U I/FLOAT (NBLK SZ ) 4 UM 
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UI2M « SU 12/FLOAT (NBIKSZ) + UI2M 
PUIM * SPUI/FLOAT(NBLKSZ) ♦ PUIM 
PRINT + . * THE .II* *TH BLOCK CALCULATION OVER' 

11 CONTINUE 

UROF « UREF/FLOAT (NBLK) 

QINF * 0.5+RH0*UREF*UREF 
PM * PM/FLOAT (NBLK) 

PP2M * PI2M/FL0AT (NBLK ) - PM*PM 
UM * UM/FLQAT(NBLK } 

UP2M « UI2M/FL0AT (NBLK ) - UM*UM 
PUM « PUIM/FLOAT (NBLK) - PM*UM 
RMSP * SORT(PP2M J 
RATIO - RMSP/OINF 
RMSU « SQRKUP2M) 

TU * RMSU/UREF 

CORPU » PUM/CRMSP+RMSU) 

PRINT *, 'UREFIM/S) '.UREF. ' QINF(PA) * .QINF 
PRINT *» ' U*AVG (M/S ) * .UM» ' P-AVC (PA J • ,PM 
PRINT *, 'RMS-P(PA) '.RMSP.' RM S-U (M/S )' .RMSU 
PRINT +, 'PU-AVG'.PUM.' CORCOF-PU ' .CORPU 
WRITE (8*205) UREF ,Q INF , PM, RMSP. RAT 10 .UM, 

1 RMSU, TU, PUM, CORPU 
CLOSE (3) 

RETURN 

500 PRINT *, 'RMS VALUE CALCULATION OF WALL PRESSURE' 
UREF * 0.0 
PI2M * 0.0 
PM » 0.0 

DC 21 II * 1 .NBLK 

SUREF « 0.0 

SPI * 0.0 

SPI2 * 0.0 

DO 22 I * 1.NBLKS2 

K ■ I ♦ (II - U+NBLKSZ 

VOAT * FLOAT (C AL3ATA (1 ,K) ) *5.0/2048.0 

POAT * FLOAT (C AL3ATA ( 2 , K) } +5.0/2 04 8 .0 

PI(I) * PDAT*1000.0*CF3/(GF*S> 

SUREF * SUREF ♦ VDAT 
P 12 * PI(I)*PI(I) 

SPI - SPI ♦ PI(I) 

22 SPI2 * SPI 2 ♦ PI 2 

URCF * SUREF /FLOAT (NBLK SZ ) * UREF 
PM * SPI /FLOAT (NBLKSZ ) « PM 
PI2M * SPI 2/FLOAT (NBLKSZ) ♦ PI2M 
21 CONTINUE 

UREF * UREF/FLOAT (NBLK) 

IF (UREF .LE. 0.0) UREF * ABS(UREF) 

UREF * SQRT(2.0*UREF*CC+CF2/RH0) 

QINF * 0 • 5*RH0*UREF+UREF 
PM * PM/FLOAT (NBLK ) 
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PP2M * PI2M/FL0AT(NBIK) - PM*PM 
RMSP * SORT (PP2M) 

RATIO * RMSP/Q IMF 

PRINT*. 'UREF (M/S) • '.UREF,‘ Q INF ( PA 1 : • .0 INF 
PRINT ♦ , P-AVG(PA) : '.PM 

PRINT *» 'RMS-P(PA) :* ,RMSP.‘ RMS-P/Q INF :• .RATIO 
WRITE (3*207) UR EF»5INF 
WRITE (8 *206 ) PH »RHSP .RATIO 
2C6 FORMAT (5X.*P“MEAN(PA)®' »F10 »4//5X » 'P-RMS ( PA ) * * . 

1 F10.4,5X.‘P-RMS/QINF*‘,E15.7//) 

201 FCRMATI//5X. ‘SAMPLING RATE* • ,F8 . 1 ,10X , 

1 • SAMPLE S/CHANNEL* * » I6//J 

202 F0RMAT(5X . ’BLOCK SIZE* • .15 ,5X , ‘NO. OF BLOCKS* •, 13// ) 

203 FORMAT (/2X , ' ATM • PRESS. (IN. HO * • ,F7.3 ,4 X , 

1 ' ABS • TEMP. (OEG. K)*‘ . F10. 3//2X , »GAI N FACTOR**. FS.l , 

2 10X, ‘TRANSDUCER SENSITIV1 TY (MV/PS! ) * • .F7.3//2X . 

3 ‘DIFF. TRANSDUCER SENSIT I V ITY ( IN ./VOLT ) * ' ,F7.4// 

A 2X . ' A IR 0ENSITY(KG/M*43)*» ,F7.4/) 

2 04 FORMAT (/5X, ‘A* • ,F10 .4 ,5X, »B * • .F10 .4 ,5X. ‘N* • ,Fl0.4/> 
205 FORMAT (/5X. 'U-REF (M/S ) * ' ,F7 .4 ,11X, ‘Q I NF (PA ) = » .F8.4// 

1 5X.'P-MEAN(PA)*',F10.4,8X,‘P-RMS(PA)*‘ .F8.4//5X, 

2 'P-RM S/Q1NF* ‘.E15.7//5X, *U-ME AN( M/S )* * .F10.4.7X* 

3 ‘U-RMS(M/S)*» .FI0.4//5X, ‘TL=* .E12.5//5X. 

4 ‘PU-MEAN(N/M-S)*‘ .E15.7.5X, 

5 ‘CORR. COEFF** .E15.7//) 

200 FORMAT (5X .A) 

207 FORMAT (/5X»*U-REF(M/S)*',F6.4,5X.‘QINF(PA) =‘,F8.4/) 
CLOSE (8) 

RETURN 

END 
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BJ2 Computer Program for the Calculation of Power Spectra of 

Fluctuating Wall Pressure 

The procedure for calculating the power spectra of fluctuating 
wall pressure using Fast Fourier Transform (FFT) technique was pre- 
sented in Chapter III (Section 3.3.2). A computer program was written 
to calculate the power spectral density function of fluctuating wall 
pressure, at discrete frequencies, from its Fourier components ob- 
tained by a standard FFT subroutine. The calculations of power spec- 
tral density functions at discrete frequencies were obtained in blocks 
containing equal number of digitized samples. The averaging carried 
out over the number of blocks selected, then gave the final power 
spectral density function at discrete frequencies. The frequency band- 
width could be varied, if required, by changing the size of the block. 
The steps to calculate the intensity of wall pressure fluctuation, from 
the power spectral density distribution obtained, were also included in 
the program. A FORTRAN listing (PROGRAM 2) of the computer pro- 
gram is included here. 
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PROGRAM 2 


INCLUDE • /US F/ INCLUDE /MR. F * 

PROGRAM TO CALCULATE POWER SPECTRAL DENSITY FUNCTION 
CF WALL PRESSURE FLUCTUATIONS USING FAST FOURIER 
TRANSFORM TECHNIQUE 

PROGRAM REACS DIGITIZED DATA STORED ON DISC FILES 
UNDER NAME 'DAT . I JK • 

PROGRAM WRITTEN IN USER INTERACTIVE MODE 
NCH- TOTAL NO OF SIMULTANEOUSLY DIGITIZED CHANNELS 
TNSAMPL- TOTAL NO. OF SAMPLES COLLECTED 
RFRQ- SAMPLING RATE 

BRFRQ- THROUGHPUT RATE CF A/S UNIT = 1MHZ 
MDP- NO. OF DATA SETS 
RAWDATA- DIGITIZED SAMPLES 

NCN- CH. NO. CN WHICH FFT ANALYSIS TO BE DONE 

INTEGER*!! RAWDATAUOOOOO) 

INTEGER NCH, TNSAMPL .NCN .NU 
REAL RFRQ, BRFRQ 
CHARACTER*20 FLNM 

PRINT +, 'ENTER FILE NAME WHERE DATA ARE STORED* 

READ 1001, FLNM 
1001 FORMAT (A) 

PRINT *,’ ENTER CHANNEL NO. FOR FFT* 

READ *, NCN 

PRINT *, ’ENTER EXPONENT FOR LENGTH OF DATA ARRAY* 

READ ♦, NU 
NL * 2**NU 

PRINT ♦.•CHANNEL*' , NCN, 'LENGTH OF DATA ARRAY = • ,NL 
CPEN (3, FORM* 'UNFORMATTED* . STATUS* 'OLD* .FILE =FLNM) 
REWIND 8 

REA0(8) NCH, TNSAMPL. RFRC, BRFRQ, MDP 
PRINT *, NCH, TNSAMPL. RFRQ, BRFRQ, MDP 
NSAMPL = TNSAMPL/NCH 
NFFT * INT (NSAMPL/NL J 

111 PRINT *, 'ENTER NO. OF 2**NU 3LCCKS TO AVG. FFT* 

READ +, NF1 

IF INF 1 .GT. NFFTI THEN 

PRINT * , ’REDUCE NO. OF BLOCKS < OR = TO’, NFFT 

GO TO 111 

ENDIF 

PRINT ♦. ’ENTER GAIN FACTOR OF AMPLIFICATION* 

READ *, GF 

PRINT *, ’ENTER TRANSDUCER SENSITIVITY IN MV/PSI • 

READ *, SF 

CALL D I SC RE ADC RAWDATA, NSAMPL ,NCH ,NCN ,NU,RFRO ,NF 1 ,GF ,SF) 

STOP 

END 
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SUBROUTINE DISCREAO 

SUBROUTINE DISCREAD(RAWDATA .NSAMPL»NCH»NCN»NU»RFRQ,NF »GF »SF ) 

READS DIGITIZED DATA AND PROCESS IT TO OBTAIN 
POWER SPECTRAL DENSITY FUNCTION DISTRIBUTION 
N- LENGTH OF DATA ARRAY FOR FFT * 2+<‘NU 
NU- EXPONENT FOR FFT (SHOULD HOT 3E > 12 * » ALSO 
CONTROLS FREQUENCY 2ANDW1DTH AT A GIVEN 
SAMPLING RATE 

INTEGER+2 RAWDATA (NCH.NSAMPL) 

DIMENSION DAT ARY (4096 ) , X IMA G( 4096 ) .SPC(4096) 

DIMENSIGN SMDB ( 40S6) 

CHARACTER+30 TEXT 
CHARACTER*20 OUTFILE 
READ (3) TEXT 
PRINT 1002. TEXT 
1002 FORMAT (2X .A) 

READ (3) RAWDATA 
REWINO 8 

PRINT ♦.•DATA REAOING FINISHED* 

CLOSE (8) 

PRINT ♦, ’ENTER FILE NAME FOR OUTPUT AS DATOUT.IJK' 

READ 1301. OUTFILE 
1001 FORMAT (A) 

PRINT 1002. OUTFILE 

OPEN ( 3 .FORM* * FORMATTED • .FILE “OUTFILE ) 

REWIND 3 
WRITE (3*1) TEXT 

1 FORMAT (2X .A ) 

N = 2+*NU 

WRITE (8, 2) NCN.RFRQ.NSAMPL.NU.N.NF.GF.SF 

2 FORMAT (2X, 'CHANNEL NO. FOR FFT:»I2//2X. 

1 ’SAMPLING RATE *’,F10.1,5X. 

2 • SAMPLES/CHANNEL *'»I &//2X, 

2 ’EXPONENT FOR FFT *',I3//2X, 

3 ’LENGTH QF DATA ARRAY FOR FFT *»,I5//2X, 

A ’NO. OF BLOCKS TO AVERAGE FFT *»,I4//2X. 

5 ’AMPLIFICATION GAIN FACTOR »'»F7.1//2X, 

6 ’TRANSDUCER SENSITIVI TY(MV/PSI I *'.F7.3/J 

C C FI * CONVERSION FACTOR FOR 1 PSI TO PASCAL 

CF1 * 6894.76 
DO 133 I = l.N 
SFC(I) * 0.0 
133 SMDB(I) = 0.0 
NMl * N/2 ~ 1 
DC 34 IT * 1 .NF 
DO 33 I * l.N 
II * I * (IT - 1 J*N 
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DATA * FLOAT (RAWDAT A (NCN»I 1 ) )*5 . 0/2048.0 
DATARY ( I ) * DATA+1000 «0*CF 1/(GF+SF) 

XIMAG (I) =0.0 

33 CONTINUE 
SMSPC * 0.0 

PRINT ♦»»TiiE», IT.'TH FFT STARTS* 

CALL FFT (DAT ARY»XIMAG»N »NU ) 

REFDB * 20 .0*1 .0E-06 
DO 99 J = 1,N 

SP = IDATARY(J)**2 ♦ XIMAG ( J)**2 )*2/ (RFRQ*FLOAT IN) ) 
DB = 20.0* ALOGlOtSQRT(SP) /REFDB) 

SPC(J> * SP ♦ SPC(J) 

SMDB(J) = OB ♦ SMDBt J) 

FREQ = FLOAT ( J )*RFRQ/FLCAT ( N ) 

99 CONTINUE 

DO 199 JJ * 2 ,NMl 

SMSPC * SMSPC * (SPC(JJ) ♦ SPC(JJ ♦ 1))/IT 
199 CONTINUE 

PSQ = 0 .5*SMSPC+RFRQ/FLCAT (N) 

RP = SQRTCPSQ) 

PRINT *, *P-RMS * »RP 

34 CONTINUE 
WRITE (3»3 J 

3 FORMAT {/5X , * I • , 4X , *F (HZ ) * . 7X , * SPC { PA**2/HZ ) * , 

1 5X » * SPC CDB ) • / ) 

M * N/2 
DC 35 I = l.M 

FREQ * FLOAT (I )*RFRQ/FLCAT (N) 

SPC(I) * SPC ( I ) /FLOAT (NF) 

SMOB(I) = SMDB ( I J/FLCAT INF ) 

35 WRI TE (8 *4 ) I .FREQ , SPC < I ) .SHOE ( I) 

4 FORMAT (2X, 14 ,2X ,F10 .3 » 3X»E 1 3 .5 »5X , F7 .2/) 

WRITE (8»5) RP 

5 FORMAT (5X»*P-RMS(PA) = *,F10.4> 

RETURN 

END 

SUBROUTINE FFT 

SUBROUTINE FFT (XREAL .XIMAG ,N,NU) 

XREAL- REAL PART OF THE INSTANTANEOUS SIGNAL 
XIMAG- IMAGINARY PART OF THE INSTANTANEOUS SIGNAL 

C IMENSION XREAL(N).XIHAGtN) 

N2 = N/2 
NU1 = NU - 1 
K = 0 

DO 100 L * 1 »NU 
102 DO 101 1=1 *N2 
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o n o 


P * I BITR (K/2**NU1 iNU) 

ARG = 6.283185»P/FL0AT(NI 
C * CQS(ARG) 

S = SIN(ARG) 

K1 * K ♦ 1 
K1N2 * K1 ♦ N2 

TREAL * XREAL(K1N2]*C ♦ XINAG(K1N2)*S 
TIMAG = XIMAG(K1N2)*C - XREAL(K1N2 )+S 
XREAL (K1N2 } = XREAL(Kl) - TREAL 
X IMAG (K1N2 ) * XIMAG(Kl) - TIMAG 
XREAL(Kl) * XREAl(Kl) ♦ TREAL 
XIMAGIK1) * XIMAG(Kl) + TIMAG 
101 K « K ♦ 1 
K * K ♦ N2 

I F ( K .LT. N) GO TO 102 
K = 0 

NL'l = NU1 - 1 
100 N2 * N2/2 

00 103 K = 1 »N 

1 = IBITR(K - 1 »NU ) + 1 
IFl I .LE. K) GO TO 103 
TREAL * XREAL(K) 

TIMAG * XIMAG(K) 

XREAL (K) = XREAL(I) 

XIMAG(K) * XIMAG(I) 

XREAL ( I ) * TREAL 
XIMAG(I) * TIMAG 
103 CONTINUE 
RETURN 
END 

FUNCTION IB1TR 

FUNCTION IBITR(J.NU) 

J1 = J 
IBITR = 0 
CO 2C0 I * l.NU 
J2 = Jl/2 

IBITR * IB ITR+2 + (J1 - 2*J2) 

200 J1 * J2 
RETURN 
END 
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APPENDIX C 


YAW SENSITIVITIES OF THE SURFACE STATIC AND THE 

PRESTON TUBE 


The yaw characteristics of the static tube resting on the surface 
(termed as the surface static tube) and the Preston tube were experi- 
mentally investigated in two- and three-dimensional flows. The results 
obtained on the yaw characteristics for both the tubes are presented 
and discussed here. 

C.1 Yaw Sensitivity of the Surface Static Tube 

Fig. C.l shows the variation of (P s (a) - P s (ao)], normalized by the 

free stream dynamic pressure, with the yaw angle for three different 
flow conditions, where P s (a) is the static pressure read by the surface 
static tube at an arbitrary yaw angle a and P s («o) is the static pressure 

recorded by the surface static tube when aligned in the local flow di- 
rection ao. The value of the angle ocq corresponds to the maximum 

output recorded by the surface static tube. It is quite evident, from the 
different flow conditions investigated, that the static pressure ob- 
tained with the surface static tube is less if it is not aligned in the local 
flow direction. Also, the static tube resting on the surface, similar to 
the static tube away from the solid wall, is insensitive to the yaw angle 
of ± 5°. However, the error in measured static pressure increases with 
the increase in the yaw angle. It is also noted that the surface static 
tube is comparatively less sensitive to the yaw angle in the presence of. 
adverse pressure gradient. 
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Different Flow Conditions 



In a three-dimensional flow with high degree of skewness, the 
error in measured static pressure using the surface static tube due to 
its misalignment with the local flow direction is shown in Fig. C.2. It is 
clearly observed that the smaller value of the static pressure will be 
obtained if the static tube is not aligned in the local flow direction. 
Moreover, the magnitude of the error in static pressure due to mis- 
alignment is small (< 1%) if the error is expressed as a fraction of free 
stream dynamic head. It is important to mention here that, in general, 
the normalization of static pressure with the free stream dynamic pre- 
ssure is used to evaluate error in the calculated velocity. In Fig. C.2(b) 
the error in measured static pressure due to misalignment is norm- 
alized with the static pressure recorded by the surface static tube 
when it is aligned in the local flow direction. It is quite evident, that if 
the static tube were not aligned in the local flow direction it would 
have recorded approximately 20% smaller static pressure. 

CJ2 Taw Sensitivity of the Preston Tube 

Fig. C.3(a) shows the variation of [Pp(a) - Pp(ao)]. normalized by 
Pp(cco), with the yaw angle for three different flow conditions, where 
Pp(a) is the total pressure read by the Preston tube at an angle a with 
the local flow direction and Pp(ao) is the total pressure recorded by 

the Preston tube when it is aligned in the local flow direction. The 
Preston tube is observed to be insensitive, the error being less than 
2%, to the yaw angle of ± 5®. However, the pressure recorded by the 
Preston tube reduces rapidly beyond the yaw angle of 5^. 

Fig. C.3(b) shows error in the total pressure recorded by the 
Preston tube, in a three-dimensional flow with the high degree of 
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Fig. C.2 Error in Surface Static Tube Pressure Due to its Misalignment 
(Note: a = 0 Implies that the Centerlines of the Surface Static 
Tube and the Flat Plate are Parallel) 
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(a) Yaw Sensitivity of Preston Tube in Different Flow Conditions 
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(b) Error in Preston Tube Pressure Due to its Misalignment With 
the Local Flow Direction- A Typical Error Distribution in 
Three-Dimensional Flow 

(Note: a = 0 Implies that the Centerlines of the Preston Tube 
and the Flat Plate are Parallel) 

Fig. C.3 Yaw Sensitivity of the Preston Tube 
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skewness, due to its misalignment. The value of a = 0 corresponds to 
the normal setting, i.e., parallel to the undisturbed flow direction, of 
the Preston tube. It is observed, that if the Preston tube were not 
aligned in the local flow direction, it would have recorded approxi- 
mately 15% less total pressure. For the wall shear stress measure- 
ment, this observed error in the Preston tube output due to misalign- 
ment would cause approximately 10% error in the calculated value of 
the wall shear stress. 
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